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A Elisa, Irene y Carmita 




En la presente obra se proponen una serie de problemas ajustados a la Termodinámica 
de Gibbs que suele corresponder a la segunda etapa que se imparte en los cursos de 
Termodinámica en las Licenciaturas de Física, Química o Ingeniería. 

Se proponen un conjunto de problemas resueltos adaptados a la metodología particu¬ 
lar de la obra. Su resolución se ha planteado de manera que permita al alumno no sólo 
resolver los problemas, sino comprender los conceptos utilizados así como analizar cuá¬ 
les son los objetivos que se pretenden cubrir con su resolución dentro de cada tema. 

Con ello se facilitan herramientas que permiten al alumno adquirir una mejor for¬ 
mación en la solución de cualquier problema en general, así como tener una visión 
más concreta de algunos sistemas y procesos particulares aplicando el formalismo 
general de la Termodinámica. 

Los dos primeros capítulos incorporan un resumen de los conceptos teóricos bási¬ 
cos utilizados a lo largo de la obra, con objeto de centrar la atención en los aspectos 
más relevantes del desarrollo de la Termodinámica, lo cual sirve además como repa¬ 
so de los mismos. 

Asimismo se incluyen comentarios, anotaciones históricas o puntos de vista sobre 
problemas particulares que pretenden llamar la atención del lector y hacer que la lectura 
de la obra resulte más amenta y agradable. 
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Capítulo 1 

Termodinámica de los 
sistemas cerrados 


Una de las formas habituales de desarrollar la Termodinámica del Equilibrio, comúnmente 
expresada como Termodinámica, es seguir un esquema conceptual bastante próximo a 
lo que fué su desarrollo histórico, cimentado a mitad del siglo XIX en los trabajos de 
James Joule, William Thomson (Lord Kelvin) y Rudolf Clausius, con las aportaciones 
posteriores de Josiah W. Gibbs. Este esquema general, sistematizado por Max Planck 
a finales de dicho siglo, es lo que se conoce como Fundamentación Clausius-Kelvin de 
la Termodinámica del equilibrio de los sistemas cerrados. 

Si bien se supone que el lector ya ha seguido en un curso teórico este desarrollo 
y está familiarizado con el mismo, resumiremos a continuación sus principales defi¬ 
niciones y resultados con el objeto de centrar las ideas que van a ser utilizadas a 
continuación a lo largo del presente libro. 

• Sistema Termodinámico. 

— Un sistema termodinámico es todo sistema macroscópico limitado por 
una superficie cerrada real (pared) o ficticia. 

• Contactos Termodinámicos. 

— Cuando la variación de las propiedades físico-químicas de un sistema Si 
puede inducir variaciones en las de otro X 2 se dice que ambos entran en 

interacción. 

— Las interacciones entre sistemas termodinámicos reciben también el nombre 

de contactos termodinámicos. 

— Un sistema S se dice que está aislado, cuando no puede existir ningún 
contacto entre él y su medio exterior. 
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• Estados de equilibrio. 


— Se denomina estado de un sistema al conjunto de los valores que to¬ 
man, en un instante dado, los parámetros macroscópicos que describen sus 
propiedades físico-químicas. 

— Experimentalmente se observa que en todo sistema termodinámico existen 
unos estados particulares denominados estados de equilibrio, tales que 
no sólo son estables en el tiempo, sino que lo continúan siendo cuando se 
aísla el sistema o cualquier subsistema del mismo. 

• Ligaduras. 

— Se denominan ligaduras de un sistema, a las condiciones restrictivas a las 
que éste se encuentra sometido y que restringen a determinados intervalos 
o a un valor concreto los parámetros que lo describen. 

• Procesos Termodinámicos. 

— Un proceso termodinámico es el paso del sistema desde un estado de 
equilibrio hasta otro estado también de equilibrio. 

— Un proceso espontáneo es el proceso inducido en el sistema por la elimi¬ 
nación de una o varias ligaduras internas. 

— Un proceso en el que las variables del sistema cambian en una cantidad 
infinitesimal se denomina proceso infinitesimal. 

— Se denomina sucesión de procesos infinitesimales a un proceso finito que 
tiene lugar mediante sucesivas etapas infinitesimales. 

— Un proceso se dice que es cíclico, cuando el estado final coincide con el 
inicial. 

• Trabajo Termodinámico. 

— Se denomina trabajo realizado sobre un sistema termodinámico o trabajo 
termodinámico, en un proceso dado, al trabajo realizado por las fuerzas 
que el medio exterior ejerce sobre el sistema. 

— El análisis de los distintos sistemas termodinámicos ha permitido consta¬ 
tar la existencia, para cada sistema E, de un conjunto de variables termo¬ 
dinámicas extensivas {X,}™ =1 , denominadas coordenadas de trabajo, y 
un conjunto de variables intensivas {Y¿}™ =1 , denominadas variables con¬ 
jugadas tales que la variación dXi de cada una de las variables extensivas 
en un proceso infinitesimal cualquiera, implica la realización de un trabajo 
sobre E, dado por 

n 

sw = J2 Y ^x t ( 1 . 1 ) 

i= 1 

— Trabajo disipativo Wdi s es el trabajo realizado sobre el sistema no aso¬ 
ciado a la variación de las coordenadas de trabajo. 
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— En un proceso infinitesimal en general se tiene 


sw = Y, YidXi + 6Wdi ° ( 1 . 2 ) 

2=1 


• Proceso cuasiestático. 

— Se denomina proceso infinitesimal cuasiestático a un proceso infinitesimal 
en el que no se realiza trabajo disipativo, y proceso cuasiestático a toda 
sucesión de procesos infinitesimales cuasiestáticos. 

• Postulado de las variables de estado. 

— Se denominan variables de estado de un sistema a un conjunto cualquie¬ 
ra de variables independientes que permite especificar completamente sus 
estados de equilibrio. 

— Se constata empíricamente, que los estados de equilibrio de un sistema ho¬ 
mogéneo cerrado con n coordenadas de trabajo, quedan completamente es¬ 
pecificados por n - 1-1 variables independientes: ai, a 2 ,...,a ra +i- 

• Fuente de trabajo. 

— Un sistema se dice que es una fuente de trabajo para otro, respecto a 
una determinada coordenada de trabajo, cuando en todo contacto mecánico 
entre ambos relativo a dicha coordenada, la correspondiente variable con¬ 
jugada de la fuente permanece constante. 

— Por lo tanto, en todos los estados de equilibrio mecánico entre un sistema 
E y una fuente de trabajo Ey-., la variable Y¿ conjugada de X t satisfacerá 
la condición 

Yi = YP (1.3) 

donde Yf ; es el correspondiente valor de la variable conjugada para la 
fuente. 

• Contacto térmico. 

— Se dice que dos sistemas están en contacto material cuando pueden in¬ 
tercambiar materia entre si. 

— Se dice que dos sistemas están en contacto mecánico cuando pueden 
realizar trabajo entre ellos. 

— Se constata experimentalmente, que aún cuando un sistema E esté sometido 
a ligaduras que impiden los contactos material y mecánico con el medio 
exterior es todavía posible inducir, desde el exterior, procesos en él. 

— Las interacciones de este tipo reciben el nombre de interacciones o contac¬ 
tos térmicos. 
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• Pared adiabática. 


— Se constata experimentalmente la existencia de unas paredes impermeables 
que impiden el contacto térmico, denominadas paredes adiabáticas. 

— Todo proceso que tiene lugar en un sistema rodeado de una pared adiabática 
se denomina proceso adiabático. 

• Equilibrio térmico. 

— Se observa que al poner en contacto dos sistemas aislados Ej y £2 mediante 
la eliminación de determinadas ligaduras, permaneciendo el sistema total 
E — Ei + £2 aislado, los sistemas £ a , a = 1,2 abandonan su estado inicial 
para alcanzar cada uno de ellos un nuevo estado de equilibrio, diciéndose 
entonces que ambos sistemas están en equilibrio mutuo. 

— Se dice, asimismo que estos estados finales son estados de equilibrio mutuo 
respecto de los contactos existentes entre y 

— En particular, dos sistemas están en equilibrio térmico cuando se en¬ 
cuentran en equilibrio mutuo respecto del contacto térmico. 

• Principio cero. 

— Dos sistemas cualesquiera en equilibrio térmico con un tercero se encuen¬ 
tran a su vez en equilibrio térmico entre si. 

• Isotermas. 

— Se denomina isoterma al conjunto de estados de un sistema en equilibrio 
térmico con uno dado 

— Los estados de una isoterma constituyen una curva del espacio termo- 
dinámico del sistema. 

— Se denominan isotermas correspondientes de dos sistemas a las curvas 
en las que se sitúan los estados de equilibrio térmico de uno y otro sistema. 

• Temperatura empírica. 

— El principio cero permite definir en cada sistema termodinámico una función 
de estado 9 (ai,..., a n +i) que se denomina temperatura empírica, la 
cual permanece constante en una isoterma 

0(ai,...,a„+i) = 9 0 (1.4) 

— Definir una escala de temperaturas es adoptar un convenio en el que 
la temperatura empírica adopte el mismo valor numérico en las isotermas 
correspondientes de los distintos sistemas. 
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• Primer principio. 

— A lo largo del desarrollo de la Termodinámica se ha aceptado que ésta puede 
cimentarse a partir de la aceptación sin demostración de unos principios o 
leyes generales enunciando evidencias no contradichas por la naturaleza. 
Estos principios establecidos por Clausius han sufrido algunas variantes 
o matizaciones con objeto de aportar un mayor grado de rigurosidad al 
esquema conceptual. Así pues el enunciado que aquí presentamos para el 
primer principio es planteado por Born y Caratheodory en los primeros 
años del siglo XX. 

— Enunciado: 

* a) Dados dos estados de equilibrio cualesquiera Bi y f ?2 de un sistema 
cerrado, siempre es posible alcanzar mediante un proceso adiabático, 
bien el estado B 2 a partir del B\ o bien el B\ a partir del B 2 . 

* b) El trabajo realizado sobre el sistema en un proceso adiabático sólo 
depende de los estados inicial y final. 

• Energía interna. 

— A partir de este enunciado del primer principio, se introduce el concepto de 
energía interna U, definida como una variable termodinámica tal que la 
diferencia de energía interna entre dos estados B\ y B 2 del sistema es 

U 2 — U\ = W$ (1.5) 

(o U 2 — U\ = —W 21 ) 

donde Wf 2 (o W 2 { 1 ) representa el trabajo realizado sobre el sistema en 
cualquier proceso adiabático que lo lleve del estado B 1 al B 2 (o del B 2 al 
Bi). 

— Con esta definición, la energía interna es una función potencial para el 
trabajo realizado sobre un sistema en un proceso adiabático. 

• Calor. 

— Una vez definida la energía interna, se introduce el concepto de calor ab¬ 
sorbido QI 2 por el sistema, en un proceso c del estado B 1 al , como 

Qí 2 = U 2 -U 1 - W¡2 ( 1 . 6 ) 

siendo Wf 2 el trabajo realizado sobre el sistema en dicho proceso y U 2 — U\ 
la correspondiente variación de la energía interna. 

— Según la definición de calor absorbido por el sistema, en un proceso adiabático 
el calor absorbido es nulo. 

— En un proceso infinitesimal 


dU = SQ + óW 


(1.7) 
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• Fuente de calor. 


— Un sistema Ex¡, es una fuente de calor para otro sistema E, cuando en 
las interacciones térmicas entre E y Ej¡, la temperatura de Ejj, permanece 
constante. 

— De este modo en los estados de equilibrio mutuo entre Ex¡, y E la tempe¬ 
ratura de este último será siempre igual a To. 

• Proceso reversible. 


— Se dice que un sistema cerrado E experimenta un proceso reversible si, 
una vez finalizado el proceso, es posible devolver a E a su estado inicial 
de forma que los alrededores E a ; recuperen asimismo su estado inicial. Es 
decir, si es posible completar un proceso cíclico en E de manera que tam¬ 
bién los alrededores, y por tanto el sistema total E+E tt /, experimenten un 
proceso cíclico. Cuando esto no es posible se dice que se trata de un proceso 
irreversible. 


• Motor térmico y máquina frigorífica. 


— Se concibe un motor térmico como un sistema de evolución cíclica cuyo 
objetivo es realizar trabajo sobre los alrededores a costa de absorber una 
cierta cantidad de calor. 

— Se define el rendimiento 77 de un motor térmico que absorbe un calor Q 2 
de una fuente E t 2 a temperatura T 2 , cede un calor Q\ a una fuente Er, a 
temperatura Ti, con T 2 > Ti, y realiza sobre los alrededores E a i un trabajo 
W', como el cociente entre el trabajo producido y el calor absorbido 


W_ 

Q2 


( 1 . 8 ) 


— Se denomina máquina frigorífica a un sistema de evolución cíclica cuyo 
objetivo es extraer calor de un sistema más frío que sus alrededores a costa 
de recibir una cierta cantidad de trabajo. 


— Se define la eficiencia £ de una máquina frigorífica que absorbe en cada 
ciclo un calor Qi de una fuente fría a temperatura Ti y cede un calor Q' 2 
a una fuente caliente a temperatura T 2 , T 2 > Ti, recibiendo un trabajo W, 
como el cociente entre el calor que absorbe de la fuente fría y el trabajo 
realizado sobre la máquina 



(1.9) 
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• Segundo principio. 


— Se proponen dos enunciados equivalentes: 

* Enunciado de Clausius 

— No existe ningún sistema termodinámica cuyo único efecto al evolucionar 
cíclicamente sea extraer calor de un sistema y cederlo a otro que se encuen¬ 
tra a mayor temperatura que el anterior. 

* Enunciado de Kelvin Planck 

— No existe ningún sistema termodinámica cuyo único efecto, al evolucionar 
cíclicamente, sea extraer calor de un determinado sistema y realizar la mis¬ 
ma cantidad de trabajo. 

• Teorema de Carnot 

A partir del segundo principio se demuestra el Teorema de Carnot, que junto 
con el correspondiente corolario permite establecer que: 

— El rendimiento de un motor térmico reversible que interacciona con sólo 
dos fuentes de calor es independiente de la naturaleza del sistema que actúa 
como motor y sólo depende de la temperatura de ambas fuentes, siendo este 
rendimiento mayor que el de cualquier otro motor que opere irreversible¬ 
mente entre ambas fuentes de calor. 

• Temperatura termodinámica. 

— El teorema de Carnot permite introducir la Temperatura Termodinámica 
como 



( 1 . 10 ) 


donde Q y Q' s son, respectivamente, los calores absorbido y cedido, por 
el sistema que recorre un ciclo de Carnot, de una fuente de calor Ey a la 
temperatura T y una fuente de calor de referencia E t s a la temperatura 
Ts, cuyo valor se fija por convenio. 

• Teorema de Clausius. 

— Del segundo principio se sigue, asimismo, el Teorema de Clausius: 

- Cuando un sistema cerrado experimenta un proceso cíclico en el que entra 
en contacto con un conjunto de fuentes de calor a las temperaturas termo¬ 
dinámicas Ti, la suma algebraica de los cocientes entre los calores absorbidos 
de cada fuente y las temperaturas de las mismas es menor o igual que cero, 
verificándose la igualdad si y sólo si el proceso es reversible. 



( 1 . 11 ) 
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• Entropía. 


— El teorema de Clausius permite introducir la entropía de un sistema, defi¬ 
nida como una variable tal que la diferencia de entropía entre dos estados 
del sistema B i y B 2 es 


5 2 -Si 



( 1 . 12 ) 


siendo c un proceso reversible cualquiera entre los estados Bi y B 2 , SQ el 
calor absorbido por el sistema en cada etapa infinitesimal del mismo y T 
su temperatura en esa etapa. 


— En un proceso reversible infinitesimal 


6Q = TdS 


(1.13) 


— La definición de la entropía permite establecer que: 

* a) La variación de entropía en un proceso adiabático reversible de B\ 
a B 2 es cero. 

S 2 = S 1 (1.14) 

* b) La variación de entropía en un proceso adiabático irreversible de B\ 
a B 2 es positiva. 

S 2 > Si (1.15) 


• Principio de entropía máxima. 

— Se constata experimentalmente que los estados de equilibrio de un sistema 
aislado están determinados por el principio de entropía máxima: 

— Un estado de equilibrio de un sistema aislado es un estado de máxima en¬ 
tropía entre todos los posibles estados compatibles con las ligaduras del sis¬ 
tema que tengan su misma energía interna. 


Este principio puede enunciarse de forma equivalente como un principio de energía 
interna mínima. 

• Principio de energía interna mínima. 

— Un estado de equilibrio de un sistema termodinámico aislado es un estado 
de mínima energía interna respecto de todos los estados compatibles con las 
ligaduras del sistema y que tengan su misma entropía. 
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Capítulo 2 

Termodinámica de los 
sistemas de composición 
variable. 


• Variables de estado de un sistema abierto. 

— Es posible definir variaciones de energía interna y de entropía para 
un sistema que intercambia materia con los alrededores (sistema abierto) 
aunque los conceptos de calor y trabajo no estén perfectamente definidos 
en este caso 

— Los estados de equilibrio de un sistema homogéneo de composición varia¬ 
ble (bien sea por ser abierto o porque en el mismo puede tener lugar una 
reacción química) están unívocamente determinados por los valores de su 
energía interna, coordenadas de trabajo, y número de moles de las espe¬ 
cies químicas presentes en el sistema. El mismo resultado puede aplicarse 
tomando la entropía en lugar de la energía interna. 

— En el caso de que tenga lugar una reacción química, los números de moles 
de las especies químicas no son variables independientes. Sus variaciones 
están sujetas a unas proporciones determinadas por los coeficientes este- 
quiométricos Vj de cada especie química y pueden expresarse en términos 
de las variaciones de una única magnitud, el grado de avance £ 

dNj = Vjd(( (2-1) 

— La entropía de un sistema heterogéneo es la suma de las entropías de cada 
uno de los subsistemas que lo componen. La misma propiedad se verifica 
para la energía interna. 

— Para un sistema heterogéneo es posible definir un conjunto de variables 
independientes que definen su estado. Estas variables pueden elegirse entre 
las variables que definen cada uno de los subsistemas homogéneos, con las 
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relaciones que establecen las condiciones de ligadura a que están sometidos 
cada uno de los subsistemas. 

• Expresión diferencial de la energía interna de un sistema homogéneo 
abierto. Potencial químico. 

— La expresión 


dU = TdS - PdV +...+/i c dN c (2.2) 


expresa la variación infinitesimal de la energía interna en términos de las 
variaciones infinitesimales de entropía, volumen (o coordenadas de trabajo) 
y número de moles de las especies químicas constituyentes del sistema. La 
variable 




( dU 
\dNj 


(2.3) 


sv 


que representa la variación relativa de la energía interna frente al número 
de moles de la sustancia j-ésima es una magnitud intensiva que se denomina 
potencial químico de dicho componente. 


• Ecuaciones de Euler y de Gibbs-Duhem. 


— La energía interna es una magnitud extensiva y por tanto es una función 
homogénea de primer orden con respecto a las variables extensivas inde¬ 
pendientes del sistema S , V, { Nj }. Por tanto verifica la relación de Euler 
para las funciones homogéneas 


U = TS - PV + (2-4) 

3 = 1 

— La relación de Euler establece que existe una relación entre las variacio¬ 
nes infinitesimales de las magnitudes extensivas denominada ecuación de 
Gibbs Duhem 

-SdT + VdP - Ndn = 0 (2.5) 


• Condiciones de equilibrio 

— El principio de entropía máxima permite establecer condiciones de equi¬ 
librio entre los subsistemas de un sistema heterogéno. Asimismo permite 
establecer una condición de equilibrio para un sistema en el que puede tener 
lugar una reacción química. 

* En el caso de dos fases en equilibrio (sistemas separados por paredes 
diatermas y permeables), estas condiciones son 

• Igualdad de temperaturas 

• Igualdad de presiones 

• Igualdad de potenciales químicos en ambas fases para todos los 
componentes. 
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* En el caso de equilibrio osmótico (sistemas separados por paredes 
diatermas y semipermeables) 

• Igualdad de temperaturas 

■ Igualdad de potenciales químicos para las sustancias capaces de 
atravesar la membrana 

* En el caso de un sistema en el que pueda tener lugar una reacción 
química (equilibrio químico) 

■ La suma de los potenciales químicos ponderada por los coeficientes 
estequiométricos (con su signo correspondiente) debe ser igual a 
cero. 

• Regla de las fases. 

— Para un sistema con c componentes en el que coexisten f fases el número 
máximo de variables intensivas que pueden ser variadas independientemente 
(grados de libertad) es 

l = c + 2 — f (2.6) 

— El número máximo de fases que pueden coexistir en equilibrio para un sis¬ 
tema con c componentes es 

f<c + 2 (2.7) 


• Ecuación de Clausius-Clapeyron. 


— Para una sustancia pura la condición de equilibrio entre fases permite es¬ 
tablecer que 

* El equilibrio de fases tiene lugar a una sola presión para cada tempe¬ 
ratura (o viceversa) 

* La variación de presión y temperatura manteniendo el equilibrio de 
fases está ligada por la ecuación (Ecuación de Clausius-Clapeyron) 


dP s^l — s*- 1 ) 
dT = u( 2 ) - lá 1 ) 


( 2 . 8 ) 


donde s^, s^ son las entropías molares y v'- 1 ', los volúmenes mo¬ 
lares de las respectivas fases. 


• Calor latente. 


— Se define calor latente como el calor absorbido cuando un mol de sustancia 
pasa reversiblemente de una fase (1) a otra (2). Este calor latente está dado 
por 


¿12 = T (s (2) - s (1) ) 


(2.9) 
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• Punto triple 

— Para una sustancia pura en un diagrama P-T, la coexistencia entre dos 
fases está representada por una curva, y la coexistencia de tres fases por un 
punto, denominado punto triple. 

• Punto crítico. Ley de los estados correpondientes. 

— La curva de equilibrio líquido vapor termina en un punto de temperatura 
T c y presión P c característica de cada sustancia denominado punto crítico 
por encima del cual (es decir a T > T c o P > P c ) no es posible observar la 
transición de líquido a vapor, sino que un aumento de temperatura o presión 
implica una variación continua de las características físicas del fluido. 

— Muchas propiedades de los fluidos toman valores similares cuando sus varia¬ 
bles se normalizan relativas a los valores en el punto crítico. Estas variables 
se denominan coordenadas reducidas y a esta circunstancia ley de los es¬ 
tados correspondientes. 

• Ecuación fundamental. 

— Cualquier propiedad termodinámica puede ser expresada como una relación 
entre las 2c+5 variables (en general 2n+2c+3 si el sistema tiene n coorden- 
das de trabajo ) U, S , V, Ni, ..., N c , T, P, fii, ..., ¡i c o expresiones que puedan 
ser derivadas de las mismas. De estas variables c + 2 son independientes, 
por lo que existe un máximo de 2c 4- 3 relaciones independientes entre ellas. 

— De la ecuación 

U = U{S,V,N 1 ,...,N C ) (2.10) 

se pueden deducir el conjunto máximo de c + 3 relaciones independien¬ 
tes por lo que de la misma pueden deducirse todas las propiedades termo¬ 
dinámicas del sistema, y por lo tanto se le denomina ecuación fundamen¬ 
tal. (Ecuación fundamental en la representación de la energía o energética) 

— La ecuación 

S = S(U,V,Ni,...,N c ) (2.11) 

es equivalente a la anterior y por tanto también es una ecuación fundamen¬ 
tal. (Ecuación fundamental en la representación de la entropía o entrópica) 

— Si introducimos la variable F = U — TS (Energía libre o función de 
Helmholtz) y consideramos la relación F = F(T,V,Ni,...,N c ), o bien 
la variable G = U — TS + PV (Energía libre o función de Gibbs) y 
consideremos la relación G = G(T, P, Ni,..., N c ), ambas relaciones tienen 
también el carácter de ecuación fundamental del sistema. Asimismo la 
variable H = U+PV (Entalpia) cuando se considera H = (S, P, Ni, ..., N c ) 
es también una ecuación fundamental. 

— Cualquier transformada de Legendre de la ecuación fundamental energética 
o entrópica expresada en términos de sus variables canónicas o naturales es 
asimismo una ecuación fundamental. 
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• Potenciales termodinámicos 


— La función de Helmholtz cuya expresión diferencial se escribe como 

C 

dF = -SdT - PdV + Y N dN j (2.12) 

3 = 1 

en un potencial para el trabajo que es posible extraer de un sistema en 
un proceso reversible isotérmico. Este es el máximo trabajo que puede 
extraerse del sistema en condiciones isotérmicas. 

— La entalpia, cuya expresión diferencial se escribe como 

C 

dH = TdS + VdP + Y N dN (2.13) 

3=1 

es un potencial para el calor absorbido por un sistema cerrado en un proceso 
reversible isobárico. 

— La función de Gibbs cuya expresión diferencial está dada por 

C 

dG = —SdT + VdP + Y /V A j (2.14) 

3=1 

es un potencial para el trabajo químico realizado sobre un sistema cerrado 
en un proceso reversible a temperatura y presión constante. 

• Relaciones de Maxwell. 

— La aplicación del teorema de Schwartz de igualdad de las derivadas segun¬ 
das cruzadas de una función diferenciable a los potenciales termodinámicos 
conduce a un conjunto de relaciones de gran utilidad entre las derivadas 
de las variables termodinámicas denominadas relaciones de Maxwell. 
Algunos ejemplos son 


d 2 u \ _ 
dVdS , 


a 2 u 

OSdV 


(dT \ __ _ (dP\ 

V dV ) S,N \ OS ) V,N 


d ¿ F 
dVOT , 


ít_f 

dTOV 


(dS_\ 

V dV ) T,N 


(9P) 

\dT ) V,N 


d 2 H 
dSOP , 


d 2 H 

dPdS 


(dV\ 

V OS ) P.N 


(dT\ 

V dP ) S,N 


d 2 G 

dTOP 


d 2 G 
dPOT , 


(dV\ =-(9S) 

\dTjp,N \dPJT,N 
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• Relaciones de Gibbs-Helmholtz 


— Las definiciones de los potenciales termodinámicos conducen a algunas re¬ 
laciones de utilidad entre ellos, así por ejemplo las relaciones de Gibbs 
Helmholtz entre entalpia y función de Gibbs. 

H = G - T {§)„ (215) 

a = -r 2 (2i6) 


• Equilibrio estable. 

— Un sistema homogéneo a una determinada temperatura y presión se dice 
que se encuentra en un estado de equilibrio estable cuando dicho estado es 
un mínimo absoluto de la función de Gibbs entre todos los estados de igual 
temperatura y presión. Si existe un sistema en un estado que no verifica 
estas características se dice que se encuentra en un estado metaestable. 

• Mezcla de gases ideales. 

— En una mezcla de gases ideales el volumen de la mezcla es igual a la suma 
de los volúmenes de cada uno de los componentes a la misma temperatura 
y presión que la mezcla 

C 

V(T, P, {IV,}) = £ Vf(T, P , Nj) (2.17) 

3 = 1 


— Se define la presión parcial de un componente en una mezcla de gases 
ideales como la presión a la que estaría sometido dicho gas si él sólo ocupase 
el mismo volumen que la mezcla en las mismas condiciones de temperatura. 


NjRT 

V 


(2.18) 


— La presión parcial está relacionada con la presión total de la forma 


Pj = XjP ( 2 - 19 ) 

donde xj es la fracción molar de dicho componente en la mezcla. 

— La presión de una mezcla de gases ideales es la suma de las presiones par¬ 
ciales de cada uno de los componentes (Ley de Dalton). 

C 

P = J2 P 3 ( 2 - 20 ) 

3 = 1 
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La energía interna de una mezcla de gases ideales es la suma de las 
energías internas de cada componente a la misma presión y temperatura 
que la mezcla 

C 

U(T, P , {Nj}) = U°(T, P , Nj) (2.21) 

3 = 1 

Como la energía interna molar de un gas ideal sólo depende de la tempera¬ 
tura, la energía interna molar de una mezcla de gases ideales sólo depende 
de la temperatura. 

La entropía de una mezcla de gases ideales es la suma de las entropías 
de cada uno de los componentes a la misma temperatura y volumen que la 
mezcla 

C 

S(T, V, {Nj}) = S°( T > V, Nj) ( 2 . 22 ) 

3 = 1 

esta igualdad no se verifica cuando la entropía es expresada en términos de 
temperatura y presión. 

La Función de Helmholtz de una mezcla de gases ideales es la suma de 
las Funciones de Helmholtz de cada uno de los componentes a la misma 
temperatura y volumen que la mezcla 

C 

F(T, V,, {Nj}) = V.] Nj) (2.23) 

3 = 1 

por tanto es posible en esta representación obtener una ecuación fundamen¬ 
tal de la mezcla de gases ideales como suma de la ecuación fundamental de 
cada uno de ellos. 

El potencial químico de un componente en una mezcla de gases ideales 
es igual a su potencial químico a la misma temperatura y volumen que la 
mezcla 

l ij (T,V,{N j })=n°j(T,V) (2.24) 

El potencial químico de un componente en una mezcla de gases ideales a 
una determinada temperatura y presión es igual a su potencial químico a 
la presión parcial y misma temperatura que la mezcla. 

(T, P, {Nj}) = fj,j (T, Pj) (2.25) 

El potencial químico de un componente en una mezcla de gases ideales a 
una determinada temperatura y presión es igual a su potencial químico en 
estado puro a la misma temperatura y presión que la mezcla más el término 
RT ln Xj 

IXj (T, P , {Nj}) = n°(T, P) + RTlnxj (2.26) 

donde Xj es su fracción molar en la mezcla. 
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— Como la función de Gibbs verifica 


G(T, P, {Nj }) - £ N {T, P, {Nj}) (2.27) 

3 = 1 

La función de Gibbs de una mezcla de gases ideales es la suma de las 
funciones de Gibbs de cada uno de los componentes en estado puro a la 
presión parcial de cada componente en la mezcla 

C 

G(T, P, {Nj}) = G°(T, Pj , Nj) (2.28) 

3 = 1 


• Disolución ideal. 

— Se define una disolución como un sistema homogéneo pluricornponente. Se 
define disolución perfecta como aquella en la cual se verifica la condición 

fjtj(T, P , {N^) = Hj(T, P ) + RT\nxj (2.29) 

en todo el rango de concentraciones. 

— Una disolución se dice diluida cuando existe una sustancia que denomi¬ 
namos disolvente para la cual su fracción molar es próxima a la unidad 
X\ ss 1 mientras que para las restantes que denominamos solutos x s ~ 0. 
Se define una disolución ideal cuando en un determinado rango de con¬ 
centraciones el potencial químico depende linealmente con el logaritmo de 
la concentración, siendo la pendiente de dicha dependencia el término RT. 
Es decir 

fij(T, P , {JV ¿ }) = ii°/{T, P) + RT ln Xj (2.30) 

En el caso del disolvente es posible tomar lim^-u y por tanto /.x" (T, P) 
corresponde a su potencial químico en estado puro, y por tanto 

Nj}) - fi° 1 {T,P) + RT\nx 1 (2.31) 

• Disolución no ideal. Actividad y fugacidad. 

— En una disolución no ideal se define la actividad a 3 del componente j-ésimo 
de manera que se verifica 

Hj (T, P, {Nj }) = P ) + RT ln a 3 (2.32) 

— En una mezcla de gases no ideales se define la fugacidad de un componente 
en la mezcla f 7 de manera que 


n(T, p, {Nj}) = iT/T, P = 1) + RT ln fj (2.33) 
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Equilibrio de un líquido con una mezcla de gases poco solubles. 


— La presión parcial de un vapor en equilibrio con su líquido en una mezcla 
de gases ideales poco solubles en el mismo es igual a la presión de equilibrio 
líquido vapor de dicho componente en estado puro a la misma temperatura. 

Humedad relativa. 


— Denominamos presión de saturación como la presión de equilibrio líquido 
vapor para una sustancia pura a una determinada temperatura. 

— Se define la humedad relativa como el cociente (en tanto por ciento) 
entre la presión parcial del vapor de agua en la atmósfera y la presión de 
saturación a la misma temperatura. 

H r = 100§ (2.34) 

S 

Disolución de una mezcla de gases en un líquido. Ley de Henry. 

— A una temperatura dada la concentración de un gas que se disuelve en un 
liquido es proporcional a la presión parcial de dicho componente en la fase 
gaseosa (Ley de Henry) 

Disolución ideal en equilibrio con su vapor. Ley de Raoult. 

— La presión parcial de un componente volátil en equilibrio con una mezcla de 
gases es proporcional a su fracción molar en la fase líquida, siendo la cons¬ 
tante de proporcionalidad igual a su presión de saturación (Ley de Raoult) 

Pj = P°Xj (2.35) 

Descenso de la presión de vapor. 


— La reducción relativa de la presión de vapor de un disolvente cuando en el 
mismo se disuelve una determinada sustancia es proporcional a la concen¬ 
tración de ésta. 



(2.36) 


Equilibrio líquido vapor para un sistema multicomponente 


— Para un sistema formado por dos o más componentes la temperatura de 
ebullición a una presión dada no es única sino que depende de las concen¬ 
traciones relativas de los distintos compuestos. 

— Cuando para un sistema multicomponente existe equilibrio entre la fase 
líquida y la gaseosa, esta última es más rica en el componente más volátil 
(de menor punto de ebullición). 

— Existen mezclas en las cuales el equilibrio líquido-gas es posible para las 
mismas concentraciones en el gas y en el líquido. Estas mezclas se denomi¬ 
nan azeotrópicas. 



• Aumento ebulloscópico. 


— El aumento de la temperatura de ebullición de una disolución ideal 
es proporcional a la concentración de las sustancias disueltas en la misma. 
La constante de proporcionalidad depende de la temperatura de ebullición 
del disolvente (al cuadrado), la constante de los gases y el calor latente del 
disolvente. 



(2.37) 


• Descenso crioscópico 

— El descenso del punto de congelación de una disolución ideal es propor¬ 
cional a la concentración de las sustancias disueltas en ella. La constante 
de proporcionalidad depende de la temperatura de congelación del disol¬ 
vente (al cuadrado), la constante de los gases y el calor latente de fusión 
del disolvente. 



(2.38) 


• Ley del reparto de Nernst. 

— Cuando un soluto se reparte entre dos fases líquidas no miscibles la rela¬ 
ción entre las concentraciones está en relación inversa de las correspondien¬ 
tes constantes de Henry relativas a ambas disoluciones.(Ley del reparto de 
Nernst) 


• Presión osmótica. 

— La presión osmótica de una disolución es proporcional a la concentración 
de las sustancias disueltas. La constante de proporcionalidad depende de 
la temperatura, la constante de los gases y el volumen molar del disolvente. 



(2.40) 
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• Calor de dilución. 


— Se denomina calor de dilución o calor de disolución a la diferencia entre 
la entalpia de la disolución y la entalpia de los componentes en estado puro 
a una determinada presión y temperatura. 

C 

Q sol = H(T, P, {Nj}) - ]T H°(T, P , Nj) (2.41) 

3 =1 

• Equilibrio químico para una reacción a T y P constante. 


— En un sistema homogéneo que se mantiene a temperatura y presión cons¬ 
tante en el que tiene lugar una reacción química el estado de equilibrio se 
alcanza cuando la función de Gibbs es mínima. En el caso en que dicho 
equilibrio se alcance con una coexistencia entre reactivos y productos éste 
tiene lugar cuando 

C 

E^ = 0 ( 2 - 42 ) 

3 = 1 


• Afinidad. 


— Se define la afinidad A como 

“-(IL 


(2.43) 


de manera que dG = — Ad £ en un desplazamiento infinitesimal de la reac¬ 
ción 

* La afinidad también suele denotarse A = — A G y al término A G se le 
denomina energía libre de reacción. 

— La afinidad verifica 

C 

A=(-AG) = -5>^ (2.44) 

3 = 1 

Se define la afinidad propia A° como 

c 

A° = (-AG°) = -5>^ (2.45) 

3 = 1 

donde ¡j,? es el potencial químico de la sustancia j-ésima en estado puro. 

• Reacción química entre fases puras. 


— En una reacción química entre fases puras la función de Gibbs es una función 
lineal del grado de avance 


G 



(2.46) 
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por lo tanto el equilibrio se produce en uno de los extremos del intervalo 
£ € [0,1] siendo la reacción completa (£ = 1) cuando (Xq=i < 0 o 

imposible (£ = 0) cuando (Xq=i > 0 

— En términos de la afinidad propia una reacción entre fases puras será com¬ 
pleta cuando ésta sea positiva A° > 0 e imposible cuando es negativa 
A° < 0. 

• Reacción química en una disolución ideal 

— En una reacción química en una disolución ideal en el equilibrio se verifica 

C 

]Jx?=K(T,P) (2.47) 

3 =1 

donde la función K(T, P) denominada constante de equilibrio está dada 
por 

K(T, P ) = e A °/ RT (2.48) 

siendo A° la afinidad propia. 

— En una reacción química en una disolución si 

* A° » 0 entonces el equilibrio se produce de manera que la fracción 
molar de los productos es grande frente a la de los reactivos. 

* A° « 0 entonces el equilibrio se produce de manera que la fracción 
molar de los productos es pequeña frente a la de los reactivos. 

• Entalpia de reacción. Calor de reacción. 

— Se denomina entalpia propia de reacción o calor de reacción a la diferencia 
de entalpia entre el estado correspondiente a los productos y a los reactivos 
ambos en estado puro 

C 

AH° = i; 3 h" (2.49) 

3=0 

— Atendiendo al signo del calor de reacción, para una reacción que tiene lugar 
a temperatura y presión constante 

* Si A H° < 0 la reacción es exotérmica, es decir libera calor en el 
transcurso de la misma. 

* Si A H° > 0 la reacción es endotérmica, es decir absorbe calor en el 
transcurso de la misma. 

• Entropía de reacción. 

— Se define la entropía propia de reacción AS° como la diferencia entre las 
entropías molares de los reactivos y productos en estado puro ponderadas 
con los respectivos coeficientes estequiométricos. 

C 

AS° = 5>*¿ (2-50) 

3 = 1 


30 


• Cálculo de las afinidades propias. 


— Las afinidades propias A° = —A G° pueden calcularse a partir de los calores 
de reacción y las entropías propias de reacción 


(2.51) 


A G° = AH° - TAS° 


• Ley de Hess. 

— El calor de reacción de una reacción combinación lineal de otras es asimismo 
combinación lineal de los calores de reacción de éstas (Ley de Hess) 

— Toda reacción puede ponerse como combinación lineal de las reacciones de 
formación de cada uno de los compuestos que intervienen en la misma, 
por tanto el calor de reacción puede ponerse como combinación lineal de 
los calores de formación de dichos compuestos. Esta misma propiedad se 
verifica para la entropía de reacción y la afinidad propia. 

• Variación del calor de reación con la temperatura. 

— La variación de la entalpia propia de reacción con la temperatura es igual 
a la diferencia de las capacidades caloríficas de reactivos y productos pon¬ 
derados con los correspondientes coeficientes estequiométricos. 



(2.52) 


• Variación de la afinidad propia con la temperatura. 

— La relación de Gibbs Hclmholtz permite relacionar la variación de la afini¬ 
dad propia con la temperatura con el calor de reacción 



(2.53) 


— Como A° = RT ln K(T, P) la variación de la constante de equilibrio con la 
temperatura está dada por 



(2.54) 


• Variación de la afinidad propia con la presión. 

— La variación de la afinidad propia con la presión está dada por 



(2.55) 


— En una reacción entre gases ideales 



(2.56) 
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y por lo tanto si 

* ^ 2j—i v j^j <0 la reacción se ve favorecida por un aumento de presión. 

* (e;- -,) < 0 la reacción se ve desfavorecida por un aumento de 
presión. 

• Principio de Le Chatelier. 

— Un aumento de temperatura provoca un desplazamiento hacia la derecha 
(mayor formación de productos) en el caso de una reacción endotérmica y 
hacia la izquierda (menor formación de productos) en una reacción exotérmica. 

— Un aumento de presión provoca un desplazamiento de la reacción en el 
sentido en el que disminuye el volumen total del sistema. 

— Por tanto podemos decir que: 

En el caso de que en un sistema se introduzca una perturbación exterior 
el sistema reacciona en el sentido que tiende a anular dicha perturbación 
(Principio de Le Chatelier). 

• Estado estándar 

— Se considera estado estándar de una sustancia aquel de temperatura T = 
298.15 K (25°C) y P = 1 bar. Las magnitudes referidas a este estado se 
denominan magnitudes estándar y suelen encontrarse tabuladas para las 
reacciones de formación de las distintas sustancias. 

• Reacción en una pila electroquímica. 

— En una reacción que tiene lugar en un pila electroquímica la variación 
de la energía libre o función de Gibbs cuando 1 mol de sustancia es inter¬ 
cambiado en uno de los electrodos está dada por 


(2.57) 


tx 2 — C\ — — n e J~S 


donde n e es el número de electrones de valencia del elemento que se in¬ 
tercambia en la reacción, T es la constante de Faraday y ¿i es la fuerza 
electromotriz de la pila. 

— En el mismo proceso 



(2.58) 


— y por tanto 



(2.59) 
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• Enunciado de Thompsen y Berthelot. 

— Thompsen y Berthelot plantearon erróneamente el enunciado 

* Cualquier reacción química que tiene lugar sin la realización de trabajo 
externo tiene lugar en el sentido en que produce la máxima cantidad de 
calor, o sea una reacción es espontánea sólo si es exotérmica 

AH < 0 (2.60) 


• Enunciado de Gibbs. 

— El planteamiento de Gibbs permite enunciar que 

* Cualquier reacción química en condiciones de temperatura y presión 
constantes evoluciona en el sentido en que hace mínima la función de 
Gibbs. 

AG<0 (2.61) 

— Como A G = AH — TAS, el principio de Thompsen y Berthelot es equi¬ 
valente al de Gibbs en aquellos casos en los que AS > 0 o bien AS < 0 
pero TAS es pequeño frente a AH lo cual ocurre en un gran número de 
reacciones. 

• Teorema del calor de Nernst. 


— Nernst postuló en su teorema del calor que en las proximidades del cero 
absoluto 

dAG 8AH 

lim —r=- = lim ^ = 0 (2.62) 

t^o dT T — >o dT y ’ 

• Tercer principio de la Termodinámica. 


— Planck enunció este resultado de la forma 

* En el cero absoluto un proceso isotérmico es asimismo isoentrópico. 

— El Tercer principio de la Termodinámica establece que 

* En toda sustancia pura en equilibrio termodinámica estable y tal que si 
se encuentra en fase sólida forma un cristal perfecto la entropía puede 
tomarse como nula en el cero absoluto 


lim S = 0 (2.63) 

T —>0 

— O bien en forma de un principio de inaccesibilidad 

* Es imposible reducir la temperatura de un sistema al cero absoluto me¬ 
diante un número finito de operaciones. 
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• Contrastaciones experimentales del tercer principio. 

— El coeficiente de dilatación isobárica a = y (§§^) P se anula en el cero ab¬ 
soluto. 

— La pendiente de la curva de equilibrio de fases se anula en el cero absoluto 
— Los calores específicos se anulan en el cero absoluto. 

• Mecanismos de transmisión del calor. 

— Un cuerpo puede intercambiar calor con los alrededores a través de los 
siguientes mecanismos: 

* Conducción. Contacto directo con otro cuerpo sin que tenga lugar 
movimiento de materia en el entorno. También puede darse entre dis¬ 
tintas regiones de un mismo sistema 

* Convección. Contacto con un fluido en movimiento que contribuye 
de forma conjunta al transporte de masa y energía. 

* Radiación. Emisión (/absorción) de radiación electromagnética sin 
contacto material directo con los alrededores. 

• Ley de Fourier 

— El flujo de calor a través de una placa de un sólido homogéneo e isótropo 
de espesor Aa: y área A, cuyas paredes se mantienen a las temperaturas Ti 
y To respectivamente está dado por 

$ = (2.64) 

donde k es una constante característica del material denominada conduc¬ 
tividad térmica. 

• Ecuación de Fourier 


— El flujo de calor que atraviesa una superficie orientada que rodea una por¬ 
ción de volumen V de un sólido homogéneo e isótropo está dado por 

<I> = - í kVT-dA (2.65) 

Ja 

donde VT es el vector gradiente de temperatura VT = y 

gL 4 el vector elemento de superficie (vector de módulo la superficie gL 4 y 
orientación normal al elemento de superficie). 

— Este flujo de calor puede ponerse utilizando el teorema de la divergencia 
como 

i> = - f kV 2 TdV (2.66) 

Jv 
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— El calor perdido por unidad de tiempo por un volumen V de un sólido 
homogéneo e isótropo como consecuencia de la variación de temperatura en 
su interior está dado por 



(2.67) 


— La ecuación de distribución de temperaturas en un sólido homogéneo e 
isótropo (Ecuación de Fourier) está dada por 



( 2 . 68 ) 


• Transmisión de calor por radiación. 

— Se llama emitancia 77 a la energía radiada por un cuerpo por unidad de 
tiempo y superficie. 

— Se denomina irradiación Tí a la energía de la radiación que incide sobre 
un cuerpo por unidad de tiempo y superficie. 

— Se define la absortividad como la fracción de la energía que incide sobre 
un cuerpo que es absorbida por éste. 

• Cuerpo negro. 

— Se define cuerpo negro como una sustancia con absortividad igual a la 
unidad. 

• Ley de Kirchoff. 

— La irradiación dentro de una cavidad de paredes de temperatura uniforme 
una vez alcanzado el equilibrio es independiente de la naturaleza de los 
materiales que la forman y es únicamente función de la temperatura. 

— La irradiación en una cavidad de paredes de temperatura uniforme T es 
igual a la emitancia de un cuerpo negro a la misma temperatura. 

— La emitancia de un cuerpo negro es únicamente función de la temperatura. 

— La emitancia de un cuerpo cualquiera a temperatura T es igual al producto 
de su absortividad por la emitancia de un cuerpo negro a la misma tempe¬ 
ratura (Ley de Kirchoff). 

• Ley de Stefan Boltzmann 

— El calor absorbido por radiación por un cuerpo a temperatura T cuyos alre¬ 
dedores se encuentran a una temperatura T p es proporcional a la diferencia 
de las potencias cuartas de las temperaturas de los alrededores y del cuerpo. 


SQ oc (T^ - T 4 ) 


(2.69) 
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— La emitancia de un cuerpo negro es proporcional a la potencia cuart.a de su 
temperatura (Ley de Stefan Boltzmann) 

1Z b {T) = oT 4 (2.70) 

donde a la constante a se le denomina constante de Stefan Boltzmann. 

• Transmisión de calor por convección. 

— El flujo de calor entre una superficie sólida a temperatura T s y un fluido 
a temperatura T a es proporcional a la diferencia de temperaturas entre el 
sólido y el fluido 

$ = hA(T s - T a ) (2.71) 

donde A es la superficie y h es un parámetro denominado coeficiente de 
convección que depende de las características del sólido y del fluido, así 
como del estado de movimiento de éste último. 

• Radiación electromagnética como sistema termodinámico. 

— La radiación electromagnética encerrada en una cavidad de volumen 
V de paredes de temperatura uniforme T puede considerarse como un 
sistema termodinámico en equilibrio. 

— Su volumen está definido como el volumen de la cavidad, su temperatura 
como la de las paredes de ésta y su energía interna como la variación de 
energía interna que experimentaría un cuerpo que la absorbiera completa¬ 
mente. 

— La radiación electromagnética que incide sobre la superficie de un cuerpo 
ejerce presión sobre ésta y está determinada por un tercio de la densidad 
de radiación (u v = y) 

P = ^u v (2.72) 

— La densidad de radiación es proporcional a la potencia cuarta de la tempe¬ 
ratura 

U = aT^V (2.73) 

— La relación de Maxwell (§y) T = (fr)v P erm A e ver que I a entropía de la 
radiación en equilibrio es proporcional al cubo de la temperatura 

S(T,V) = ^aT 3 V (2.74) 

O 

— La ecuación fundamental de la radiación electromagnética en equilibrio con 
las paredes de una cavidad está dada por 

S(U,V) = f-aiu^vi (2.75) 

O 
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• Sistemas magnéticos. 

— Se define el vector de inducción magnética B de forma que la fuerza que 
actúa sobre una carga q que se mueve con velocidad v a lo largo de un 
campo magnético está dada por 

F = qv x B (2.76) 

— Se define la magnetización o momento magnético por unidad de volumen 
como la diferencia entre la inducción magnética en un punto de la muestra y 
la inducción magnética en el vacío (en ausencia de muestra) en dicho punto. 

M = — (É-B 0 ) (2.77) 

fJ-o V 7 

donde /z 0 es una constante denominada permeabilidad del vacío. 

— Se denomina intensidad del campo aplicado al término 

H = —B 0 (2.78) 

Ho 

— Se define la susceptibilidad magnética n m de manera que 

M = n m H (2.79) 

• Propiedades magnéticas de los materiales. 

— Los materiales con susceptibilidad magnética negativa se denominan dia¬ 
magnéticos. En ellos el campo magnético externo se atenúa en el interior 
del material. 

— En las sustancias con susceptibilidad magnética positiva pude ocurrir 

* La susceptibilidad no varía apreciablemente con el campo ni con la 
temperatura y la sustancia se desmagnetiza al desaparecer el campo 
exterior: Materiales paramagnéticos 

* La susceptibilidad no varía apreciablemente con el campo pero aumenta 
con la temperatura hasta una cierta temperatura de transición: Mate¬ 
riales antiferromagnéticos. 

* La susceptibilidad aumenta significativamente con el campo y la sustan¬ 
cia permanece magnetizada cuando el campo desaparece. Materiales 
Ferromagnéticos. 

• Trabajo realizado sobre un sistema magnético. 

— Para un material homogéneo en el que el campo es constante en el interior 
del mismo se define la magnetización total como 


M = VM 


( 2 . 80 ) 


donde M es el módulo de la magnetización. 
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— El trabajo realizado sobre un sistema homogéno en el que el campo es 
constante en su interior para variar su magnetización en d.M está dado por 

bW = B a dM (2.81) 

— Un sistema magnético puede considerarse como un sistema termodinámico 
con la magnetización total Ai como coordenada de trabajo y la inducción 
mágnética en el vacío B a como variable conjugada. 

• Superconductividad. 

— Kammerling Onnes observó en 1911 que la resistencia de algunos materiales 
que denominó superconductores se anulaba por debajo de una cierta 
temperatura. 

— Meissner probó en 1933 que un estado de superconductor es un estado de 
diamagnetismo perfecto es decir que el campo en el interior de un super¬ 
conductor es nulo 

5 (s) = 0 (2.82) 

• Ecuación de conservación de la energía para un fluido en régimen 
estacionario. 

— La ecuación de la conservación de la energía para un fluido que circu¬ 
la en régimen estacionario por un determinado recinto intercambiando 
calor y trabajo con los dispositivos del mismo se escribe como 

q- w ' = h 2 -h 1 + i(c¡ - c?) + (z 2 - Zi) g (2.83) 

donde q y w' representan el calor absorbido y el trabajo realizado por unidad 
de masa de fluido que circula, h \, h 2 son las respectivas entalpias específicas 
a la entrada y salida del recinto, C 2 ,ci las respectivas velocidades y z 2 ,z\ 
las alturas en el campo gravitatorio terrestre. 

• Dispositivos de intercambio de energía para un fluido que circula en 
régimen estacionario. 

— En una caldera o un intercambiador de calor un fluido absorbe una 
cantidad de calor a presión constante siendo los cambios de energía cinética 
y potencial despreciables así como el trabajo realizado por el fluido 

q = h 2 -h 1 (2.84) 

— En un compresor (si el fluido es un líquido suele llamarse bomba) el fluido 
entra a una determinada presión P\ y se comprime hasta una presión P 2 
después de haberse realizado un trabajo w sobre él, siendo despreciables la 
transferencia de calor y las variaciones de energía cinética y potencial 

w = h 2 — h\ (2.85) 
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— En una turbina se produce el efecto inverso que en el compresor 


w' = — (J 12 — h\) (2.86) 

— En un difusor un aumento de la sección del dispositivo por el que circula 
el fluido provoca una disminución de su velocidad, el intercambio de calor 
y la realización de trabajo se consideran despreciables así como la variación 
de energía potencial. 

h -2 - hi + ^(c¡ - c\) = 0 (2.87) 

— En una tobera se produce un efecto similar al del difusor sólo que en este 
caso un estrechamiento de la sección provoca un aumento de la velocidad 
del mismo. 

— En una válvula de estrangulamiento el fluido pasa de una región a alta 
presión a una región a baja presión sin realizar trabajo ni intercambiar 
calor con los alrededores siendo su variación de energía cinética y potencial 
despreciable (Efecto Joule-Kelvin) 

h 2 - h x = 0 (2.88) 



William Thomson “Lord Kelvin” (1824-1904) a la derecha de John William Strutt “Lord Rayleigh” (1842-1919). 
Ambos se destacaron por carácter polifacético realizando aportaciones en múltiples ramas de la Física y del cono¬ 
cimiento científico. 


39 





Josiah Willard Gibbs (1839-1903), con sus trabajos “Métodos gráficos en la Termodinámica de los fluidos”, “Un 
método de representación geométrica de las propiedades Termodinámicas de las sustancias por medio de superfi¬ 
cies” y “Sobre el equilibrio de las sustancias heterogéneas” publicados entre 1873 y 1875, sentó las bases del méto¬ 
do utilizado actualmente en la Termodinámica de los sistemas de composición variable, basado en la ecuación fun¬ 
damental y los principios de extremo. Su trabajo se difundió en Europa principalmente a través de la obra de James 
Clerk Maxwell (1831-1879) y fue fundamental para el amplio desarrollo de la fisicoquímica a finales del siglo XIX, 
principalmente a través de los trabajos de Vant’Hoff, Ostwald y Nemst entre otros. 
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Capítulo 3 

Variables de estado y 
ecuaciones de Euler y Gibbs 
Duhem 

3.1 Objetivos. 

• Identificar las variables de estado de un sistema abierto. 

• Comprender la descripción de un sistema en términos de: variables que lo des¬ 
criben, ligaduras a las que está sometido y variables independientes. 

• Familiarizarse con el uso de las ecuaciones de Euler y Gibbs Duhem. 

3.2 Problemas 


3-1) Describir las variables de estado de: a) El agua del grifo, b) Una 
botella cerrada de agua con gas, c) Una botella abierta de agua con 
gas, d) Una gota de agua, e) Una goma elástica. 


• a) El agua es un fluido y por tanto puede considerarse como un sistema con una 
única coordenada de trabajo, el volumen, por lo que su estado estará descrito por 
los valores de su energía interna, volumen y número de moles. Alternativamente 
podemos usar la entropía como variable independiente. 

— Las variables intensivas correspondientes serían temperatura, presión, y 
potencial químico. 
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— Sin embargo un vaso de agua recién extraído del grifo, debemos entenderlo 
como una masa de agua en equilibrio con la atmósfera, por lo que aparte 
del agua debemos contemplar como parte de su composición los gases at¬ 
mosféricos disueltos en la misma, sin perjuicio de considerar también otros 
componentes como sales minerales y sustancias utilizadas en la cloración. 

— Las condiciones de ligadura a que está sometido este sistema serían las de 
temperatura y presión constante igual a la ambiente. 

— El sistema es un sistema abierto y por tanto la cantidad de cada sustancia 
puede variar, en particular si la humedad de la atmósfera es pequeña el 
agua se evaporará después de un cierto tiempo. 

• b) Una botella cerrada de agua con gas es un sistema formado por dos fases 
homogéneas líquido y vapor. 

— Cada fase corresponde a un sistema expansivo por lo que aparece descrita 
por los valores de su energía interna, volumen y número de moles de cada 
sustancia presente. También podemos tomar la entropía en lugar de la 
energía interna como variable independiente. 

— Lo que se entiende habitualmente por agua con gas es una disolución de 
dióxido de carbono en agua por lo tanto en la fase de vapor habrá dos 
componentes CO 2 , H 2 O al igual que en la fase líquida, si bien en ésta 
un pequeño porcentaje forma un componente denominado ácido carbónico 
H 2 CO 3 por lo que en el número de moles debemos tener en cuenta Nco 2 > 
Nh 2 o, Nh 2 co 3 , sujetos a la ligadura de la existencia de la reacción química 

co 2 + h 2 o —► h 2 co 3 

y a que la masa total de los componentes es constante. 

— Las variables intensivas serán temperatura, presión y potencial químico de 
cada componente. 

— El volumen de cada subsistema es una variable ligada ya que el volumen 
total es constante + V^ = V < - bot ' ) 

— Si la botella no está aislada térmicamente la temperatura será una variable 
ligada y su valor corresponderá al de la temperatura atmosférica. 

• c)Una botella abierta de agua con gas, en el equilibrio obedece a la misma 
descripción del agua del grifo, (exceptuando las sales disueltas y los compuestos 
de cloración), es decir es una disolución de gases atmosféricos en agua. En 
el equilibrio, la proporción de CO 2 disuelto es la misma que la del agua del 
grifo (agua sin gas) sólo que necesitamos un tiempo muy largo para que dicho 
equilibrio se produzca. 

• d) Para describir una gota de agua consideremos en primer lugar una gota de 
agua en equilibrio con su vapor en un recipiente cerrado. 
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— El trabajo realizado sobre la gota está caracterizado por el producto de 
la presión (cambiada de signo) por la variación de volumen, (— PdV ) sin 
embargo cuando la gota es pequeña los efectos asociados a la tensión super¬ 
ficial son relevantes por lo que debemos incluir un segundo término (crdA) 
caracterizado por el producto del coeficiente de tensión superficial a por la 
superficie A, por tanto tenemos un sistema con dos coordenadas de traba¬ 
jo: volumen V y superficie A, y como variables conjugadas presión (-P) y 
coeficiente de tensión superficial a. 

— Por tanto los estados del sistema estarán descritos por la energía interna (o 
la entropía) el volumen, la superficie y el número de moles de agua presente 
en la gota. 

— Sin embargo, hemos de tener en cuenta que si la geometría de la gota es 
prefijada por ejemplo una gota esférica, el volumen y la superficie no son 
variables independientes ya que V = ^ 7 r?’ 3 ; A = 4nr 2 . Como dV = 4irr 2 dr 
y dA = 8 nrdr tenemos que dA = ^dV 

— Por tanto si tomamos como variables S , V,A,N y escribimos la expresión 
diferencial de la energía interna 


(3.1) 


dU = TdS - PdV + adA + p dN 


teniendo en cuenta la relación entre volumen y superficie la ecuación ante¬ 
rior puede escribirse como 



(3.2) 


donde ahora entropía, volumen y número de moles son variables indepen¬ 
dientes 

— En el caso de una gota de agua en la atmósfera es necesario tener en cuenta 
la presencia de los gases atmosféricos en la misma. Si bien esta presencia 
puede ser despreciable a muchos efectos puede ser un factor relevante por 
ejemplo cuando hay presente sustancias contaminantes lo cual es importante 
a la hora de explicar fenómenos como el smog producido en las ciudades 
con un alto grado de contaminación. 

• e) En una goma elástica, el factor más relevante cara a su variación de energía 
interna es la realización de trabajo estirando la misma. El trabajo realizado 
puede ponerse como el producto de la tensión r por la variación de longitud, por 
lo que puede considerarse ésta (o más en concreto su diferencia con la longitud 
en ausencia de tensión o elongación l = L — L 0 ) como coordenada de trabajo 
y la tensión como variable conjugada. En muchos de los problemas el volumen 
de la goma puede considerarse constante así como su composición por lo que la 
variación de energía interna está dada por 


dU = TdS + rdl 


(3.3) 
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3-2) Obtener las derivadas de la entropía con respecto a la energía interna, 
volumen y número de moles a partir de las correspondientes derivadas 
de la energía interna. 


* Método: utilizar las relaciones entre las derivadas de una función entre 
tres variables ligadas de la forma z = z(x,y). 




(3.4) 


• Para obtener (§§) VN hacemos uso de la primera de las relaciones expuestas 
(3.4) por lo que tenemos 



(3.5) 


que teniendo en cuenta que 




(3.6) 


• Para obtener (§^) UN hacemos uso de la segunda de las relaciones (3.4) obte¬ 


niendo 



que al ser 



= — P queda escrita como 



(3.8) 


• Análogamente para obtener ( JrM -) planteamos 

\ aly 3 J UV(N i 


UV{N) 



(3.9) 



) = Un nos lleva a 

aN i J SV(N) 


) 



(3.10) 
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3-3) Probar que 


c) (§s\ = c SL . 

> (dT)VN T ' 


— al ( du \ — C ■ hl (QS-) —C P . 

A ) V dT JVN — ^ v ’ u ' (dT/PN T ’ 

- d) í^) =-T(W\ (dT) 

U 1 (dT )u \dS )u (dT) S\dU) S 


*Método: Comparar la definición de capacidad calorífica con la expresión 
diferencial de la energía interna (o entropía). En el último apartado uti¬ 
lizar las relaciones entre las derivadas entre tres variables ligadas de la 
forma z = z(x,y). 

• a) Consideremos un proceso reversible infinitesimal a volumen constante en el 
que la temperatura varía dT. La variación de energía interna está dada por 

dU = 6Q- PdV (3.11) 

y el calor absorbido por el sistema está dado por 

6 Q = C v dT (3.12) 


Ahora bien, por ser el proceso a volumen constante, el término —PdV es nulo 
por lo que 

dU = C v dT (3.13) 


Por otra parte, en general en un proceso infinitesimal 

dU 




dT 


VN 


du\ „ r x 

7w) iv + T. 

TN j = i 


dN, 


dN¿ 


(3.14) 


3 / TV (N) 


Si el sistema es cerrado (dNj =0) y el volumen permanece constante (dV = 0) 
tenemos 

dU =(fff) dT (3.15) 


VN 


de donde deducimos comparando con (3.13) 


=c 

dTj VN v 


(3.16) 


• b) Si consideramos el mismo proceso que en el apartado anterior el calor absor¬ 
bido se escribe como SQ = TdS por lo que según (3.12) 


TdS = C v dT 


(3.17) 
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En general, en un proceso infinitesimal 



(3.18) 


que análogamente al caso anterior para un proceso en un sistema cerrado a 
volumen constante se reduce a dS — (§t)vn^ P or O 116 com P aran< io con 
(3.17) tenemos 


(a,s_\ = a 

W)vn t 


(3.19) 


• c) En el caso de un proceso reversible infinitesimal a presión constante 


SQ = CpdT 


(3.20) 


y por un razonamiento similar al del caso anterior 



(3.21) 


que para un sistema cerrado a presión constante se reduce a dS = (§t) pn ^ y 
teniendo en cuenta que SQ = TdS llegamos a 



(3.22) 


• d) Para probar la igualdad {jjf) s (§p) s tenemos en cuenta 

por las propiedades del cálculo infinitesimal que 


(dV 

\dT 


u 




(3.23) 


ahora bien 

fdV\ (l)y r (9V\ 
\9SJu (f£) 5 \dUj s 


(3.24) 


con lo que el resultado queda probado escribiendo la última relación en la forma 


(av\ _/av;\ / ot\ 
[diy) s ~W)s\9ü)s 


(3.25) 
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3 - 4 ) Obtener la ecuación de Gibbs-Duhem en la representación de la en¬ 
tropía 

Ud(^) + Vd(^)-jp / N j d(^ r ) = 0 . ( 3 . 26 ) 

3 = 1 


* Método: Diferenciar la ecuación de Euler y utilizar la expresión diferencial 
de la entropía. 

• Energía interna, volumen y número de moles constituyen un conjunto de va¬ 
riables extensivas independientes que describen de forma única el estado de un 
sistema homogéneo. Cualquier magnitud extensiva, y en particular la entropía, 
es una función homogénea de primer orden respecto de dichas variables indepen¬ 
dientes 

XS = S(XU, XV, XNi ,..., XN C ) 

y por tanto verifica el teorema de Euler 

s^Lu + ^v-Y^Nj 

JP 'JP / J ^ d 

3 =1 

• Si diferenciamos esta expresión obtenemos 

dS = Ud(^) + i dU + Vd(^) + ¿ N 0 d(^) ¿ ^dNj (3.29) 

3 = 1 3=1 

que teniendo en cuenta la expresión diferencial de la entropía dS = TpdU +-^dV— 
YPi-i dNj simplificando términos se llega a 

Ud^ + Vd^-jpNAjt)^ o ( 3 - 3 °) 

3 = 1 


( 3 . 27 ) 

( 3 . 28 ) 
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3-5) Obtener la forma del potencial químico de un gas ideal a partir de 
sus ecuaciones térmica PV = NRT y energética de estado U = cNRT 


* Método: Integrar la ecuación de Gibbs Duhem sustituyendo los valores de 
Uy V 


• La ecuación de Gibbs Duhem para un gas (sistema homogéneo) monocomponente 
se escribe como 


Ud(^) + Vd(^)-Nd(^) = 0 

(3.31) 

que podemos poner de la forma 


„/i, u 1. V P 

T ~~ ~Ñ T + iV T 

(3.32) 


Si sustituimos los valores obtenidos de las ecuaciones térmica y energética 

% = cRT de estado obtenemos: 


d(f) = cRTd(^) + (3.33) 

que al ser integrada entre los valores T 0l P 0 y los valores T, P conduce a: 

f-f= cri 4 +k KíI) (3 - 34) 

que podemos escribir como: 

»(T,P) = ^T + TcRln^+RTln(^^\ +RT\nP (3.35) 

Los primeros sumandos sólo dependen de la temperatura y corresponden al valor 
del potencial químico para una presión igual a la unidad. Por tanto la expresión 
anterior suele denotarse como 

¡i(T,P) = n(T,P =l) + RT\nP (3.36) 

o bien 

P) = rj'(T) + RTlnP (3.37) 


• Ejercicio adicional: 

Sustituyendo en la ecuación los valores de la ecuación térmica y energética (3.34) 

— Obtener el potencial químico en términos de T,V,N 
— Obtener el potencial químico en términos de U,V y N 
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3-6) Obtener la ecuación U — U(S, V, N) para un gas ideal a partir de las 
ecuaciones térmica PV = NRT y energética de estado U = cNRT y la 
ecuación ¡j,(T,V,N) obtenida en el problema anterior. 


* Método: Sustituir en la ecuación de Eider de la representación entrópica 
los valores de y ^ en términos de U, V,N. 


• De las ecuaciones térmica y energética podemos obtener 

J_ _ cNR p _ NR 

TU TV 


(3.38) 


• La expresión de ^ como función de U,V podemos obtenerla de (3.34) utilizando 
la ecuación térmica dada en (3.38) obteniéndose 


d 

T 


do 

T 

-■*- n 


= cR ln 



+ R ln 



(3.39) 


Sustituyendo (3.38),(3.39) en la ecuación de Euler 


S=^U+^V-^N 


(3.40) 


se tiene 


S = —jj—u 


NR, 

ir 


do 


v ~ I ™ +cRln > Wü 


U 0 N 




N (3.41) 


que simplificando se reduce a 


S = (c+l)NR- N 



R ln 





(3.42) 


• Teniendo en cuenta que de la ecuación de Euler y de (3.38) deducimos que 


&o / i 1 \ n do 

Wr iC+l)R -T. 


(3.43) 


la ecuación queda escrita como 


S_ 

N 


S ° 




ín_\- { c+1) 

W 


(3.44) 


• Ejercicio adicional: 

Obtener la forma de la ecuación en la representación de la energía U = U{S 1 V 1 N). 
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3-7) A 520° C y 500 atmósferas de presión una mezcla estequiométrica de 
N 2 y H -2 según la reacción 

N 2 (g) + 3H 2 (g) <=> 2 NH 3 (g) (3.45) 

alcanza el equilibrio con un 50% de H 2 . Determinar el grado de avance 
de la reacción y las fracciones molares de los restantes componentes 


* Método: Relacionar el número de moles finales con el grado de avance y 
tener en cuenta que la fracción molar de hidrógeno en el equilibrio es 0.5 


• El grado de avance está definido de manera que su variación corresponde al 
cociente entre la variación del número de moles de cada componente y el corres¬ 
pondiente coeficiente estequiométrico d£ = —— por lo que en un proceso finito 
en el que el número de moles de cada sustancia varía de un estado inicial N° a 
un estado final Nj tenemos para cada sustancia: 

Nn 2 -n^ 9 _ c c n h 2 -n^ 9 _ c c JV W h3-jvj Hí _ t r 

_1 S so _3 S So 2 S So 


• Asumiendo condiciones estequiométricas, es decir el número de moles inicial es 
igual a los coeficientes estequiométricos (Nf = Vj) para los reactivos y cero para 
los productos y tomando grado de avance cero (£ 0 = 0) para el estado en el que 
sólo tenemos reactivos se obtiene 

Nn 2 — 1 t Nh 2 ~ 3 t _ £ 

_1—S _3~S 2 ^ 

de donde 

N N2 = — f + 1 N H2 = —3£ + 3 N NH3 = 2^ 


(3.46) 

(3.47) 


• El número de moles totales será N = + 1 — 3£ + 3 + 2£ = —2^ + 4 por lo que 

las fracciones molares en el equilibrio están dadas por 


_ _ n n 2 

%N 2 — ]y 


-£+1 

4-2Í 


%H 2 


n h 2 

N 


— 3g+3 
4-2J 


x N h 3 =- 


JV wf 


N 


4-2? 


• La solución de la ecuación que iguala la fracción molar de H 2 dada por el enun¬ 
ciado del problema nos da 

XH2 = = °' 50 (3 - 48) 

de donde deducimos un valor del grado de avance £ = 0. 5 lo que nos da para las 
restantes fracciones molares 

x N2 = « 0.166 x NH 3 =^Bt ~ 0- 33 (3.49) 

es decir, el equilibrio se alcanza con un 50% de H 2 un 17% de N 2 y un 33% de 
NH 3 
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Capítulo 4 

Condiciones de equilibrio y 
equilibrio de fases 

4.1 Objetivos 

• Comprender la aplicación de las condiciones de equilibrio de fases. 

• Comprender la aplicación de la regla de las fases. 

• Familiarizarse con la utilización de la Ecuación de Clausius-Clapeyron. 

• Comprender la resolución de procesos en la región de saturación. 



Benoit-Paul-Emile Clapeyron (1799-1864). Ingeniero de minas francés que se especializó en el diseño de locomo¬ 
toras y ferrocarriles. En 1824 publicó su “Memoria sobre la potencia motriz del calor”, en la que expresaba de una 
forma analítica las ideas desarrolladas por Sadi Camot diez años antes y que fue una obra clave para el nacimien¬ 
to de la Termodinámica como rama de la Ciencia. 


51 



4.2 Problemas 


4-1) Dos cilindros conteniendo un determinado fluido se encuentran co¬ 
nectados por un pistón diatermo de la forma que indica la figura. 
Suponiendo despreciables los efectos del entorno (el sistema se en¬ 
cuentra aislado) determinar las condiciones de equilibrio del sistema. 



* Método: Plantear las condiciones de equilibrio teniendo en cuenta la li¬ 
gadura que la presencia del pistón supone para los volúmenes del fluido 
contenido en ambos cilindros. 

• Consideremos un estado de equilibrio del sistema, el fluido en cada cilindro estará 
caracterizado por su energía interna, volumen y número de moles. Denotemos 
con superíndice (1) las variables del fluido en el primer cilindro y con superíndice 
(2) las variables del segundo. 

• Variables que describen el subsistema 

— Energía interna, U W 

— Coordenadas de trabajo: Volumen V 1 ' 1 ' 1 

— Número de moles JV' 1 * 

— Análogamente el subsistema estará descrito por su energía interna 
U*' 2 ),volumen V^ 2 ) y número de moles 

• Las ligaduras a las que está sometido el sistema son: 

— El sistema total está aislado 


[/(i) + u&) = cte 


— El número de moles en cada cilindo es constante 

N W = cte N^ = cte 


(4.1) 


(4.2) 
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— Los volúmenes aparecen ligados por la restricción de que la longitud total 
de ambos cilindros más el pistón es constante. Si denotamos por L\ y L -2 
las longitudes de la región de ambos cilindros ocupada por el fluido y por 
L la longitud del pistón, la correspondiente ligadura estaría escrita de la 
forma 

L\ L *2 L = cte (4-3) 

— En términos de los volúmenes y teniendo en cuenta que en cada cilindro 
V W = L¿A W siendo su sección, escribimos 


yW 1 v (2) 

+ AP) 


( 4 . 4 ) 


• En un proceso virtual desde el estado de equilibrio la variación primera de la 
entropía dS = dS (1) + dS^ en términos de las variables de cada uno de los 
subsistemas está dada por 

dS =té) iU(i) + w> dVÜ) -^ dN(i) (4 - 5) 

3 =1 

ahora bien, por las condiciones de ligadura 

dUM = -dUW dNW = dNW = 0 ( 4 - 6 ) 


Por tanto podemos elegir como variables no ligadas la energía interna y el vo¬ 
lumen del primer sistema ya que las del segundo sistema pueden ponerse en 
términos de las del primero. 

• En estas condiciones la variación de entropía en términos de las variables inde¬ 
pendientes está dada por 


dS = 



rw) + 



A ( 2 ) p( 2 ) \ 

AWT( 2 )) 


dV^ 


( 4 . 7 ) 


Como en virtud del principio de entropía máxima en un proceso virtual desde 
el equilibrio dS = 0 independientemente de los valores de dU ^ y dU^ 2 \ los 
coeficientes de dichas variaciones deben ser nulos por lo que 


1 

1 i 

Tí 1 ) 

Tí 2 ) ) 

” p(U 

J 4Í 2 )p( 2 ) 

Tí 1 ) 

.4(!)T( 2 ) 


= 0 


p(!) = p(2) 

j\(l)p( 1 ) _ j\(2)p(2) 


• Es decir en el equilibrio existe igualdad de temperaturas, mientras que las pre¬ 
siones están en relación inversa a sus secciones, es decir si un cilindro tiene mitad 
de sección que el otro su presión será el doble. 
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4-2) Obtener las condiciones de equilibrio para una gota de agua en equili¬ 
brio con su vapor, teniendo en cuenta el efecto de la tensión superficial. 
(Supóngase la gota esférica). 


*Método: Plantear las condiciones de equilibrio teniendo en cuenta la con¬ 
tribución a la energía interna debida a la tensión superficial y la ligadura 
entre el volumen y la superficie de la gota. 

• Consideremos un estado de equilibrio de un sistema aislado constituido por una 
gota de agua en equilibrio con vapor de agua. 

— Las variables que describen a la gota de agua (subsistema J3 1 ) son: 

* Entropía SW 

* Coordenadas de trabajo: Volumen V^ 1 *, Superficie AW 

* Número de moles 

— Las variables que describen al vapor de agua (subsistema ^ 2 ) son: 

* Entropía S*' 2 ) 

* Coordenadas de trabajo: Volumen V^ 2 * 

* Número de moles N^ 2 '. 

— Aplicaremos al sistema las condiciones del principio de energía interna 
mínima. 

”El estado de equilibrio del sistema es un estado de mínima energía interna 
de entre todos los estados de igual entropía compatibles con las restantes 
ligaduras a que esté sometido el sistema. ” 

— Condiciones de ligadura: 

£-(i) + £(2) = cte V W + y( 2 ) = cte + N& = cte 


El área de la gota no es una variable independiente del volumen, la relación que 
existe entre ambas viene dada por: 


A = 4nr 2 

= 17rr 3 

(4.8) 

es decir 

dA = 87r rdr 

dV^ = 47r r 2 dr 

(4.9) 

con lo que 

dA = 

-dV^ 

r 

(4.10) 
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• El principio de energía interna mínima establece que en todo proceso virtual a 
S = cte desde el estado de equilibrio el primer término en el desarrollo de A U, 
es nulo (dU = 0) por lo que se tiene 

dU = TWdSW - P (1) dl/ (1) + odA + ^dN (1 > + T^dS^+ 

_p(2)y(2) + Ai (2) dAr (2) = Q ( 4 .n) 

y atendiendo a las condiciones de ligadura 

dSW = -dS « dVW = -dV<Ü dA = 2 dV W dN<Ü = -dN^ 

r 

con lo que la expresión (4.11) se reduce a 


dU = (T (1) - T (2) ) dS (1 ) - ^P (1) - P (2) - y ^ dE (1) + 

+ (a¿ (1) - {i (2) ) dN (1) = 0 (4.12) 

siendo variables independientes. Por lo tanto las condiciones de 

equilibrio son 

T (1) = P (2) (4.13) 

P( i) = P( 2 > + — (4.14) 

r 

p (1) = a¿ (2) (4.15) 


• Observamos que la temperatura es la misma en el interior de la gota que en el 
exterior, pero no ocurre lo mismo con la presión. La igualdad de los potenciales 
químicos debemos interpretarla de la siguiente forma: 

/r( 1 )(T, p( 2 ) + ^) = f¿ 2 \T,pW) (4.16) 


Observamos que: 

— A medida que aumenta el radio de la gota la diferencia de presiones exterior 
e interior se hace menor y el efecto de la tensión superficial es despreciable. 

— Fijada una presión exterior existe un valor r del radio de las gotas de 
líquido que hace que éste se mantenga en equilibrio. 

• Ejercicio adicional: 

— Discutir cómo varía el potencial químico de la gota con su radio, y teniendo 
en cuenta que la materia fluye de la región de mayor a la de menor potencial 
químico discutir que ocurre si una gota tiene un radio mayor o menor que 
el radio de equilibrio. 
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4-3) En las condiciones del problema anterior encontrar el valor del radio 
de una gota de agua que se mantiene en equilibrio con su vapor a la 
temperatura de T = 20°C' con una presión parcial del vapor de agua 
p( 2 ) = 20 mmHg siendo la presión de saturación para dicha temperatura 
p(°) = YJ.^mmHq y el coeficiente de tensión superficial del agua a = 
0.0728 Jm~ 2 


* Método: comparar el valor obtenido en el equilibrio gS l \T, P^ + ^r) = 
/i( 2 )(T, P^ 2 )) con el valor obtenido para una cantidad grande de agua r —» oo 


• El estudio de las condiciones de equilibrio de una gota de agua con su vapor a 
una presión exterior P* 2 ) teniendo en cuenta los efectos de la tensión superficial 
nos lleva a que la presión en el interior P^ 1 ' está dada por P 1 - 1 ' 1 = P* 2 ) + — y 
que entre los potenciales químicos se tiene la relación: 


/i¡ ¡) (T, P (2) + y) = 

^\T,PW) 

(4.17) 

Para obtener una forma explícita del radio de la gota podemos comparar el caso 
anterior (r finito) con el caso r —*• oo es decir una cantidad grande de líquido, 
para el que se obtiene 

$(T,PW)=Ú' 

;) (T,P (0) ) 

(4.18) 

donde P*- 0 ' es la presión de equilibrio líquido vapor cuando los efectos de la 
tensión superficial son despreciables. 

Restando ambas expresiones 



H^(T, PW + -)- P (0 >) = a4 1 

r 

’\t, pP^-zi^r,?* 0 )) 

(4.19) 

y teniendo en cuenta que al ser 



’=(£ 

), 

(4.20) 

podemos escribir 

r p' 


p{T,P)-p{T,P') = ^ 

1 v(T, P)dP 

(4.21) 

Por tanto la expresión (4.19) se puede poner como 


p(2)_|_ 

/ r v {l \T 1 P)dP = 

Jp(°) j 

,p< 2 ) 

/ V ( v \T, P)dP 

lp( 0) 

(4.22) 


donde y son los respectivos volúmenes molares del líquido y del vapor. 
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• Para la fase líquida el volumen molar i>0) (T, P) no varía demasiado con la presión 
por tanto puede considerarse constante, y para la fase de vapor podemos tomar 
la aproximación de gas ideal 


V W(T,P) = — 


con lo que (4.22) queda 


( p(2) + V 


p(°)^ v V)(T,pW) = i?Tln^j 


(4.23) 


(4.24) 


de donde podemos despejar r en términos de la presión P^ 2 \ la presión PÍ°) 
de saturación para una gran cantidad de líquido a esa temperatura, el volumen 
molar del líquido v^ l \ el coeficiente de tensión superficial a, y la temperatura T 
obteniendo 


<JV 


(0 


v (l)(p(2) _p(0)) 


(4.25) 


• Numéricamente 
cr = 0.0728 Jm" 2 

pi l ) = lgcm~ 3 x 10 ~ 3 kgg~ 1 x 10 6 cm 3 m -3 = 10 3 kgm~ 3 

= 1 - * « 

P< 2 > =20mmHgx ¿,^,,,-1 * 101.3 x itfPaatm- 1 =2665.8 Pa 
p(°) = 17.5 mmHg x atm -1 x 101-3 x 10 3 PaatmT 1 = 2332.6 Pa 

T = (20 + 273.15) = 293.15 K 

u(p( 2 ) - P ( °)) = 1.8 x 10 _5 m 3 mo/ _1 x 333.2 Pa = 5.9976 x 10~ 3 JmoP 1 
RT ln ^ = 8.314Jmor 1 P:- 1 x 293.15 K x ln 2 ^; 3 g° = 325. 42.Jm.oP 1 

• Por tanto para el radio se obtiene 


9 0.0728 Jm 2 x 1.8x10 5 m 3 mol 9 

"" 5.9976 x 10~ 3 J mol~ l — 325.42 J moP 1 


es decir, el radio en equilibrio de la gota de agua es de unos pocos nanómetros. 
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4-4) El punto de ebullición normal del butano es de T a = —0.5 °C. Teniendo 
en cuenta que el calor de vaporización es de l v = 23.4/cJ mol" 1 estimar 
la presión de vapor de una botella de butano a la temperatura de 25°C. 
Un mechero de butano se lleva en el bolsillo adquiriendo un valor de 
la temperatura próximo a la temperatura corporal 37 °C. Calcular la 
presión en el interior del mechero. 


* Método: Aplicación directa de la ecuación de Clausius-Clapeyron. 




La ecuación de Clausius-Clapeyron para el caso de la vaporización se escribe de 
la forma 


dP _ l v P 
dT ~ Kfi 


(4.27) 


donde l v es el calor de vaporización y i? la constante de los gases. 


• Si suponemos que el calor latente es constante e integramos entre dos estados 
sobre la curva de equilibrio de fases B 0 (P 0 , T a ) y B(P, T ) llegamos a 


ln 



¿12 (l 
R \T 



(4.28) 


• Con los datos numéricos aportados por el problema 


¿12 
R 


23.4 x 10 3 J mol- 1 
8.314 J moP 1 K~ l 


-2814.5 K 


(4.29) 


¿12 1 
R T 0 


2814.5 1< 
272.65 A' 


10.323 


(4.30) 


• Por tanto podemos establecer que la variación de la presión de vapor del butano 
con la temperatura puede escribirse de la forma 


P = P 0 exp ^-2814.5 A'i + 10.323^ (4.31) 

que nos ofrece para los casos T = 298.15A' (25°C) 

P = 1 atm x exp (—2814. 5 K 298 1 15K + 10.323) = 2.42 atm 

y para T = 310.15 1\ (37 °C) 

P = 1 atm x exp (—2814.5 K 310 \ 5K + 10. 323) = 3.48 atm 
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4-5) El carbonato cálcico se presenta en la naturaleza en dos fases sólidas 
diferentes: la calcita y la aragonita de densidades respectivas p ca idta = 
2.71 g/cm 3 ; Paragonita = 2.93 g / cm 3 . Medidas de las capacidades ca¬ 
loríficas nos permiten encontrar la entropía de ambas fases s ca idta = 
0.222 calg^K^ * 1 y s ara gonita — 0.212 calg^ 1 K~ 1 . A la temperatura T = 
25 °C ambas fases coexisten a la presión P = 3737 bar. ¿Cuál es el au¬ 
mento de presión por cada grado de aumento de temperatura a lo 
largo de la curva de coexistencia?. 


* Método: Aplicación directa de la Ecuación de Clausius Clapeyron 


• El problema es una aplicación directa de la ecuación de Clausius-Clapeyron 


dP s^ 2 ) — sW 
dT = u( 2 ) - vA) 


(4.32) 


• Si tomamos como fase (2) la calcita y por fase (1) la aragonita y teniendo en cuen¬ 
ta la relación entre volumen molar y densidad v = -y los respectivos volúmenes 
molares serán 

(2) _ íoo g mol _ 36.9 cm 3 mo¡ _1 x 10 -6 m 3 cm -3 = 36.9 x 10 ~ e m 3 mol~ 1 

2.93 g cm á 

i A 1 ) = i°0ff moi = 34.1 cm 3 mol~ 1 x 10 -6 m 3 cm -3 = 34.1 x 10 _6 m 3 moZ _1 

z.tlg cm á 

y las correspondientes entropías 

= 0.222 cal g~ l K~ x x 4.18 J cal~ l x 100 gmoT 1 =92.8 Jmol^K -1 
s (■*■) = 0.212 cal g^K^ 1 x 4.18 J cal~ l x 100 g moT 1 =88.6 Jmol^K -1 


• Sustituyendo en la ecuación de Clausius-Clapeyron 


dP 92.8 J mol 1 K 1 — 88.6 J mol 1 K 1 

dT 36.9 x 10 ~ 6 m 3 mol~ 1 — 34.1 x 10 ~ e m 3 mol~ 1 


1.5x10 6 PaK~ 1 (4.33) 


donde hemos tenido en cuenta que J m 3 = N m 1 = Pa. Por tanto, como 

1 Pa = 10 _5 6ar 

dP 

— = 15 bar K- 1 (4.34) 

• Es decir la presión aumenta 15 bar por cada grado, por lo que si la temperatura 

aumenta 1°C (= 1A') la presión aumenta de 3737 a 3752 bar. 
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4-6) Se discute si entre los diferentes factores que afectan a la profundidad 
de la capa de hielo en los mares polares la variación de la temperatura 
de fusión con la presión debido a la profundidad es un factor relevante. 
Analizar este hecho calculando la variación de la temperatura de fusión 
del agua con la profundidad teniendo en cuenta que la densidad del 
hielo a 0 °C es de 0.917 g/crrc’y la del agua es 0.9998 g/cm 3 y el calor 
latente de fusión del hielo es 80 cal/g. 


* Método: Calcular la variación de presión con la profundidad y aplicar 
la ecuación de Clapeyron. 


• Asumiremos que la presión a una determinada profundidad está dada por la 
expresión hidrostática P = P 0 + pgz donde P 0 es la presión en la superficie 

• Así pues, la variación de la presión con la profundidad está dada por 

f = « <“ 5 > 


• La variación de la temperatura de fusión con la presión nos viene dada por la 
ecuación de Clausius-Clapeyron 


dT T(i/ 2 )-u«) 

dP T f 


(4.36) 


donde son los respectivos volúmenes molares del agua y del hielo, y lf 

es el calor latente de fusión. Por tanto la variación de la temperatura de fusión 
con la profundidad podemos escribirla como 


dT _ dT dP 
dz dP dz 


(4.37) 




Analicemos en primer lugar la variación de la temperatura de fusión con la 
presión a partir de la ecuación de Clausius-Clapeyron. La relación entre volumen 
molar y densidad la escribimos de la forma v = -S = w = — donde M es la 

J N M P 

masa molecular. 


dT TM {^ - ^y) 
dP l f 


(4.38) 


El calor latente en esta fórmula es el calor latente molar, ahora bien, el enunciado 
del problema nos da el calor latente específico (por unidad de masa), la relación 
entre ambos está escrita de la forma 


lf(calg 1 )xM(gmol 1 )=lf{calmol 3 ) 


por lo que en la ecuación anterior podemos poner 

dT _ T {-^Y - 717 ) 
dP l f 


(4.39) 


(4.40) 
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• Numéricamente el calor latente 

lf = 80 cal g -1 = 80 cal g~ l x 4.18 J cal~ l x 1000 g kg~ x = 334.4 x 10 3 J kg~ 1 

• Las respectivas densidades 

p = 0.917 g cm~ 3 x 10 ~ 3 kgg~ 1 x 10 6 c?n 3 m~ 3 = 917 kgm~ 3 
= 0.9998 gcrrC 3 x 10~ 3 kgg^ 1 x 10 6 cm 3 m -3 = 999.8 kgm~ 3 


• Así pues, sustituyendo datos 

dT _ T{ 999 g fc gm -3 ~ 917 kgm~ 3 ) 

dP 334.4 x 10 3 J kg~ 1 

— = ( " 9 - 8fc9m ~ a ~ 917fc g^ a) = (—2. 7007 x lO-^Pu- 1 ) x T 
¿P 334.4 x 10 3 J fcg -1 ' ' 

• Para el punto triple esta expresión nos ofrece una pendiente 

dT 

— = -7.38 x 10~ 8 A' PíP 1 = -7.47 x lO^Katm^ 1 
dP 

que también podemos poner como 

^ = -1.3555 x 10 7 Pu A' -1 = 133.8 aíro AT _1 


(4.41) 

(4.42) 

(4.43) 

(4.44) 


• Suponiendo que esta pendiente no varía demasiado con la temperatura, la varia¬ 
ción de la presión con la profundidad ^ = pg nos viene dada por 

dP 

— = 999.8 kgm~ 3 x 9.81 ms” 2 = 9808 ParrC 1 (4.45) 

por tanto 

— = — — = -7.38 x 10 -8 A' PíA 1 x 9808Pam _1 = -0.7 x lO^A'm^ 1 
dz dP dz 

Es decir, que la variación predicha de la temperatura de fusión es de 0.7 K por 
cada kilómetro de profundidad. 
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4-7) La regla de Trouton establece que la variación de entropía en un proce¬ 
so de vaporización para un fluido en el que las interacciones molecula¬ 
res no son demasiado intensas está dada por S M_ S (0 =88 J K~ 1 mol~ 1 . 
Estimar el calor latente de vaporización y obtener una ecuación para 
la curva de equilibrio líquido vapor para una sustancia que verifique 
esta regla. 


* Método: Relacionar el calor latente con la variación de entropía e integrar 
la ecuación de Clausius-Clapeyron. 

• El calor latente de vaporización está dado por 



donde T a es la temperatura de ebullición y «W y s(0 son las entropías molares 
del vapor y del líquido respectivamente. 


• Si aplicamos la regla de Trouton nos queda 


l v —88 J K 1 mol 1 x T a 


(4.47) 


• Sustituyendo esta expresión en la ecuación de Clausius-Clapeyron e integrando 
llegamos a: 



(4.48) 


y tomando el valor de R = 8.314 J mol 1 K 1 y P 0 = 1 atm llegamos a 



(4.49) 


donde la presión está dada en atmósferas. 


• En el caso del agua T 0 = 373.15 K por lo que 


3949.8. 

P = exp(10. 585 -——) 


(4.50) 


siendo T la temperatura en Kelvin 
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4-8) En el caso de la sublimación y la vaporización, la variación del calor 
latente con la temperatura puede ponerse de la forma 


í dl v 
\dT 


sat 




donde c p y Cp 1 son las respectivas capacidades caloríficas a presión 
constante del vapor y de la fase condensada (líquido o sólido) y () t 
denota que la derivada tiene lugar sobre la correspondiente curva de 
equilibrio de fases. En el punto triple del agua (T = 273.16 K, P = 
4.58 mmHg) el calor latente de vaporización l v = 2.5 x 10 6 J kg~ 1 y las 
correspondientes capacidades caloríficas son c ( p l) = 4218 J kg~ 1 K^ 1 y 
Cp^ = 1850 J kg~ 1 K~ 1 . Encontrar una ecuación para la dependencia 
de la presión de vapor del agua con la temperatura considerando la 
variación del calor latente con la temperatura. 


* Método: Obtener una dependencia lineal del calor latente con la tempera¬ 
tura y sustituir en la Ecuación de Clausius-Clapeyron. 

• Si la derivada del calor latente está dada por 

(451) 


suponiendo que las capacidades caloríficas son constantes en el intervalo de tem¬ 
peraturas objeto de estudio integrando la expresión anterior podemos poner 

lv(T) = l v (T a ) + (4”) - 4'))(r - To) (4.52) 

Es decir el calor latente tiene una dependencia con la temperatura de la forma 


l v (T) = A + BT (4.53) 

donde 

A = l v (T 0 )-{c^ - C V)T 0 (4.54) 

B=(cW-cU) (4.55) 


• Si sustituimos esta expresión en la ecuación de Clausius-Clapeyron e integramos 


llegamos a 



1 dP 
PdT' 


dT 



A + BT 
RT 2 


dT 



Ai 1 1 \ B T 

R\T~Y 0 ) + R ln Y 0 


(4.56) 


(4.57) 
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que podemos desglosar como 


ln 



A 1 
RT 


B A 

R ln + ~KT 0 


I 1 " 3 ’- 


(4.58) 


En nuestro caso con los valores numéricos del enunciado del problema 

c£° - 4 l) = (1850 J kg~ 1 K~ 1 - 4218 J kg^R- 1 ) = -2368.0 Jkg^R- 1 

A = 2.5 x 10 6 Jkg- 1 - (-2368.0 J kg~ 1 K~ 1 ) 273.16A = 3146 x 10 3 JfccT 1 

A = 3146 x 10 3 J kg x 18 g moR 1 x 10 ~ 3 kg moR 1 = 56.64 J moR 1 

i? = (1850 J kg~ x K^ 1 — 4218 J kg~ 1 K~ 1 ') x 18 g moR 1 x 10~ 3 fc(/ = —42.624 J mol 

| = 6813 A' 

|=5.1268 

^ | lriT 0 = 24.941+ 28. 757 = 53.69 

Por tanto la ecuación buscada, dada la presión en mmHg corresponde a: 


P = 4.58 exp 


-6813- -5.12lnT+53.69 


(4.59) 


Ejercicio adicional. 

— Contrastar la expresión obtenida con el valor tabulado de P = 17.5 mmHg 
a T = 20 °C y de 23.75 mmHg a 25°C, y con los valores que se hubieran 
obtenido en el supuesto de que se considerara el calor latente constante e 
igual al valor en el punto triple. 



Las propiedades de la condensación del vapor, así como otro buen número de magnitudes experimentales, entre 
ellas las desviaciones del comportamiento de gas ideal, fueron ampliamente estudiadas por el francés Henri-Victor 
Regnault (1810-1878) quien además fue uno de los pioneros en la técnica de la fotografía. 




4-9) Vapor de agua saturado a la temperatura de T\ = 150°C y presión 
P\ = 4.758 6a?’ se expande adiabática y reversiblemente en una turbina 
hasta que su presión desciende a la presión atmosférica P -2 = 1.0136a?’ 
(T = 100°C'). Calcular la fracción de vapor que se condensa. En el esta¬ 
do inicial la entropía del vapor saturado s^(150°C) = 6.8379 kJ kg~ 1 K~ 1 
y en el estado final .s^' 1 (100°C) = 7.3549 kJ kg~ l K^ 1 mientras que pa¬ 
ra el líquido saturado s^(150°C) = 1.8418 kJkg^K -1 , S2* *(100 o C) = 
1.3069 kJkg-'K- 1 


* Método: Observar que el proceso es isoentrópico y calcular las variaciones 
de entropía. 

• El proceso de expansión es isoentrópico por lo tanto en el estado final la entropía 
será igual a la entropía del estado inicial. Observamos que a 100°C' la entropía 
del vapor saturado es mayor que la entropía del estado inicial (vapor saturado 
a 150°C) por tanto en éste parte del vapor debe estar condensado y por tanto 
hay una coexistencia de líquido y vapor. 

• Si denotamos por m ^ la fracción de masa de vapor que se condensa en la ex¬ 
pansión, la entropía en el estado final debe ser 

S 2 = m^s® + ^1 — rn^ j (4.60) 

la cual debe ser igual a la entropía en el estado inicial, que corresponde al estado 
de vapor saturado 



• Así pues planteamos 

+ ^1 — s 2 ’^ (4.62) 

• Numéricamente 

6.8379 kJkg^R- 1 = m w 1.3069 kJkg^R- 1 + (1 - m m ) 7.3549 kJkg~ 1 I\~ 1 
que nos da la solución 

m (l) = 0.09 (4.63) 

• Es decir el 9% del vapor de agua se condensa en la expansión. 

• Ejercicio adicional. 

— Comparar el valor de la entropía del líquido saturado a la temperatura T^ 
con el valor calculado en términos del correspondiente a la temperatura 
T ( 1) suponiendo que el calor específico del agua saturada es constante c v = 
1 cal g” 1 /! -1 
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4-10) Se conecta una bomba de vacío sobre un recipiente adiabático que 
contiene inicialmente 500 g de agua a 100°C. Mediante la acción de 
la bomba el descenso de presión provoca la evaporación del agua y 
por tanto la masa que permanece en el recipiente se enfría. Teniendo 
en cuenta que el calor latente del vaporización del agua es l v (cal/g ) = 
798. 21 — 0. 7T donde T es la temperatura Kelvin, determínese la masa 
de agua que queda en el calorímetro cuando la temperatura alcanza 
el punto triple i = 0.1°C. 


*Método: Considerar que el proceso de evaporación es isoentrópico y rela¬ 
cionar las correspondientes variaciones de entropía en cada etapa infinite¬ 
simal. 


• El proceso de bombeo podemos modelizarlo de la forma siguiente: Un pistón en 
contacto con el líquido se desplaza aumentando el volumen disponible y provo¬ 
cando la vaporización de una cierta cantidad de líquido. El vapor es extraído 
del recipiente y se repite de nuevo el proceso. 


• Supongamos que en cada movimiento del pistón se evapora una cantidad dm 
de líquido. Si el proceso es adiabático la variación de entropía es nula. Dicha 
variación de entropía corresponde a la variación de entropía del líquido dS ^ = 
mcp^f- donde m es la masa, c p el calor específico a presión constante y dT la 
variación experimentada de temperatura. Por otra parte la variación de entropía 
del líquido que pasa a vapor está dada por dS^ = —dmijt donde l v es el calor 
latente específico y — dm representa la fracción de masa vaporizada (dm es la 
variación de masa del líquido) por lo que podemos escribir la ecuación 


dT j 

mc p —— dm 


lv(T) 

T 


= 0 


que podemos reescribir de la forma 


dT dm 

Cp UT) = 


(4.64) 


(4.65) 


• La ecuación en variables separadas después de integrada conduce a 


f Tf dT 


lv(T ) 


ln 


mj_ 

mi 


(4.66) 


/» 273.16 rj-i 

1Ca,S " K "L» ( 798.21 - 0 . 700 ^»-. 

de donde obtenemos que la masa final de líquido presente en 
m/ = 419.72 g. 


ny_ 

5005 


(4.67) 


el calorímetro es 
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Capítulo 5 

Ecuación fundamental y 
Potenciales Termodinámicos. 

5.1 Objetivos 

• Familiarizarse con el concepto de ecuación fundamental. 

• Obtener a partir de la ecuación fundamental distintas ecuaciones de estado 

• Obtener a partir de la ecuación fundamental en una representación la ecuación 
fundamental en otras representaciones. 

• Familiarizarse con el uso de las relaciones entre magnitudes termodinámicas y 
en particular las relaciones de Maxwell. 

• Estudiar las magnitudes asociadas a distintos procesos particulares a partir de 
la información contenida en la ecuación térmica y las capacidades caloríficas. 



Hermann von Helmholtz (1821-1894) fue un médico alemán interesado desde muy pronto por los problemas físi¬ 
cos. Su trabajo “Sobre la conservación de la fuerza (Kraft)” publicado en 1847 es la primera formulación explícita 
del principio de conservación de la energía. 
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5.2 Problemas 


5-1) La ecuación fundamental de un gas ideal U = U(S, V, N) está dada por: 


U = U a + 




Cg+i) 

exp 





(5.1) 


• a) Obtener su ecuación energética de estado U = U(T,V, N). 

• b) Obtener la ecuación térmica P = P(V,T,N) 


*Método: a) A partir de U = U(S,V,N) obtener por derivación T = 
T(S, V, N) y eliminar S entre ambas ecuaciones b) A partir deU = U(S, V, N) 
obtener por derivación T = T(S, V, N ) y P = P(S, V, N) y eliminar S de 
entre estas últimas ecuaciones para obtener P = P(T , V, N) 


• a) Derivando la ecuación fundamental (5.1) respecto de la entropía T = (j]^ )v. n 
obtenemos T = T(S,V,N) de la forma: 



La entropía podemos eliminarla dividiendo entre sí las ecuaciones (5.2) y (5.1) 
con lo que obtenemos 

\ = NcR (5.3) 

o lo que es lo mismo 

U = NcRT (5.4) 

• b) Derivando la ecuación fundamental (5.1) respecto del volumen P = 
llegamos a 


P = U o — 

cV n 




c+1 

exp 





(5.5) 


Para obtener 
siones (5.5) y 


la ecuación térmica P = P(T,V,N) dividimos entre si las expre- 
(5.2) resultando 


P _NR 
T — ~V~ 


(5.6) 


o lo que es lo mismo 

NRT 
~ V 


(5.7) 
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5-2) La radiación electromagnética en una cavidad en equilibrio con las 
paredes de ésta o radiación del cuerpo negro puede considerarse como 
un sistema expansivo (gas de fotones) de ecuación fundamental: 

U = aS 4/3 V- 1/3 (5.8) 

donde a es una constante dada por a = (|) 4 ^ 3 (¿) 1 ^ 3 ? en I a que c es I a 
velocidad de la luz y a la constante de Stefan Boltzmann. Obtener sus 
ecuaciones energética U = U (T, V ) y térmica P = P(T , V) de estado. 


*Método: a) Obtener T = T(S, V) por derivación de U = U(S, V) y 
deducir U = U(T,V ) eliminando S entre ambas ecuaciones, b) Obtener 
T = T(S,V) y P = P(S,V) por derivación de U = U(S,V) y deducir 
P = P(T , V) eliminando S, U entre ambas ecuaciones 


• Para obtener U = U(T, V) obtenemos T = T(S,V) derivando (5.8) respecto de 
S 

T ^ { ^ )v = í aSl/3V ~ 1/3 (5 - 9) 

Para eliminar S podemos despejarla de (5.8) 

S = a- 3 / 4 P 1/4 t/ 3/4 (5.10) 


y sustituirla en (5.9) 



T =^a (a- 3 / 4 P 1/4 l/ 3/4 ) V3 V- 1/3 

(5.11) 

llegando a 

T = -a 3/4 V- 1/4 U 1/4 

3 

(5.12) 

de donde 




C/=(j) a->VT‘ 

(5.13) 


sustituyendo el valor de la constante a 3 =(f) 4 (¿) 4 se obtiene 


de donde vemos que la densidad de energía es proporcional a la potencia cuarta 
de la temperatura. 
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• Para obtener la ecuación térmica P = P(T, V ) derivamos (5.8) respecto de V 
para obtener P = P(S, V ) 

p = " O » = -° s4/3 (- 5 l/ ' 4/3 ) = 5“ s4/3 (' / "‘' 3 ) P-is) 

Para eliminar S despejamos en (5.9) S = (|) 3 a -3 V T 3 y sustituimos en (5.15) 
llegando a 



que queda como 

p= KD a ~ 3r4 (5 ' 17) 

que teniendo en cuenta que a -3 = (|) 4 y se llega a 

P=- -T 4 (5.18) 

3 c 

la presión de la radiación del cuerpo negro es únicamente función de la tempe¬ 
ratura. 

• Ejercicio adicional. 

— Teniendo en cuenta que la constante de Stefan Boltzmann está dada por 
a = 5.67x 10 -8 J s _1 m” 2 K~ 4 y la velocidad de la luz c = 299.8x 10 3 fems -1 
calcular la densidad de energía y presión de radiación para los siguientes 
casos: Radiación de fondo del universo T = 2.72 A',Temperatura media 
aproximada de la superficie de la tierra T = 15 °C, Temperatura en la 
superficie del sol T = 5.77 x 10 3 K 
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5-3) La ecuación fundamental de un gas ideal está dada por: 


U = U„ 




(C+1) 

exp 




So\ 

No] 


(5.19) 


• a) Obtener las funciones S = S(T, V,N),S = S{P, V,N),S = S(T, P , N), 

• b) Obtener las capacidades caloríficas a volumen y presión cons¬ 
tante y probar la relación de Mayer c p — c v = R 


*Método: Combinar entre si las derivaciones de la ecuación fundamental. 


• Para obtener S = S(T, V., N) tenemos en cuenta que por derivación de la ecuación 
fundamental T = (yyt¡) VN obtenemos T = T(S,V,N) de la forma 


T = 


RcN 


U n 




(5.20) 


que teniendo en cuenta que para U = U 0 ,V = V „, N = N a , T = jf- podemos 
simplificar como 


T = T a 




Sólo nos resta despejar S 


S_ 

Ñ 


^r + Rcl n 
N n 



+ R In 




(5.21) 


(5.22) 


• De forma similar, para obtener S(P, V, N ), podemos deducir P = P{S, V, N) por 
derivación de la ecuación fundamental P = — (^) SN 


U n 


1 


cK \V 0 


V 


N 

ñ: 


exp 


1 


5 


Re \N N 0 


que de forma análoga al caso anterior podemos poner como 

c+l __ c+1 

exp 


P = Pnl £- 


Vo 


y despejando S llegamos a: 
5 S, 


N 


1 


s So 


Re \N N c 


— - -f = Re ln — + Rln — — 


N N, 




c+l 


Vo 


N 


N n 


-(c+l)' 


(5.23) 


(5.24) 


(5.25) 
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Para obtener S = S(T, P, N) tenemos en cuenta que al dividir (5.21) y (5.24) se 
tiene 


T _ IV 
P _ RÑ 


(5.26) 


con 

RT„ 

Po 


lo que si sustituimos en la expresión (5.22) y tenemos en cuenta que jf- = 
llegamos a 


S_ 

Ñ 


+ Re ln 
N n 



+ Pin 



(5.27) 


que podemos escribir como 


S_ 

N 


k +E(c+ 1 )ln (B +Rln (T ) 


(5.28) 


• Para obtener las capacidades caloríficas a presión y volumen constante tenemos 
en cuenta que Cy = yC'p=T(|^) p , lo cual podemos obtener a partir 

de las ecuaciones S = S(T , V, N) y S = S(T. P, N) deducidas anteriormente. 

Así en primer lugar derivando S(T, V, N) en (5.22) llegamos a 

Cy= T (^pj v =NR¿ (5.29) 

y derivando S(T,P,N) en (5.28) obtenemos 

Cp=t (ñP) p =nr{c+i) (5 - 3o) 

de donde sigue trivialmente la relación de Mayer c p — c v = R y la siguiente 
expresión para el índice adiabático 


Cp _ c + 1 
Cy c 


(5.31) 


• Con estos resultados las ecuaciones deducidas en apartados anteriores pueden 
ser reescritas de la siguiente forma 


s = s 0 +ln^^^ + Pin 



(5.32) 



(5.33) 

(5.34) 


donde s, v ,c p , c v representan las magnitudes molares entropía, el volumen, y 
capacidades caloríficas respectivamente. 
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5-4) Utilizar los resultados del problema anterior para obtener los poten¬ 
ciales termodinámicos para un gas ideal. 


*Método: Para obtener F = F(T, V, N ) tener en cuenta F = U(T, V, N) — 
TS(T, V, N). Seguir idéntico razonamiento para H y G en función de sus 
respectivas variables naturales. 

• Para obtener la función de Helmholtz F = U — TS reparamos en que de los 
problemas anteriores hemos obtenido 

u = c v T 

s = So + hi^jr^j +Rln(^f^j 

por lo que podemos escribir la función de Helmholtz molar / = u — Ts como 

f(t,v) = c v T-T ^s 0 + \n(^-^j +jRln (¡r)) ( 5 - 35 ) 

Como para T = T a , v = v 0 se tiene / = c v T a — T 0 s 0 , deducimos = c v — s 0 por 
lo que reagrupando términos llegamos a 

/«,„) = /„!- ri„(0''-i i ri»(£) (5.36) 

• Para obtener la función de Gibbs G = G(T, P, N), dada en forma molar por 
g = u — Ts + Pv reparamos en que 

u = c v T 

s = s a + ln (jf- j P + iíln (^) 

Pv = RT 


por lo que podemos escribir 

g(T,P) = c v T-T(s 0 + ln(JQ ” 

que teniendo en cuenta que c v + R = c p ponemos 

g(T,P) = c p T-T(a 0 +]n(^j ” + i?ln0^) (5.38) 

y por un razonamiento similar al caso anterior c p — s a = y por tanto 

g(T,P) = Y o T ~ Tln (yQ P+RT ln(¿) ( 5 - 39 ) 


i?ln(^ 


RT 


(5.37) 
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• Para obtener la entalpia H = H(P, S, N) dada por H = U + PV necesitamos 
expresar U = U(P, S, N) y V = V(S } P,N) sin embargo, podemos acortar el 
procedimiento expresando primero U y PV como función de T y luego expresar 
T = T(S, P, N) 

u = c v T 
Pv = RT 

con lo que la entalpia en función de la temperatura está dada por 

h = u + pv = c v T + RT = c p T (5.40) 


Ahora bien de la expresión obtenida anteriormente S = S(T, P , N ) dada por 


s = s 0 + ln(^pj P +R\n 


Po 

P 


podemos despejar T llegando a 


T = T ° AP exp — (s - s 0 ) 


que sustituida en (5.40) conduce a 


h(s, P) = c p To ( ) exp < — (s — s 0 ) 


(5.41) 


(5.42) 


(5.43) 


• Ejercicio adicional: 

— Contrastar estos resultados con los que podrían haber sido obtenidos utili¬ 
zando las expresiones deducidas de la ecuación de Euler 

F — PV + ¡iN 
H = TS + pN 
G = pN 

y la forma del potencial químico obtenida en el problema del capítulo 3. 


74 


5-5) Obtener la ecuación fundamental en la representación de la función de 
Helmholtz para un gas de fotones de ecuación fundamental energética 

U = aS' 4/3 y" 1/3 (5.44) 

donde a es una constante dada por a = (|) 4 ^ 3 (jj) 1 c es I a velocidad 
de la luz y a la constante de Stefan Boltzmann. 


*Método: Obtener T = T(S,V) por derivación deU = U(S,V ) y eliminar 
S y V entre estas dos ecuaciones y la definición de Función de Helmholtz 


F = U-TS. 


• Si derivamos (5.44) con respecto a la entropía obtenemos 

T = ^aS 11 / 3 ^- 1 / 3 (5.45) 

O 

por lo que nuestro problema se reduce a eliminar S,U de las ecuaciones 

U = aS^V- 1 / 3 
T = faS* 1 / 3 !/- 1 / 3 
F=U -TS 


Esto podemos hacerlo despejando S = S(T, V) en (5.45) lo que nos da 


S = 



(5.46) 


• Para obtener U = U(T, V) sustituimos S = S(T, V) en la ecuación fundamental 
lo que nos lleva al resultado ya obtenido en problemas anteriores 


U = 



(5.47) 


• La función de Helmholtz se obtiene sustituyendo F(T, V) = U ( T , V) — TS(T, V) 

F =U -TS = ^ a~ 3 VT A - T ^ a" 3 7T 3 -| ^ a~ 3 VT 4 (5.48) 

que podemos poner teniendo en cuenta el valor de a como 

F=--VT a (5.49) 

3 c 
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Ecuaciones térmicas de estado. 


• Las relaciones que se deducen de los potenciales termodinámicos y las relaciones del 
cálculo diferencial permiten expresar distintos coeficientes termodinámicos en términos 
de las capacidades caloríficas y coeficientes deducibles de la ecuación térmica de estado. 
La ecuación de estado de un fluido puede ser determinada de diferentes maneras. A 
baja presión la mayoría de los gases aparecen bien descritos por la ecuación del gas 
ideal 

Pv = RT (5.50) 

donde v es el volumen molar. 

• La ecuación de Van der Waals describe aparte de las propiedades próximas al compor¬ 
tamiento de gas ideal, la transición de gas a líquido y la existencia del punto crítico 

(p + ( v ~ b ) = RT (5.51) 


Las constantes o y b son características de cada sustancia, si bien pueden referenciarse 
a las magnitudes del punto crítico de la forma: 


27 R 2 T 2 

64 P c 

(5.52) 

, KT C 

i= w r 

(5.53) 


• Una ecuación de estado de particular interés es la denominada desarrollo del virial 
la cual expresa el producto de la presión por el volumen molar en términos de un 
desarrollo en serie de potencias de 4 o de P. 


o bien 


Pv = RT( 1 + 


MT) 

V 



+ ....) 


Pv = RT{ 1 + B' 2 {T)P + B' 3 (T)P 2 + ....) 


(5.54) 

(5.55) 


Los coeficientes B¿(T) o B'\T) son funciones de la temperatura específicas de cada 
sustancia y pueden determinarse experimentalmente o bien a partir de modelos mole¬ 
culares. 

Con frecuencia estos desarrollos pueden truncarse a partir del primer orden dando lugar 
a ecuaciones del tipo 

Pv = RT{ 1 + ^11) ( 5 . 5 6 ) 

v 

Pv = RT(Í + B' 2 (T)P) (5.57) 
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• Puede probarse que: 


— Los coeficientes del virial B- 2 (T) } B%(T) pueden obtenerse a partir de los B^T), 
B' 3 (T) y viceversa. (Despejar P en el primer desarrollo y sustituir en el segundo, 
e igualar ambos desarrollos) 

— La ecuación de Van der Waals puede ponerse como un caso particular de los 
desarrollos del virial (Escribirla de la forma P = ^ y desarrollar el 

término 1 b ) 

• Existen una amplia variedad de ecuaciones térmicas para fluidos algunas menos usadas 
como las de Berthelot, Clausius y Dietireci, y otras de uso especializado como las de 
Redlich Kwong, Beattie-Bridgeman, Benedict-Webb-Rubin,etc. 


• Una goma elástica considerada como un sistema mono dimensional con la elongación 
L — L 0 (donde L es la longitud y L a es la longitud en ausencia de tensión), como 
coordenada de trabajo y la tensión r como variable conjugada puede describirse en 
términos de la ecuación térmica: 


r = cT 




(5.58) 


• Una interfase (superficie de separación entre dos fases) en la cual el área es la coorde¬ 
nada de trabajo y el coeficiente de tensión superficial es la variable conjugada puede 
describirse en términos de la ecuación térmica de Eótvos 

a = v 2 ' 3 k(T c - T) (5.59) 

donde v es el volumen molar del fluido, T c es la temperatura crítica y k una cons¬ 
tante característica del mismo. Para el agua k = 1.03 er g/K] T c = 374 K\ v = 
18 cm 3 /mol 



La consideración de las fuerzas entre las moléculas que forman un fluido llevaron a Johannes van de Waals (1837- 
1923) a proponer su ecuación de estado. En la figura, de pie junto a Heike Kamnrerling Onnes (1853-1926) prime¬ 
ro en la licuefacción del helio y descubridor de la superconductividad. 
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5-6) Joule y Gay Lussac comprobaron en respectivos experimentos que en 
una expansión expansión adiabática hacia un volumen inicialmente 
vacío (expansión libre) no se produce variación de la temperatura. 
Sin embargo este resultado sólo es exacto si el gas puede considerarse 
ideal. Evaluar la variación de temperatura que se produciría si el 
comportamiento del gas fuera el de un gas de Van der Waals. 

(p + (v - b) = RT (5.60) 

Analizar cual sería la variación de temperatura si 1 mol de nitrógeno 
(a = 0.136Pa x m 6 moI” 2 ; b = 38.5 x 10 _6 m 3 moZ _1 ) a 25°C en un recipiente 
de 5 litros duplica su volumen sobre un recipiente vacío. Supóngase 

c v = | R. 


*Método: Expresar la temperatura en términos de las variables conocidas 
a lo largo del proceso, en este caso volumen y energía interna (permane¬ 
ce constante). Reducir la correspondiente derivada de la temperatura a 
magnitudes calculables a partir de la ecuación térmica. 


• El problema se reduce a analizar la variación de temperatura en un proceso en 
el que se produce una cierta variación de volumen. Como el gas se expande 
adiabáticamente hacia un volumen vacío la variación de energía interna es cero 
y por tanto el proceso tiene lugar a energía interna constante. 

• Podemos pues plantear el problema tomando la temperatura como función del 
volumen 

T = T(V ,...) (5.61) 

la segunda variable de la que tenemos información es la energía interna, que 
permanece constante por lo que consideramos 


T = T(V, U) 

En una etapa infinitesimal del proceso se tiene 

dT 8T 

dT=( w ) u dV + ( w ) v dU 


(5.62) 


(5.63) 


y como éste tiene lugar a energía interna constante (dU = 0) esta expresión se 
reduce a: 

dT = (^)udV (5.64) 
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• Podemos pues estudiar el proceso analizando el coeficiente (fpO¡y que tratare¬ 
mos de reducirlo a magnitudes que puedan determinarse a partir de la ecuación 
térmica y a las capacidades caloríficas. 

• Como la derivada (§f)u involucra un potencial termodinámico que es la energía 
interna, la expresamos en términos de derivadas de la misma: 


W ( f v)t 

0V )U (§%) V 


(5.65) 


• Como V y T no son las variables naturales de la energía interna para evaluar las 
derivadas (fpOr y {%r)v es necesario considerar U(T,V) de la siguiente forma: 



U = U(S,V) = U{S{T,V),V) 

(5.66) 

por tanto 

.dU. .du. .OS. .dU. 

{ W )T ^ { dS )v{ dV )T + { W )s 

(5.67) 

y 

.du. BU. f dS. 

( ar )v - ^dS' v ^dT’ v 

(5.68) 

de donde 

(—)t = T{—) t - P 
y dV >T y dV ,T 

(5.69) 


( dU ) T( dS ) 

dT V T dT* v 

(5.70) 

En el segundo de los 

casos podemos identificar la derivada con 

la capacidad 

calorífica a volumen constante 



( du ) T( dS 
y ~df’ V ~ 1 y dT )v ~ L ' v 

(5.71) 


y en el primer caso debemos obtener 

• Para determinar esta derivada reparamos en que V y T son las variables naturales 
de la función de Helmholtz y como S = — (§fr)y corresponde por tanto a una 
de las relaciones de Maxwell 


/dS\ ( d 2 F \ / d 2 F \ ídP\ 

\dV) T - ~ [dVdfjy - ~ V dTdVjy ~ \df) V 

con lo que el coeficiente analizado puede escribirse como 

( dT T{%) v -P 

{ dV )u= CV 


(5.72) 


(5.73) 


donde el coeficiente 


(wt) v' P ue de determinarse a partir de la ecuación térmica. 
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El resultado del problema puede obtenerse integrando la ecuación 


dT = — 


T (dE) _ p 
1 KdTJV r 


Cy 


-dV 


(5.74) 


sin embargo, a efectos de poder integrar, las magnitudes (f^) y y P deben ser 
expresadas en función de T y V. 


Así pues ponemos 

con lo que 

por lo que 


T 


dP 

OT 

dP 

OT 


v 


RT 
+ (v - b) 

(5.75) 

R 


(v-b) 

(5.76) 

~ P = ~2 

(5.77) 


por tanto la variación de temperatura (5.74) podemos ponerla como 

a 


dT - 


-dv 


(5.78) 


Integrando entre un estado inicial [T 0 ,v 0 ) y un estado final (T,v) llegamos a 


T-T 0 = 

Si v = 2v a entonces tenemos 


f 

J V n 


-dv 


ai 1 1 


V Vr 


_ al 

T — t o = _ 

Cy 2 Vq 


(5.79) 


(5.80) 


Vemos que este resultado reproduce el comportamiento para un gas ideal donde 
a = 0 y por lo tanto el enfriamiento es nulo en este caso. 

Analicemos este resultado numéricamente para el nitrógeno 

a = 0.136 Pam e mol~ 2 

c v = | R = |8.314 JmoZ _1 A' _1 

v„ = 5 / mol ” 1 x 10~ 3 m 3 Z _1 = 5 x 10 -3 m 3 mol -1 

Por tanto 

0.136 N m~ 2 m 6 mol~ 2 1 


T - T n = - 


§8.314 Nmmol- 1 R- 1 2 x 5 x 10- 3 m 3 mo¿- 1 


= -0. 65 K 








5-7) La atmósfera es un sistema complejo y en los procesos que tienen lu¬ 
gar en la misma intervienen múltiples efectos no siempre sencillos de 
analizar. Uno de los factores que intervienen en las partes bajas de 
la misma son las variaciones de temperatura asociadas al movimien¬ 
to de convección vertical de las masas de aire. El proceso se puede 
modelizar como un proceso de expansión adiabática en el que una 
masa de aire se expande desde la presión a nivel de superficie hasta 
la presión a una determinada altitud la cual varía según la expresión 
de la hidrostática dP = —pgdz . Supuesto comportamiento para el aire 
de gas ideal de calor específico a presión constante c. p = 1004 J kg~ 1 K^ 1 
obtener la variación de la temperatura con la altitud debida a este 
fenómeno. 


*Método: Considerar la relación entre la variación de entropía, tempera¬ 
tura y presión y tener en cuenta que en un proceso adiabático la variación 
de entropía es nula. 


• Una expansión adiabática implica una variación nula de la entropía dS = 0. La 
variación de la entropía en función de la temperatura y la presión está dada por 

' ¡s =(I)/ t+ (I) t ' ¡p (58i) 

ahora bien (ff ) P = % y por las relaciones de Maxwell (§p) T = — (§r) n 

• Para el sistema que nos ocupa asumiendo comportamiento ideal — (§^) p = 
— ^ o lo que es lo mismo — (§^) p = P or t an to sustituyendo en la expresión 
(5.81) llegamos a 

Y dT + J¡P9 dz=0 ( 5 - 82 ) 

teniendo en cuenta que rn = pV podemos poner 

C p dT + mgdz= 0 (5.83) 


y dividiendo por la masa 


ÓT 

dz 


(5.84) 


donde c p representa la capacidad calorífica por unidad de masa o calor específico. 


• Numéricamente g = 9.8 ms “y c p = 1004 JA' 1 kg 1 por lo que 


dT 

dz 


9.8ms 2 


1004 J K~ 1 kg 1 


9.8 x 10~ 3 A'm 


-i 


(5.85) 


es decir la temperatura como consecuencia del movimiento vertical de convección 
desciende aproximadamente 1AT cada 100 metros. 
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5-8) Newton calculó una expresión para la velocidad del sonido interpre¬ 
tando que las expansiones y compresiones que tienen lugar al paso 
del sonido por un gas son isotermas y por lo tanto el gas está sujeto 
a la Ley de Boyle. Sin embargo los valores calculados subestimaban 
apreciablemente los resultados experimentales. Laplace revisó la ex¬ 
presión de Newton interpretando que las expansiones y compresiones 
son adiabáticas obteniendo una expresión más acorde con los resulta¬ 
dos experimentales. En general, la interpretación de Laplace conduce 
a que la velocidad de propagación del sonido por un fluido puede po¬ 
nerse de la forma 



(5.86) 


donde p es la densidad del fluido. 

Obtener una expresión para la velocidad del sonido a partir de la 
expresión anterior asumiendo que el comportamiento del gas es ideal. 


*Método: Expresarla como una derivada respecto del volumen y reducir 
la derivada resultante a expresiones calculables a partir de las capacidades 
caloríficas y la ecuación térmica. 

• Escribimos 



(5.87) 


• Por una parte 



(5.88) 


ahora bien 




(5.89) 


• Además podemos ver que 



(5.90) 


por tanto 



(5.91) 
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• Por otra parte la densidad p = y podemos ponerla P = = y donde v es el 

volumen molar y M la masa molecular por lo que 


(dv\ __M 
\®Pj s~ P 2 


que podemos escribir como 



• Así pues 


Cp (dP\ v 2 _ 
C v { dv ) T M 


• En términos del coeficiente de compresibilidad isotérmica Kr 
del índice adiabático 7 = y- podemos escribir 

2 1 v 

° _ 7 AtM 


(5.92) 


(5.93) 


(5.94) 



(5.95) 


• En el caso de un gas ideal (§y) T = con lo que 



(5.96) 


• Ejercicio adicional: 

— Comparar los resultados para el argón (Ma t = 40 gmol^ 1 ) y para el aire 
(21%0 2 , 79%N 2i Mq 2 = 32gmol- 1 ,M N2 = 28 gmor 1 ). 


Newton publicó sus resultados en la primera edición de los Principia (Philosopiae 
Naturalis Principia Mathematica) en 1687. El valor calculado por Newton fué de 
295m S -1 (968/ís -1 ) que comparó con los resultados de 183 ms ” 1 obtenido 
por Roberval, y 450?ns _1 obtenido por Mersenne y con sus propias medidas 
que oscilaban entre 280 y 330 ms" 1 . Al tiempo de la segunda edición de los 
Principia en 1713 Halley y Flamsteed habían obtenido un valor experimental 
más preciso de 348TOS -1 un 17% mayor que el valor calculado. Newton rehizo 
sus cálculos incorporando valores más recientes de las constantes atmosféricas 
obteniendo 298 m S -1 y atribuyó la discrepancia al hecho de que el cálculo sólo 
tiene en cuenta la velocidad del sonido en el espacio entre las partículas de aire 
y no a través de las propias partículas, así como a la presencia de otros vapores 
en el aire los cuales disminuyen su densidad por lo que el valor experimental, 
según Newton, debía ser mayor que el calculado. Ambas hipótesis quedaron 
descartadas cuando P. S. Laplace en 1816 (Ann. Chim. Phys. 3,238 (1816)) 
incorporó el factor multiplicativo del índice adiabático. 


83 




5-9) La ecuación térmica de estado de una goma elástica (elastómetro ideal) 
de longitud L a la temperatura T y sometida a la tensión t está dada 
por 



(5.97) 


Probar que: 


• i) A tensión constante, la goma se contrae ante un aumento de 
temperatura. 

• ii) A temperatura constante la entropía aumenta con la longitud 

• ii) La goma se calienta cuando se estira adiabáticamente. 

• iii) A temperatura constante, la energía interna no depende de la 
longitud. 


*Método: Utilizar las relaciones entre las derivadas parciales de L(T,t ), 
S(T,L),T(L,S) 


• i) La variación de la longitud con la temperatura manteniendo constante la 
tensión viene determinada por L(T,t) 


dL = 




(5.98) 


pero al ser r = cte el término dr = 0 con lo que 


dL = 



(5.99) 


Ahora bien 

(»l\ (Mi, 

V®7 r (ft) T ¿r(¿ + 2g) 


l (£ - ¿j) 

T(L‘+2L'¿) 


(5.100) 


Por tanto, si dT > 0 entonces dL < 0, es decir la goma se contrae si se aumenta 
la temperatura. 




ii) La variación de entropía con la longitud a temperatura constante está carac¬ 
terizada por S = S(T,L) con lo que con un razonamiento análogo al apartado 
anterior llegamos a 


dS = 



dL 
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Para calcular (ff ) T tenemos en cuenta que S = — (ff y P or tanto usando 
las relaciones de Maxwell 


dS\ _ í d 2 F \ 
~dL) T -- [dLdf) 

con lo que 


f d 2 F \ 
[dfdLj 




LoL* 


< 0 


(5.102) 


(5.103) 


es decir la entropía disminuye con la elongación. 


• iii) La variación de temperatura con la elongación en un proceso adiabático 
podemos caracterizarla por T = T(L , S) con lo que 


dT - 


dT 

~dL 


dL 


ahora bien las relaciones entre derivadas parciales nos llevan a 

(M) 

\dLjT 


dT\ 

~dí) : 


(SS) 

\dT) i 


Para calcular (§j-) L tenemos en cuenta que : 


dS 

dT 


Cl 

T 


(5.104) 


(5.105) 


(5.106) 


por lo que utilizando el valor de (ff ) T calculado en el apartado anterior llegamos 


a: 


dT 

dL 




Cl 

T 


(L 3 -Ll) 


T 


L 0 L 2 C l 


> 0 


(5.107) 


es decir, en una elongación adiabática reversible (a entropía constante), la tem¬ 
peratura aumenta con la elongación. 


• iv) La dependencia de la energía interna con la temperatura y elongación po¬ 
demos analizarla a partir de U — U(S, L) = U{S{T, L),L) de donde deducimos 
que la variación con la elongación a temperatura constante 


ídU\ _ fdU\ fdS\ fdU\ 
{dL) T -{dSj L {dL) T + {dLj s 


(5.108) 


con lo que 



(5.109) 


Es decir la energía interna de un elastómero ideal es función únicamente de la 
temperatura. 
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5-10) J. P. Joule y W. Thomson (Lord Kelvin) realizaron un experimento en 
el que un fluido a alta presión es conducido adiabáticamente a través 
de un tabique poroso a una región de más baja presión. La variación 
de energía interna del fluido puede establecerse como «2 — íq = P 2 V 2 — 
P\Vi o lo que es lo mismo ¡q + P 2 V 2 = tq + PiV\ donde Uí,Ví,P¿ son la 
energía interna y volumen molar y la presión del gas antes y después 
de atravesar dicho tabique. Esto indica que dicho proceso, conocido 
como expansión Joule-Kelvin, es un proceso que tiene lugar a entalpia 
(h = u + pv ) constante. 

Analizar como es la variación de la temperatura con la presión en una 
expansión Joule Kelvin para un gas de Van der Waals. 

(p + ( v ~ b ) = RT ( 5 ' n °) 

Analizar para que estados no se produce variación de temperatura 
(temperatura de inversión) y particularizar los resultados al caso del 
gas ideal. 


*Método: Relacionar la temperatura con la presión y la entalpia y reducir la 
derivada correspondiente a magnitudes deducibles de la ecuación térmica. 


• El proceso a analizar es la variación de temperatura con la presión en un proceso 
a entalpia constante. Por tanto será necesario poner 

T = T(P, ..) (5.111) 


como la segunda variable que se conoce es la entalpia, la cual permanece cons¬ 
tante ponemos: 

T = T(P,H) (5.112) 


• E 11 un proceso infinitesimal se tendrá 


dT = 



dP + 

H 



que como el proceso es isoentálpico (dH = 0) se tiene 


(5.113) 


iT= (^)„ iP <5114) 

ecuación que deberá ser integrada para analizar el comportamiento en una etapa 
finita. 


• Por tanto el coeficiente a analizar es (§p) H que trataremos de expresar en 
términos de magnitudes que puedan derivarse de la ecuación térmica y de las 
capacidades caloríficas. 


86 






En primer lugar expresaremos esta derivada en función de derivadas de la en¬ 
talpia: 


fdT\ (gf) r 

\ dp ) H m) P 


(5.115) 


Como T y P no son las variables naturales de la entalpia, ponemos ésta de la 
forma 


por lo que 



que expresamos como 


H(T,P)=H(S(T,P),P) 




(5.116) 

(5.117) 

(5.118) 


• Observamos que en la derivada (§p) T , la presión P y la temperatura T son las 
variables naturales de la función de Gibbs, y S es una de sus derivadas primeras, 
por lo que podemos aplicar una de las relaciones de Maxwell 


(dS\ /3 2 G\ / d 2 G \ fdV\ 

\dPj T ~ \dPdTj- \dTdP)~ \dT) p 


(5.119) 


con lo que la derivada buscada queda de la forma 


= _r(S + v 


dP) T \dTy p 


(5.120) 


• Por otra parte la derivada 


d H 
~dT 


dH\ fdS\ 


dS J p \dTjp 


(5.121) 


que podemos poner como 


dH 


dTjp \dTyp 


T(1 


(5.122) 


la cual corresponde a la capacidad calorífica a presión constante 




(5.123) 


• Así pues el coeficiente buscado podemos escribirlo como 


8T 

dP 


H 


- T (w)p + v 

Cp 


(5.124) 
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• Para el caso particular del gas de Van der Waals 


(P + ~^)(v - b) = RT 


(5.125) 


la derivada (f^) P es complicada de evaluar, pero sin embargo podemos obtener 
(§ v) p = (av \ < l ue deducimos de la ecuación térmica como 

\ar )p 


dT 


1 , 2a 


~ “ b ) + d( P + - 2 ) 


R v 3 


R' 


que teniendo en cuenta de (5.125) que 


queda 


por lo que 


de donde 


— (P+ —) = ——— 

P v 2 (v — b) 


dT\ !/ 2o,. 

-I = í ¡(- í j)<»-'>) + 


dv 


(v~b) 




M-§s)(” - » + iA, 


— v 


dT\ v f TRv 2 b - 2a(v - b) 2 \ 

dP] c p \ -2 a(v - b) 2 + TRv 3 ) 


(5.126) 

(5.127) 

(5.128) 

(5.129) 

(5.130) 


Vemos por tanto que este coeficiente depende de la temperatura y del volumen 
molar (es decir de la densidad). 

• Asimismo vemos que en el caso particular de un gas ideal (a = 0, b = 0),este 
coeficiente es nulo y por lo tanto no se produce variación de temperatura. 


• La temperatura de inversión corresponde a aquella en la que el coeficiente de 
Joule-Kelvin cambia de signo (§j¡) lf = 0 por lo que de (5.130) podemos ver que 
está dada por la solución de la ecuación 


TR.v 2 b - 2a(v - fe) 2 = 0 

que conduce a 

2a(v - b) 2 
R,v 2 b 


(5.131) 

(5.132) 


• Como en condiciones ordinarias el volumen molar v » b el término 

(y - b ) 2 _ , 

v 2 ~ 

por lo que podemos escribir 


que nos ofrece una buena estimación de la temperatura de inversión. 


(5.133) 

(5.134) 
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• Ejercicio adicional: Estimar los valores de la temperatura de inversión para 
el Nitrógeno (a = 0.136 Pam e mol~ 2 ; b = 38.5 x 10 _6 m 3 mo¿~ 1 ), el Hidrógeno 
(a = 0.0248 Pa m e mol ~ 2 ; b = 26.6 x 10 _6 m 3 mo/^ 1 ) y el Helio (a = 0.00346 Pa 
m e mol~ 2 ; b = 23.7 x 10 ~ e m 3 mol~ 1 ). 


Joule y Thomson publicaron el experimento que lleva su nombre en 1853 en 
un trabajo titulado "Sobre los efectos térmicos de los fluidos en movimiento” 
publicado en las ” Philosophical transactions of the Royal Society (London)” 143 
(1853) 357-365 y continuaron sobre el mismo tema en un segundo trabajo con 
el mismo título aparecido en dicha revista en 1854. Su interés inicial estaba 
centrado en la determinación de la función de Carnot. Esto es, la función que 
relaciona el cociente entre el trabajo realizado por un motor de Carnot y el calor 
absorbido por el mismo con la diferencia de temperaturas entre las que opera. 

El efecto Joule Thomson fue utilizado con fines prácticos por Cari von Linde 
para licuar el aire e incorporado posteriormente a los diseños de los refrigeradores 
comerciales. 



James Prescott Joule (1818-1899), hijo de un cervecero de Manchester comenzó investigando el efecto de un pri¬ 
mitivo motor eléctrico, para, después de una larga y meticulosa serie de experimentos llegar a la conclusión de que 
el calor y el movimiento (trabajo) son interconvertibles siempre en la misma proporción, resultado decisivo para el 
establecimiento del principio de conservación de la energía. 
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5-11) Las variaciones de las magnitudes termodinámicas sobre la curva de 
equilibrio de fases pueden analizarse teniendo en cuenta que sobre 
dicha curva la presión es únicamente función de la temperatura. Ana¬ 
lizar la variación del calor latente con la temperatura. 


*Método: Expresar el calor latente en términos de la entalpia y ésta como 
función de la temperatura y la presión H(T,P(T)) teniendo en cuenta que 
P(T) representa la dependencia entre estas magnitudes sobre la curva de 
equilibrio. 


• El calor latente está dado como la diferencia de entalpias molares entre ambas 
fases 

¿12 = h, 2 - h i (5.135) 

por tanto evaluar la variación del calor latente sobre la curva de saturación es 
equivalente a evaluar la variación de la entalpia sobre dicha curva. Hemos de 
tener en cuenta que sobre la curva de equilibrio varían tanto la presión como la 
temperatura, pero que ambas están ligadas por la relación P = P(T ) proveniente 
de integrar la ecuación de Clausius-Clapeyron 


dP\ 1 12 

dTj sat T(vW-vW) 


(5.136) 


donde el subíndice () sat denota que la derivación tiene lugar a lo largo de la 
curva de equilibrio de fases. 


• Consideremos por tanto h(T, P(T)), en estas condiciones 


/ dh \ / dh\ / dh \ / dP\ 


(5.137) 


ahora bien temperatura y presión no son las variables naturales de la entalpia, 
por tanto debemos poner 


h = h(T,P) = h(s(T,P),P ) 


(5.138) 


para expresar ésta en términos de sus variables naturales entropía y presión. Las 
correspondientes derivadas están dadas por 


dh\ 

dTj, 


dh 

ds 


ds\ 

dTj, t 


T ( — 
1 dT 


dP 


dS 


dP 


dP J , 


p 

Sil =ÍSÍ1 ÍÜl +ÍSÍ1 =-rfih 


con lo que sustituyendo en (5.137) llegamos a 


dh 

dT 


— Cri 


-T 


dv 

dT 


dP 

dT 


(5.139) 

(5.140) 

(5.141) 
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• Por tanto la variación del calor latente podemos ponerla haciendo uso de (5.135) 
y (5.136) como 


(di 12 \ 

\dTj sat 


- r (2)_ r (l)_ 




\\ ¿12 

ar t(u( 2 )-u(P) 


• En el caso de la fusión el coeficiente de dilatación (-§f) p es pequeño tanto para 
el líquido como para el sólido por lo que puede despreciarse conduciendo a 


(di i 2 \ 

\dTj sat 





que podemos reescribir como 


■¡(VA cf-c» 1 » 

dT } T 

/ sai 


(5.142) 


(5.143) 


• En el caso de la sublimación o vaporización el volumen molar y el coeficiente de 
dilatación de la fase condensada (líquido o sólido) puede despreciarse frente al 
del vapor y asumiendo para éste comportamiento ideal se tiene 


/ dv \ R v 

[df) p = ~P = T 


con lo que el segundo sumando es nulo y por tanto 


(dl 12 \ 

\ dT J sat 


- r (2) _ r (l) 

'-'P 


(5.144) 


(5.145) 


• Ejercicio adicional. 

— En el entorno del punto triple (T = 273.15 K) los calores específicos del 
agua líquida y del vapor de agua están dados por Cp' 1 = 4200 Jkg^ 1 K^ 1 y 
¿p' > = 1850 J kg~ 1 K~ 1 respectivamente y el calor latente de vaporización 
l v = 2.5 x 10 6 J kg -1 calcular una expresión para el calor latente en el 
entorno del punto triple. Utilizando dicha expresión y teniendo en cuenta 
que la presión en el punto triple es de 611.73 Pa , obtener una expresión 
para la dependencia entre la presión y la temperatura sobre la curva de 
saturación. 
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5-12) En un diagrama P-T para una sustancia pura en el punto triple se cor¬ 
tan las curvas de equilibrio sólido-vapor (sublimación), líquido-vapor 
(vaporización), y sólido-líquido (fusión). ¿Que es mayor, la pendiente 
de la curva de vaporización o la de sublimación?. 


* Método: Utilizar la ecuación de Clausius-Clapeyron y discutir la relación 
entre los calores latentes interpretándolos como variaciones de la entalpia. 


• En el punto triple la pendiente de la curva de sublimación está dada por 


(dP\ = l s 

\dTj s _ v T(«W- C W) 

mientras que la de la curva de vaporización 

(dP\ = ly 

\dTj T(«M-«(')) 


(5.146) 


(5.147) 


donde l v es el calor latente de vaporización y l v el calor latente de sublimación 
y T la temperatura del punto triple. 

• El volumen molar del líquido y el del sólido v^ son pequeños frente al 
volumen molar del vapor v^ v ' que es el mismo en ambos casos por lo que el 
cociente entre ambas pendientes está dado por 



ls 


(5.148) 


• Ahora bien, el calor latente corresponde a la diferencia entre las entalpias molares 
entre ambas fases, por lo que en el caso de la sublimación 


l s = h {v) - h {s) 

(5.149) 

que podemos poner como 

l„ = ñ (í,) - h,® + h,® - /i (s) 

(5.150) 

de donde deducimos que en el punto triple 

l s = lf ly 

(5.151) 

y por tanto al ser positivas las tres magnitudes 

l s ^ ly 

(5.152) 


con lo que en el punto triple la pendiente de la curva de sublimación es mayor 
que la de vaporización 
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Capítulo 6 


Sistemas magnéticos. 
Enfriamiento magnético. 

6.1 Objetivos 

• Comprender el uso de las variables termodinámicas en los sistemas magnéticos. 

• Analizar algunos procesos sencillos como el enfriamiento magnético. 

El enfriamiento magnético es una técnica para obtener bajas temperaturas 
que consiste en crear un campo magnético a temperatura constante en tor¬ 
no a una sustancia y hacerlo desaparecer mientras se mantiene la muestra 
aislada adiabáticamente. Este efecto suele denominarse efecto magneto- 
calórico. Dos factores importantes son el descenso de temperatura que 
es posible conseguir con un determinado campo y el calor liberado por la 
sustancia cuando se magnetiza a temperatura constante. 



Walther Meissner (1882-1974) y Robert Ochsenfeld (1901-1993) descubrieron en 1933 que un material supercon¬ 
ductor es un diamagnético perfecto, es decir que el campo se anula completamente en el interior del mismo. Así un 
imán colocado sobre un material superconductor levita sobre éste. 
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6.2 


Problemas. 


6-1) La ley de Curie establece que la susceptibilidad magnética es propor¬ 
cional al inverso de la temperatura lo cual en términos de la magne¬ 
tización total y la inducción magnética podemos ponerlo como 

M =^Vpf ] B G ( 6 . 1 ) 

que corresponde a una ecuación térmica para un sistema magnético. 
Calcular la variación de temperatura que experimenta una sal para¬ 
magnética que obedece a la ley de Curie cuando sufre una desimana¬ 
ción adiabática desde un campo B a hasta B a = 0. 


*Método: Considerar la temperatura en función de la inducción magnética 
y la entropía y aplicar las técnicas habituales de derivación. 


• El enfriamiento podemos calcularlo describiendo la temperatura en función de 
la inducción magnética en el vacío B a y la entropía 


T = T(B 0 ,S) 


( 6 . 2 ) 


• El proceso correspondiente es un proceso isoentrópico dS = 0 en el que la induc¬ 
ción magnética varía desde un valor B 0 hasta B a = 0 . La variación infinitesimal 
de temperatura 

ír= {m) s iB - + {^) B . iS (63) 

queda reducida a 

* r = (§k) s ÍB ° (64) 

• Atendiendo a las propiedades de las derivadas parciales podemos escribir 


fdr_\ {mk) T 

\° B o)s ® flo 


(6.5) 


El denominador podemos interpretarlo en función de la capacidad calorífica Cb 0 
a campo magnético constante (relación entre el calor absorbido y el aumento de 
temperatura cuando el campo es constante) como 


(dS\ _ C Bo 
\ dT ) Bo T 


( 6 . 6 ) 


mientras que el numerador podemos reducirlo utilizando las relaciones de Max¬ 
well. 
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• En este caso B a corresponde a la variable conjugada de la coordenada de trabajo 
magnetización M . La función de Gibbs definida corno 


G = U -TS- B 0 M 
tiene como variables naturales 

G = G(T, B a ) 

y verifica que 

dS\ ( d 2 G \ ( d 2 G 

i)B„ 


\dB 0 dT J \dTdB 0 


Por tanto podemos poner 


tdM\ 

d T= - dB„ 


c Bn 

T 


Atendiendo a la ley de Curie 


M J^Vi.CB 0 


por lo que 


con lo que 


(3M 


\dT J T 2 

d T- CVB ■ 


^V,j- l B 0 


T /i 0 Cb { 


-dB n 


• Ejercicio adicional: 



(6.7) 

( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 


— Integrar la expresión anterior teniendo en cuenta que, en las proximidades del cero 
absoluto la capacidad calorífica Cb 0 depende de la temperatura de manera que 
en la integración debe incluirse dicha dependencia. Las sales paramagnéticas a 
bajas temperaturas presentan un comportamiento anómalo denominado anomalía 
de Schottky según el cual la capacidad calorífica a campo magnético nulo está 
dada por 


c Bo (T, 0 ) 


Vb 

y2 


(6.14) 


si bien a más baja temperatura los sólidos siguen la ley de Debye 


C Bo (T, 0) = cT 3 


(6.15) 


La obtención de la dependencia de la capacidad calorífica con la magnetización 
Cb„{T, A4) puede llevarse a cabo utilizando la expresión: 


dC Bo \ fd 2 M\ 

dB 0 ) T {dT>) B¡ 


(6.16) 


y utlizando la ley de Curie para la dependencia de la magnetización con la tem¬ 
peratura. 
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6-2) Determinar el calor liberado por una muestra de una sal paramagnética 
que obedece a la ley de Curie cuando sufre una imanación isoterma 
desde un campo B a = 0 hasta un campo B a . 


*Método: Calcular el calor en términos de la variación de entropía y ex¬ 
presar ésta en función de la temperatura y la inducción magnética. 


• El calor liberado lo podemos calcular en términos de la variación de entropía en 
dicho proceso. Para ello consideramos la entropía como función de la tempera¬ 
tura y la inducción magnética en el vacío: 

S = {B 01 T ) (6.17) 

En un proceso infinitesimal la variación de entropía está dada por: 

dS= (M.) T dB - + {w) B ^ ^ 


que en un proceso isotérmico (dT = 0) se reduce a 


dS = 


dS 

dB~o 


dB a 

T 


(6.19) 


• Ahora bien, utilizando las relaciones de Maxwell 


OS \ fdM\ 
dB 0 ) T ~{dTj Bi 


( 6 . 20 ) 


con lo que 


dS = 



dB 0 

B a 


( 6 . 21 ) 


que al igual que en el problema anterior en el caso en que es aplicable la ley de 
Curie 



C_ 

y2 


V,JC l B n 


( 6 . 22 ) 


y por tanto 


dS 


C_ 

" y2 


V Po 1 B 0 dB 0 


(6.23) 


• El calor absorbido en dicho proceso isotérmico estará dado por: 

Q mT AS=-^Vp~ 1 ^ (6.24) 

es decir, el calor es inversamente proporcional a la temperatura y proporcional 
al cuadrado de la intensidad del campo magnético aplicado. 
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Capítulo 7 


Sistemas pluricomponentes 


7.1 Objetivos 

• Comprender los conceptos de presión de saturación, presión parcial y humedad 
relativa. 

• Familiarizarse con el uso de las leyes experimentales deducidas a partir de las 
condiciones de equilibrio en sistemas pluricomponentes. Leyes de Raoult y de 
Henry, descenso crioscópico, aumento ebulloscópico y presión osmótica. 



Francois Mane Raoult (1830-1901), químico francés, profesor de la Universidad de Grenoble que realizó un 
exhaustivo estudio de las propiedades fisicoquímicas de las disoluciones, estableciendo leyes sobre la presión de 
vapor, el aumento del punto de ebullición, y el descenso del punto de congelación. 
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7.2 


Problemas 


7-1) ¿Qué es más denso el aire húmedo o el seco?. Calcular la densidad del 
aire a 20°C y 1 atm de presión con un 90% de humedad y compararla 
con la correspondiente a aire seco a la misma temperatura. La presión 
de saturación del agua a 20°C, P s = 17.5 mmHg. 


*Método: Calcular las densidades teniendo en cuenta la proporción en 
número de moles de aire y nitrógeno y vapor de agua. 


• Para responder a esta pregunta debemos comparar la masa de un volumen de 
aire seco con la de un mismo volumen de aire húmedo en idénticas condiciones 
de temperatura y presión. Hemos de matizar que entendemos por aire húmedo 
aire conteniendo vapor de agua, pero completamente éste en la fase de vapor. 

Si el aire está húmedo la densidad sería 

(h) _ ™ _ Nq 2 Mq 2 + jVjv 2 Mjy 2 + Nh 2 qMh 2 o 

P ~ V ~ V 1 j 

que en términos de las respectivas fracciones molares x P ' 1 podemos poner como 

(x^Mo, + x^¡M N2 + x ( £ 0 M H20 ) N 
P {h) = - - y -— ( 7 - 2 ) 

y haciendo uso de la aproximación de gas ideal 

P {h) = (x%M 02 + x%¡M N 2 + x%l 0 M H2 o) ( 7 . 3 ) 


• Suponiendo que el aire seco se encuentra constituido por un 21% de O 2 y un 
79% de IV 2 tenemos para la densidad del mismo 

m 0.211V x M 02 + 0.79./V x Mn 2 (0.21Mq 2 + 0.79AÍjv 2 ) N . 

P = V = V = V 

que podemos poner como 


= ( 0 . 2 ^+ 0 . 7911 ^)— 


(7.5) 


• Teniendo en cuenta que en el aire húmedo la proporción relativa de nitrógeno y 
oxígeno se mantiene en la relación 21 a 79 entre si, y que Xq 2 + x^ = 1 — x^ Q 

tendremos que Xq 2 = 0.21 ^1 — x^^j ; = 0.79 ^1 — x^¿ 20 \ con lo que la 

densidad del aire húmedo está dada por 

P (h) = ((0.21 M 02 + 0.79 M N2 ) (l - x^lo) + x%l 0 M H2 o) ^ (7.6) 
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• Por tanto, del cociente entre ambas densidades 

p W _ ((0.21M O2 +0.79M W2 ) (l - zg 0 ) +x { ¡¿ 0 M H , 0 ) 




(0.21Mq 2 + 0.79Mjv 2 ) 


(7.7) 


simplificando y sustituyendo las masas moleculares M H2 q = 18 g mol 1 , Mq 2 = 
32 gmol~ x , M N2 = 28gmol~ 1 y teniendo en cuenta que x^] 0 < 1 llegamos a 




( s ) ~ x hIo + 62 (! ^ííio) < 1 


(7.8) 


por lo que llegamos a la conclusión de que el aire húmedo es menos denso que el 
aire seco. 

• Para calcular la densidad del aire húmedo haremos uso de la expresión (7.6) para 
lo que tendremos que calcular la fracción molar del vapor de agua. 


AT PfloOV 

(h) _ ^H 2 Q _ RT _ fjho 

x h 2 0 ~ at ~ py_ ~ p 

RT 


(7.9) 


donde Ph 2 o denota la presión parcial del vapor de agua y P la presión total del 


aire. 


• Ahora bien, por la definición de humedad relativa H,■ = 100 p 3 ° donde P s es la 
presión de saturación del agua a la temperatura correspondiente, por lo que la 
expresión (7.9) podemos ponerla como 


( h ) H r P s 

Xh2 ° ~ 100 x P 


y en el problema que nos ocupa 
Ah) 


_ _ 90 x n.hmmHg _ 2 

100 x 760 mmHg- Z -"‘ X iU 


(7.10) 


(7.11) 


• Sustituyendo numéricamente en la expresión (7.6) 

Alo-, = 32gmol ~ 1 x 10 ~ a kgg~ 1 = 32 x 10~ 3 kgmoP 1 
M N2 = 28 gmoP 1 x 10 ~ 3 kg g~ 3 = 28 x 10 ~ 3 kg mol~ 3 
(0.21 Mo 2 +0.79M N2 ) = 28.84 x 10 ^kgmor 1 

P laimXl01.3xl0 3 Pa atrn ~ 1 _ 11 r mnl 

RT ~ 8.314J mol -1 K _1 x (273.15+20) iC — ovo mui tn 

pW = 1.1894 kgm~ 3 

• Mientras que la densidad del aire seco (7.5) sería 

p( s ) = 28.84 x 10 ~ 3 kg mol" 1 x 41.563 molmT 3 = 1.20 kgmT 3 


(7.12) 


(7.13) 
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7-2) En un recipiente de 1 m? inicialmente vacío se inyecta 1 litro de agua 
líquida. Si el recipiente se mantiene a la temperatura constante de 
25°C', calcular el estado final de equilibrio que se alcanza sabiendo 
que la presión de vapor del agua a 25 °C es de 23,75 mmHg. ¿Y si el 
recipiente hubiera estado lleno de aire seco a la presión de 760 mmHgl 
¿Y si hubiera sido aire atmosférico con una humedad relativa del 60%?. 


*Método: Calcular la cantidad máxima de agua que pasa al estado de vapor 
aplicando la ecuación de gas ideal al vapor a la presión de saturación. 


• Consideremos en primer lugar el caso en el que el agua se inyecta en un recipiente 
vacío. 


• Para ello, calcularemos la máxima cantidad de agua que puede encontrarse en 
dicho recipiente en fase de vapor. Dicha cantidad será aquella en la que la presión 
sea igual a la presión de saturación. 


N = 


PsV 

RT 


(7.14) 


• En nuestro caso 

P s = 23,75 mmHg 
por lo que 

N = 


7 «n 3,75 n m f g -i x 101.3 x 10 3 Paatm~ 1 = 3165.6 Pa 

760 mmHg atm 1 


3165.6 Pa x lm 3 


8.314Jtoo/- 1 /!" 1 x 298.15A' 


= 1.28 mol 


(7.15) 


• Ahora bien, la masa de agua m que pasa a vapor será m = NM, donde M es la 


masa molecular y el volumen de líquido vaporizado es v 


(Z) _ rn 


por tanto 


m = 1.28 mol x 18 g mol 1 = 23. 04g (7-16) 

• es decir, sólo una pequeña proporción de 23.04 x 10~ 3 litros pasa de líquido a 
vapor. 

• Por tanto el estado final de equilibrio para la temperatura de 25°C será 

— Un volumen de U^) — o. 97 litros de agua permanece en fase líquida y el 
resto se ha vaporizado. 

— La presión es igual a la presión de saturación a esa temperatura P s = 23.75 
mmHg 
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• Si en el recipiente hubiera estado lleno de aire seco el volumen de agua evaporada 
sería el mismo ya que la ecuación (7.14) seguiría siendo válida, sólo que en este 
caso P s sería la presión parcial del vapor de agua en el recipiente. 

• La presión total (despreciando la disminución de volumen disponible para el gas 
por la presencia del agua líquida) sufriría un incremento debido a la presencia 
del vapor dado por 


P = P Q + P s = 760 mmHg + 23.75 mmHg = 783. 75 mmHg 


• Si por contrario hubiese estado lleno de aire con una humedad relativa del 60%, 
indicaría que la presión parcial del vapor de agua sería -jf = 0.6 por lo que 
existiría en fase de vapor un número de moles igual a 


N v 


0.6 P S V 
RT 


(7.17) 


es decir 0.61V del valor calculado anteriormente. Sólo pasaría a vapor 0.41V moles, 
suficientes para que la presión llegue a ser la presión de saturación. En este caso 
el estado del vapor y la presión será el mismo que el de los apartados anteriores, 
pero en agua líquida permanecerían 


V {1) = 1 - 0.4 x 23.04 x 10 -3 = 0.99 litros de agua (7.18) 



Robert Boyle (1627-1691) es considerado como el primer físico experimental o como el primer químico moderno. 
En su demostración de las propiedades elásticas del aire realizó los experimentos que llevaron a enunciar la propor¬ 
cionalidad inversa entre presión y volumen a temperatura constante, resultado conocido como ley de Boyle. 
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7-3) La ley de Boyle establece que a temperatura constante la presión de 
un gas es inversamente proporcional al volumen. Boyle realizó sus 
experimentos con aire atmosférico observando que cuando el volumen 
de una porción de gas se reducía a la mitad el aumento de presión sobre 
el mismo correspondía a la altura de la columna barométrica, (es decir 
la presión se duplicaba) ¿Qué modificación sería necesario establecer 
sobre la Ley de Boyle para el aire en función de su humedad relativa?. 


*Método: Determinar a que volumen comienza la condensación del vapor 


• La humedad relativa del aire está dada por 

H r _ Py 

100 ~ P s 


(7.19) 


donde P v es la presión parcial del vapor y P s la presión de saturación a la 
correspondiente temperatura. 

La presión total será 

P = Pv + Paire (7.20) 


donde 


N V RT 

V 


siendo N v el número de moles de vapor de agua en el recipiente. 


(7.21) 


• De (7.19), (7.21) llegamos a 


P s H r N V RT 
100 “ V 0 


(7.22) 


donde V (¡ es el volumen de aire húmedo a la presión y temperatura atmosférica. 


• Cuando disminuye el volumen la presión parcial del vapor de agua aumenta y 
llegará a la de saturación para un volumen V s dado por 


Ps = 


N V RT 

V s 


(7.23) 


• Dividiendo (7.22) y (7.23) llegamos a 


Hr_ = Vs 

100 V 0 


(7.24) 


por lo que si el volumen disminuye hasta un valor V > V s = no es necesario 
establecer ninguna corrección y el aire se comportará como una mezcla de gases 
ideales y por tanto obedece a la ley de Boyle. 
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• Cuando el volumen disminuye por debajo de V s , parte del vapor de agua se 
licuará y la presión total será 

P = P aire + P s (7.25) 

• Para el aire se seguirá verificando la ley de Boyle pero la presión parcial del 
vapor permanecerá igual a la presión de saturación por lo que la presión total 
será 

P= NaireRT +ps (7.28) 

que también podemos escribir como 

(P -P„)V = cte (7.27) 

• Numéricamente a una temperatura de 25 °C y con una humedad relativa del aire 
del 70% la presión de saturación es de 23.75 mmHg. Un volumen determinado 
de aire atmosférico podrá ser disminuido a temperatura constante hasta un valor 
V s = 0.7V o sin que la Ley de Boyle se vea alterada. Por debajo de este valor la 
expresión de la misma deberíamos escribirla estableciendo una corrección de la 
forma 

(P — 23.5 mmHg)V = cte (7.28) 


7-4) En una vivienda de 100 m 2 de superficie y 2.50 m de altura se activa 
un deshumidificador. Si la temperatura ambiente es T = 15°C ¿Qué 
volumen de agua líquida es extraído cuando la humedad relativa pasa 
del 100% al 75%?.¿Cuál es la máxima cantidad de agua que puede 
extraerse?. Presión de saturación del agua a ÍS 0 ^, P S (15°C ) = 1.704fcPa 


*MétocLo: Relacionar el número de moles de agua extraída con el cambio 
de la fracción molar y relacionar éste con la humedad relativa usando la 
definición de ésta y la definición de presión parcial. 

• La humedad relativa H r está dada por 

Hr = X 100 (7.29) 


donde 

P v es la presión parcial del vapor de agua en la masa de aire. 

P S (T ) es la presión de equilibrio líquido vapor para el agua pura a esa 
temperatura (Presión de saturación). 
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• Para disminuir la humedad relativa del 100% al 75% es necesario disminuir la 
presión parcial del vapor de agua de P s a 0.75P S . En el estado inicial (humedad 
relativa 100 %) la fracción molar de vapor de agua es 

H 

II 

(7.30) 

y en el estado final será 

(2) _ 0.75 P s 

p 

(7.31) 

• Por tanto de ambas expresiones deducimos que 

x {2) = 0.75xi 1] 

(7.32) 


Ahora bien, como la fracción molar corresponde al cociente entre el número de 
moles de vapor de agua y el número de moles total esta expresión conduce a 

= 0.75 N& (7.33) 


por lo que el número de moles de agua extraída de la habitación está dado por 

N¡Ü - = 0.25 NW (7.34) 


• El problema queda reducido pues a determinar el número de moles de agua 
presente en forma de vapor en la habitación en el estado inicial. Esto puede 
llevarse a cabo teniendo en cuenta que en dicho estado la presión parcial del 
vapor de agua es igual a la presión de saturación y por tanto 


nP 


PsV 

RT 


sustituyendo numéricamente tendremos 

^ r(l} 1.704 x 10 3 Nm~ 2 x 250 m 3 


8.314 Jmol- 1 K- 1 x (273.15 + 15) K 


177.82 mol 


(7.35) 


(7.36) 


por lo que el número de moles de agua extraída al baiar la humedad relativa al 
75% será 

N( 1} - N( 2) = 0.25 x 177.82 mol = 44.455 mol (7.37) 

que corresponden a una masa de 

m = 44.455 mol x 18 g mor 1 = 800.19 g (7.38) 

y a un volumen de agua líquida 

Vi = 800.19 3 x lO^lg- 1 =0.81 (7.39) 


• El máximo volumen de agua que puede extraerse de la habitación está dado pues 
por 

E max = 177.82 mol x 18 g mor 1 x 10 g- 1 =3.21 (7.40) 
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7-5) Encontramos en una tabla que las solubilidades del oxígeno y del ni¬ 
trógeno en el agua a 20 °C y 1 atm están dadas por 1.38xl0- 3 y 6.90 xlO- 4 
moles/litro respectivamente. Calcular cuál es el valor de la constan¬ 
te de Henry para dichas sustancias. Considerando que el aire está 
compuesto por un 21% de oxígeno y un 79% de nitrógeno calcular la 
proporción de nitrógeno y oxígeno en el aire disuelto en el agua. 


*Método: Transformar la solubilidad en fracción molar y aplicar la ley de Henry. 


• La ley de Henry establece que la presión parcial de un gas es proporcional a la 
concentración de éste en el líquido en el que se disuelve. En nuestro caso si la 
presión es 1 atm y la temperatura 20°C se disolverán en el agua 1.38 x 10 -3 
moles de oxígeno por cada litro de agua. Pero 1 litro de agua en número de 
moles corresponde a 


N = 1 1 x lOOOg l 1 x 18 g ra 0 ¿-i = 55.6 mol 


por lo que la fracción molar de oxígeno cuando la presión de éste es de 1 atm 
será 


x o 2 


1.38 x 1(T 3 
5576 


= 2.482 x 1(T 5 


(7.41) 


• Por idéntico razonamiento la fracción molar de nitrógeno será 


x N2 


6.90 x 10- 4 
5576 


= 1.241x 1Q- 5 


(7.42) 


• Las constantes de Henry estarán dadas pues por 
i atm -4x10 A atm 


K o> = *t 


K 


n 2 


P No 
X N <2 


2. 482xl0- 5 
1 atm 

1.241 XlO' 5 


= 8 x 10 4 atm 


• Cuando el agua se encuentra en contacto con el aire las presiones parciales del 
oxígeno y del nitrógeno están en una relación 21 a 79 por lo que si P es la presión 
total las respectivas concentraciones son 


x 0 2 


0.21 P 
Ko 2 


Xn 2 


0.79 P 
k n 2 


(7.43) 


por lo que la relación entre ambas corresponde a 

xo 2 0.21 Kn 2 0.21 x 8 x 10 4 aím 
xn 2 0.79 Kq 2 0.79 x 4 x 10 4 aím 


(7.44) 


• Es decir que si extraemos aire del agua su concentración relativa de oxígeno y de 
nitrógeno es de aproximadamente un 33% de oxígeno y un 66% de nitrógeno, es 
decir que el aire disuelto en agua es más rico en oxígeno que el aire atmosférico. 
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7-6) El etileno-glicol (l-2etanodiol) OH—CH 2 —CH 2 —OH es un producto uti¬ 
lizado como anticongelante en los radiadores de los automóviles. Sus 
propiedades anticongelantes provienen del comportamiento no ideal 
de su disolución en agua, sin embargo, si su comportamiento fuera 
ideal, ¿Cual sería el descenso del punto de congelación por cada gra¬ 
mo disuelto en un litro de agua?. (Calor latente de fusión del agua 
A h°f = 80 cal/g) 


*Método: Evaluar la fracción molar y calcular usando la expresión del 
descenso crioscópico. 

• Para calcular la fracción molar de la sustancia disuelta tenemos en cuenta que 
un litro de agua corresponde a 


Nl = m ° 9 
18 g mo Í _1 


55. 6 mol 


(7.45) 


Por otra parte, la masa molecular del etileno glicol M eg = 62.0 gmol -1 por lo 
que 1 gramo del mismo corresponde a mol , y por tanto la fracción molar de 1 
gramo de etilenglicol disuelto en 1 litro de agua está dada por 


z 2 = -mol 


(55.6 + 



mol = 2. 9 x 10 ~ 4 


(7.46) 


• El descenso del punto de congelación podemos calcularlo como: 


A T f = 



A h ° f 


x s 


(7.47) 


• E 11 nuestro caso 


R = 8.314 J mol 1 K 1 
T° f = 273.15 A' 

A/ij = 80 calg ” 1 = 80 calg _1 x4.18 Jcal~ l x 18 gmol -1 = 6019.2 Jmol _1 
por lo que 


8.314 Jmol^K- 1 x (273.15AT) 2 , , 

A T f = - „ mQ _ f - —2- 9 x 10 - 4 = 2 .99 x 10 ~ 2 K 

6019.2 J mol- 1 
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7-7) ¿Cuál es el peso molecular de un anticongelante cuando es necesario 
añadir 1,167 kg del mismo por cada 5 kg de agua para disminuir el punto 
de congelación del agua a —7 °C1 


*Método: Aplicar la expresión del descenso crioscópico. 


En efecto, la expresión del descenso crioscópico está dada por 


A T f = 


R(Tf ) 2 


como 


N s m s /M s __ m s M\ 


N s -\-Ni m\/M\ miM s 


sustituyendo 


A T, 


/ 


R{Tf) m s M 1 


(7.48) 

(7.49) 

(7.50) 


Ahj m\M s 

• De esta expresión podemos despejar la masa molecular del soluto M s obteniendo 

m„ í \ i 


M. - 


m i 


Ah 


f 


A T, 


f 


(7.51) 


• Numéricamente, en el problema que nos ocupa 

A hR(T° f ) 2 _ 8.314J mol-'K- 1 x (273.15A') 2 x 18 gmor 1 x ÍO^kgg- 1 


A h° 


nos ofrece 


con lo que 


80 cal g^ 1 x 18 g mol 1 x 4.18 J caí -1 


MiR(Tf) 2 , 

- 2 — = 1. 855 A kg mol 

Ah° t 


(7.52) 

(7.53) 


M s = — '9 x 1.855 K kg mol 1 x — « 62 x 10 3 kgmol 1 (7-54) 

5 kg 7 A 
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7-8) Se suele entender por constante crioscópica de un disolvente el factor 
de proporcionalidad entre el descenso crioscópico y la concentración 
de dicho disolvente. Encontramos en la bibliografía que la constante 
crioscópica del agua es 1.86 kg K mol _1 . Discutir qué expresión debe 
tomarse del descenso crioscópico para obtener dicho factor de propor¬ 
cionalidad. 


*Método: Expresar la fracción molar en términos del número de moles de 
soluto y la masa del disolvente. 


• Si dicha constante es tomada en kg K mol 1 significa que las unidades de la 
concentración son mol kg~ 1 , por lo que en términos de la fracción molar 


N s 


N„ 




N 


N„ 


Ni 

JV1 Mi 


donde Ni, mi, Mi denotan el número de moles, la masa, y la masa molecular 
del disolvente respectivamente y donde hemos aproximado en el denominador el 
número de moles totales por el número de moles del disolvente ~ wf)- 


• Por tanto la expresión del descenso crioscópico la escribimos como 

R (Tf 2 

A T f 


M 

A h° t 


Mi 


1 N. 


‘/ mi 

y el factor de proporcionalidad corresponde a 

, 2 


(7.55) 


R 


Kr = 


i T ?y 


M 


A h° f 


(7.56) 


Numéricamente para el agua 

_ 8.314J mol~ 1 K^ 1 x (273.15A') 2 x 18 x I0~ 3 kg mol~ l 
c _ 6019.2 Jmol- 1 

Ejercicio adicional: 


1.86 kgmol 1 K 


— El punto de fusión normal del agua químicamente pura es 0.00 25°é7, mien¬ 
tras que el punto de fusión del agua en equilibrio con aire atmosférico a 1 
atmósfera de presión está dado por 0 °C. Esta discrepancia es atribuida a la 
disolución de los componentes del aire atmosférico en el agua. Determinar 
los moles de aire disueltos en cada litro de agua. 
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7-9) Una disolución se denomina isotónica con otra cuando ambas tienen la 
misma presión osmótica respecto del disolvente. Un caso de particular 
importancia son las disoluciones isotónicas con el suero de la sangre 
humana. Disoluciones isotónicas de cloruro sódico (suero fisiológico) 
son utilizadas con frecuencia en distintos tratamientos médicos. Sa¬ 
biendo que el suero de la sangre humana se congela a — 0.56°C, estimar 
la cantidad de CINa que hay que añadirle a 1 litro de agua para for¬ 
mar una disolución isotónica con la sangre humana. Supóngase que 
el CINa (0=35,5, Na=23) se disocia completamente en el agua y que 
lcm 3 de suero contiene lg de agua. (Temperatura ambiente 20°C). 


*Método: Calcular la presión osmótica teniendo en cuenta que la fracción 
molar del soluto aparece como un factor que puede eliminarse en las ex¬ 
presiones del descenso crioscópico y de la presión osmótica, quedando una 
relación entre estas dos magnitudes. Considerar los iones Cl~ y Na + como 
especies químicas diferenciadas para calcular la concentración. 

• La presión osmótica está dada por 


n * s 

V o 


y el descenso crioscópico por 


R(T° f ) 2 
A T f = x s 

1 A h° 


Al dividir ambas expresiones 

n TAh° 
~ATf ~ v°{Tf) 2 


(7.57) 


(7.58) 


(7.59) 


de donde la relación entre la presión osmótica y el descenso crioscópico está dado 
por 


_ AT f Ah° 
°{T°) 2 


(7.60) 


• Con los datos para el agua tenemos a 37°C 

0.56 K x 80 cal g- 1 x 4.184 J cal" 1 x 18 gmol" 1 
lcm 3 ^” 1 x 10 -6 m 3 cm~ 3 x 18 gmoC 1 x (273.15A') 2 

que resulta 

n = 7.7918 x 10 5 Pa « 7.6 atm (7.61) 
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• Sin embargo hemos de tener en cuenta que si la disolución la preparamos a 20 °C 
la presión osmótica de la misma debe ser 


n 20 = n 37 |^ = 7. 7918 x 10 3 ¡ff = 7.3648 x 10 5 Pa (7.62) 

• En términos de la concentración, la presión osmótica está dada 

H = ^x, (7.63) 


pero en este caso al estar el CINa disociado debemos tomarlo como si fueran 
dos solutos Cl~ y Na + por lo que 


„ BT. 

n — —-{x C i- + x Na +) 

V 


Ahora bien N c ¡- = N Na + = NciNa por tanto podemos poner 


TT 

n = —~2xciNa 
V o 


es decir 

J-J _ RT 2 NciNa ~ 2RT NciNa _ 2j?T mciNa/MciNa 
V o Nh 2 o + NciNa V o Nh 2 o V o m H2 o/M H2 o 

de donde podemos despejar la masa mciNa 


mciNa 


Uv 0 m H2 oM C lNa 

2RTMh 2 o 


(7.64) 


(7.65) 


(7.66) 


(7.67) 


• Sustituyendo numéricamente, suponiendo que la disolución se prepara a T = 

20 °C (293.15/t) y tomando el dato de la presión osmótica de la expresión (7.62) 
llegamos a 

v° = (lcm 3 # -1 x 10 ~ 6 m. 3 cm ~ 3 x 18<? mob 1 ) = 1.8 x 10 ~ 5 m 3 mol~ 1 
mH 2 o = 1 kg 

MciNa = 58.5 gmol^ 1 x 10~ 3 kg g~ 1 = 58.5 x 10 ~ 3 kg mol^ 1 
Mh 2 o = 18 gmol^ 1 x 10 ~ 3 kgg~ 1 

_ 7.3648 x 10 5 Pa x 1.8 x 10 - 5 m 3 moZ" 1 x 1 kg x 58.5 x I0~ 3 kg mol' 1 
mciNa ~ 2 x 8.314 J mol-iR- 1 x 293.15 A' x 18 x 10 ^kgmol- 1 

de donde obtenemos un valor numérico 

mciNa = 8.8 x I0~ 3 kg = 8.85 (7.68) 


110 









7-10) La presión osmótica del agua del mar es II = 25, 9 atm. Comparar la 
energía necesaria para desalinizar 1 litro de agua utilizando osmosis 
inversa con la correspondiente utilizando destilación (evaporándola 
y condensándola posteriormente). Suponer la destilación i) a 20°C, 
A h° = 586 cal g~ 1 y ii) a 100°C, A h° = 540 calg~ 1 . 


*Método: Calcular el trabajo necesario para hacer pasar un litro de agua 
a través de la membrana semipermeable y comparar con el calor necesario 
para provocar la ebidlición. 

• El trabajo necesario para hacer pasar un volumen V de agua a través de una 
membrana semipermeable podemos calcularlo suponiendo que dicha agua está 
encerrada en un recipiente cilindrico y que es empujada por un pistón. La fuerza 
ejercida por dicho pistón deberá ser suficiente para vencer la presión osmótica y 
por lo tanto deberá ser igual a 

F = TIA (7.69) 

donde A es la superficie del pistón. 

• Cuando el pistón se desplaza un elemento dx un volumen de agua dV = A x dx 
atraviesa la membrana semipermeable. El trabajo realizado por el pistón está 
dado por 

SW = Fdx = HAdx = I JdV (7.70) 

y por tanto cuando un volumen V atraviesa la membrana el trabajo que es 
necesario realizar sobre el mismo es 

W = UV (7.71) 

• En el caso que nos ocupa V =11 

W = 25.9 atm x 101.3 x 10 3 Paatm~ 1 x 11 x 10~ 3 m 3 Z _1 (7.72) 

lo que nos ofrece un resultado de 

W = 2.6 x 10 3 J (7.73) 

por cada litro de agua producido. 

• Para vaporizar 1 1 de agua a 20°C' mediante aporte de calor es necesario sumi¬ 
nistrarle una cantidad de energía 


Q = mL v (7.74) 

que en nuestro caso corresponde a 

Q = lkgx 58 Qcalg~ l x 4.18 Jcar 1 = 2.4 x 10 6 J (7.75) 
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por lo que vemos que existe una relación 

% « 1000 (7.76) 

W 

Es decir, el procedimiento de vaporización resulta un factor 1000 veces energéticamente 
más costoso que el de osmosis inversa. 

• Si además tenemos en cuenta la energía necesaria para pasar el agua de 20 °C a 
100°C tendremos 

Q = mc{Tf - Ti) + mL v (7.77) 

Que en nuestro caso 

c = 1 cal g^K^ 1 x 4.18 J caí -1 x 10 3 g kg~ x = 4180 J kg~ l K~ x 
L v = 540 cal g- 1 x 4.18 J cal~ l x 10 3 # kg~ x = 2.2572 x 10 6 J kg- 1 

nos ofrece 

Q = I%x4180 Jkg^K^x 80K + lkg x 2. 2572 xlO 6 J kg- 1 = 2.6x10 6 Jkg- 1 
es decir, el proceso requiere aún un mayor suministro de calor. 



Jacobus van’t Hoff (1852-1911), químico holandés que fue pionero en la aplicación de los métodos de la 
Termodinámica de Gibbs a la química. Estudió el fenómeno de la osmosis llegando a la conclusión de que las diso¬ 
luciones diluidas tienen propiedades que se asemejan formalmente a las de los gases. 
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Capítulo 8 


Reacciones Químicas. 


8.1 Objetivos 

• Comprender las relaciones entre las magnitudes termodinámicas que describen 
el equilibrio químico: grado de avance, afinidad propia, constante de equilibrio, 
calor de reacción. 

• Comprender, si desde el punto de vista termodinámico, una reacción química es 
posible o no y cual es su grado de realización en el equilibrio. 

• Analizar la energía que es posible obtener de una reacción química. 

• Analizar como influye la temperatura y la presión en el equilibrio químico. 

8.2 Problemas 


8-1) Una reacción de interés en los tubos de escape de los automóviles, así 
como en la producción industrial de hidrógeno es la denominada del 
desplazamiento del agua 

H 2 0(g) + CO(g) ^ C0 2 (g) + H 2 (g) (8.1) 

Teniendo en cuenta los datos de la tabla, analizar cuáles son las condi¬ 
ciones de equilibrio para esta reacción en condiciones estándar y cómo 
influirá una variación de la presión en las condiciones de equilibrio. 



H 2 0(g) 

CO(g ) 

C0 2 (g) 

H 2 (g) 

AH° f (kJ mol 1 ) 

-241.818 

-110.525 

-393.509 

0 

s'ftJmol-'K- 1 ) 

188.8 

197.7 

213.7 

130.7 
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*Método: Calcular la afinidad propia y la constante de equilibrio en términos 
de las entalpias y entropías de reacción. Evaluar la derivada de la afinidad 
propia con respecto a la presión. 

• Dado que toda reacción química puede ponerse como combinación lineal de las 
reacciones de formación de cada uno de sus compuestos siendo los coeficientes 
estequiométricos los coeficientes de dicha combinación, la entalpia de reacción 
puede ponerse como combinación lineal de las entalpias de formación de cada 
uno de los compuestos: 

A H° = A H° C02 + AH° H2 - A H° co - A H°^ 0 (8.2) 

A H° = (-393.509 - (-110.525 ) - (-241.818 ))kJmoG l = —41.166 kJ mol" 1 

• La entropía de reacción está dada por: 

A 0° , o o o 

Ao — s C q 2 + s H2 - s H2 q - s co — 

= (213.74 + 130.684 - 188.825 - 197.674) Jmol^R- 1 = -42.075 

• En el cálculo de la energía libre propia 

A G° = (A H° - T°AS°) (8.3) 

el término A H° = —41.166 x 10 ¿ Jmol~ 1 es negativo, mientras que el término 
—T°AS° = -298.15 A' x (—42.075 J mol~ 1 K^ 1 ) = 12.545 x 10 3 Jmol~ x repre¬ 
senta una contribución positiva, sin embargo en condiciones estándar el primer 
término es dominante dando lugar a 

A G° = -28621 J mol~ l (8.4) 

un valor negativo grande de A G° por lo que podemos intuir que la reacción 
estará muy desplazada hacia la derecha. 

• Esto podemos corroborarlo calculando la constante de equilibrio 

K = e~ AG ' RT = e 8.3i4 jmo¡-ijf-i298.i5x = 1.04 x 10° « oo (8.5) 


por tanto la reacción es prácticamente completa. 




Para analizar las variación con la presión de las condiciones de equilibrio, tenemos 
en cuenta que 


, dAG\ 

( dP ’ T ( dP 




v ó v ó 


( 8 . 6 ) 


como la reacción tiene lugar entre gases ideales Vj = igual para todos los 
gases, por lo que teniendo en cuenta que ^ v j = 0 


dAG 
8P ’ T 



(8.7) 


es decir en este caso un cambio de presión no produce variación de las condiciones 
de equilibrio. 
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8-2) Considere la reacción 


2Hg(l) + 0 2 (g)^Hg0(s) (8.8) 

Teniendo en cuenta los datos aportados en la siguiente tabla, analizar 
desde el punto de vista termodinámico las condiciones de equilibrio 
de la reacción anterior. 



AH°(kJ/mol) 

s°( J ) 

C V ( rn.nl. k) 

HgO 

-90.83 

70.29 

44.06 

o 2 

0 

205.138 

29.35 

Hg 

0 

76.02 

27.98 


*Método: Calcular la energía libre propia en términos de la entalpia y la 
entropía y suponer que éstas no varían con la temperatura para estimar a 
que temperatura la energía libre se anula. 

• La reacción propuesta reacción de formación del HgO por lo que su entalpia 
propia corresponde a la dada en la tabla 

A H° = -90.83 x 10 3 J mor 1 (8.9) 

La entropía de reacción etá dada por 

A S° = s° Hg0 -\s° 0a -s° Hg 

que sustituyendo los valores de la tabla nos ofrece 

AS* * 0 = -108.3 Jmol^K- 1 

E1 término TAS° en condiciones estándar (T = 298.15A') resulta 
TAS° = -32.29 x 10 3 J mol" 1 

• Por tanto, la energía libre A G° = AH° — TAS° está dada por 

A G° = -90.83 x 10 3 J mor 1 - (-32.29 x 10 3 JmoL 1 ) = -58.54 Jmo/’ 1 

• Vemos que es una reacción exotérmica A H° < 0 con una entropía de reacción 
negativa AS" < 0 por lo que el término — TAS° tiene una contribución posi¬ 
tiva a A G°. Sin embargo dados los valores absolutos en condiciones estándar 
el término A H° < 0 domina el resultado por lo que A G° < 0,y como se trata 
de una reacción entre fases puras podemos decir que desde el punto de vista 
termodinámico es una reacción espontánea. 


( 8 . 10 ) 

( 8 . 11 ) 

( 8 . 12 ) 


115 













• Bajo la hipótesis de que A H° y A S° varían poco con la temperatura la reacción 
se invertiría a la temperatura en la que A G° = A H° — TAS° = 0, es decir 


• El término Ac° 


Ac° = ( 44.06 


^ A H° -90.83 x 10 3 J mol- 1 coo ^ 
T ~ ~AS° ~ -Í08.3 Jmol- 1 “ 838 ‘ 7A 

( c °p) Hg o - h ( c °p)o, - 1 ( c p) ffo está dado P° r 
1 


(8.13) 


P/Oa 2 \ pj Hg 

x 29.35 - 27.98 ) J moH 1 K~ 1 


lAJmol-'K- 1 


Es una magnitud positiva pero pequeña, por lo que el calor de reacción ) p = 

A c° aumenta con la temperatura (disminuye en valor absoluto) pero en pequeña 
magnitud, igual razonamiento podemos llevar a cabo con la variación de AS°. 


• En relación a la dependencia con la presión 


<9AG°\ 

dP ) T 




(8.14) 


en este caso HgO y Hg son fases condensadas (sólido y líquido) y su volumen 
molar es pequeño, por lo que el factor dominante corre a cargo del volumen 
molar del O 2 . 


dAG°\ 
~dP ~)t 




1 RT 

2 ~P~ 


(8.15) 


• Es decir, la energía libre disminuye al aumentar la presión, por lo que un aumento 
de presión favorece la realización de la reacción. 


El óxido de mercurio HgO (rojo de mercurio) fue una sustancia clave pa¬ 
ra el descubrimiento del oxígeno, llevado a cabo por Priestley en 1774 e 
interpretado en los términos actuales por Lavoisier unos años más tar¬ 
de. Priestley calentó mercurio metálico en presencia de aire formándose el 
óxido y recogió el gas desprendido por éste cuando se calentaba a una tem¬ 
peratura superior volviendo a restituirse el mercurio. Desde el punto de 
vista de nuestro análisis interpretamos que a temperatura ambiente, si bien 
la reacción es permitida por las leyes de la Termodinámica, la cinética de la 
misma es muy lenta y por tanto ésta no se realiza en términos apreciables. 
Cuando se aumenta la temperatura calentándolo al fuego la cinética se 
acelera y por tanto se comienza a formar el óxido en la superficie. Cuando 
se aumenta la temperatura aún más (Priestley utilizó una lupa para con¬ 
centrar los rayos solares) el equilibrio se invierte y el óxido de mercurio se 
descompone en mercurio y oxígeno, un gas que Priestley interpretó como 
aire flogisticado y al que posteriormente Lavoisier identificaría no como un 
tipo de aire sino como un componente más de la atmósfera, al que bau¬ 
tizaría como oxígeno (generador de ácidos) en la creencia (que luego se 
mostraría errónea) de que era el principio generador de todos los ácidos. 
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8-3) El pentóxido de dinitrógeno N 2 O5 es un sólido que puede obtenerse por 
deshidratación del ácido nítrico HNO 3 y se descompone de la forma 

2N 2 0 5 (s) G 4^02(5) + O2(g) (8.16) 

Determinar la entalpia propia de reacción. ¿Verifica esta reacción el 
principio de Thompsen y Berthelot? 



N 2 0 5 (s ) 

no 2 

0 2 

A G° f (kJmoG 1 ) 

113.9 

51.3 

0 

s°(J mol~ 1 K 1 ) 

178.20 

240.06 

205.14 


*Método: Aplicar las relaciones entre energía libre, entalpia y entropía. 

• La energía libre propia está dada por 

AG° = A G°q 2 + 4A G° N02 2AG° N2 o 5 (8.17) 

que numéricamente nos ofrece 

A G° = (4 x 51.3-2 x 113.9) kJ moG 1 = -22.6 kJmor 1 (8.18) 

lo que nos indica que A G° < 0 y por tanto la reacción será espontánea 

• La entropía propia de reacción A S° está dada por 

A S° = Sq 2 + 4 s° N q 2 - 2s^ 2 o 5 (8.19) 

que numéricamente resulta 

A S° = (205.14 + 4 x 240.06- 2 x 178.20 )JmoG 1 K~ t = 808.98 J moG 1 Jé' 1 

• En condiciones estándar 

TAS° = 298.5A' x 808.98 Jmol^R- 1 = 2.4148 x 10 5 J moG 1 (8.20) 

• La entalpia propia de reacción A H° está dada por AH 0 = A G° + TAS° que 
numéricamente nos conduce a 

A H° = —22. 6 k J moG 1 + 2.41 x 10 1 kJ moG 1 = 218.4 kJmoG 1 (8.21) 

• Es una reacción endotérmica (A H° > 0) , siendo un ejemplo de una reacción 
espontánea (A G° < 0) que no es exotérmica y por tanto no verifica el principio 
de Thompsen y Berthelot. Esto ocurre porque en este caso el término TAS° 
domina frente al A H° es decir TAS° > AH°, por lo que A G° tiene diferente 
signo que A H°. 
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( 8 . 22 ) 


8-4) Considérese la reacción 

\ n 2 (g) + \o 2 {g) ^ NO{g) 

a T = 298 A' y 1 atm, AH° = 90250 J mol" 1 y A G° = 86550 J mol" 1 , calcular 
la constante de equilibrio y su variación con la temperatura supuesto 
gas ideal diatómico (c p = ^R) para cada componente.Evaluar el grado 

de avance de la reacción en condiciones estándar ya T = 2000A". 


*Método: Calcular la constante de equilibrio a partir de la energía libre 
propia. Calcular la variación de la constante de equilibrio en función de la 
expresión integrada de la entalpia de reacción. 


En efecto: 


K(T, P) = exp{-AG°/RT} 


(8.23) 


de donde 


A '< T ’ P > = ^ 8 .314;T4 m °'x 2 9 8 A- i = 6 - 74 X 1(r “ < 8 - 21 > 


• En términos del grado de avance 

n n 2 No 2 - \ _ N no 

_I “ _i “ 1 

2 2 

Nn 2 =No 2 = \{ 1 - 0 
Nno = f, 

por lo que la constante de equilibrio es 

xno 
1/2 1 

Vn 2 x o 2 V2V- w; • 12 
sustituyendo el valor numérico 


= £ 


K(T P) = ~^-=- - - 

’ x\l 2 xT (éa-oméa-O) 1 / 2 1(1-0 


2 ? 


6. 74 x 10~ 16 


cuya solución es 


i-É 

£ = 3.37 x 10~ 16 


(8.25) 

(8.26) 

(8.27) 

(8.28) 

(8.29) 

(8.30) 


es decir a la temperatura y presión estándar la reacción es prácticamente impo¬ 
sible. 


118 












• Para analizar la variación con la temperatura tenemos en cuenta que 

rT A H°(T,P) 


ln 


K(T,P) — ln K (T 0 , P) = f 

Jt 0 


RT 2 


-dT 


Ahora bien 


A H°(T, P) = A H°(T a , P)+ £ v 3 {c° P ) 3 dT 


(8.31) 


(8.32) 


X v Á c °p)3 - c °p no ~ ñ c °Po 2 ~ 9 C Pn 


(8.33) 


En este caso 

3 

que en la hipótesis de gas ideal diatómico para cada componente se traduce en 

X^( c p),=° (8-34) 

3 

es decir el calor de reacción es prácticamente constante, con lo que 

rT 


ln K(T, P) = ln K(T 0 , P) + A H°(T 0 , P) J ^dT 

/ 

= lnK(T 0 ,P) + AH°(T 0 ,P)- 


>T 0 

T — T a 
RTT n 


(8.35) 


A T = 2000 K 


1„ K(T, P) = -34.933 + ^0^-■ 3(20005 - 29844)^ = _ 3 ^ 93 


8.314 J mol~ 1 K~ 1 x 2000 A' x 298 A' 


de donde 


K(T, P) = exp{—3. 93} = 1.96 x 10 " 2 
y en términos del grado de avance 

2 £ 


1-e 


1.96 x 10 


-2 


que nos da 


Í = 0.9 6 x 10 ~ 2 
una concentración aún pequeña pero que es apreciable. 


(8.36) 

(8.37) 

(8.38) 


Los óxidos de nitrógeno forman una parte importante de la contaminación 
en las grandes ciudades y zonas industriales. Si bien el nitrógeno y oxígeno 
no se combinan a temperatura ambiente, su combinación es significativa a 
las temperaturas que se alcanzan en el interior del motor de los automóviles 
y, si bien cuando los productos de la combustión escapan a la atmósfera 
deberían descomponerse nuevamente en nitrógeno y oxígeno, ésta descom¬ 
posición es lenta por lo que los tubos de escape deben estar provistos de 
catalizadores que aceleran la descomposión de los mismos. 
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8-5) Estimar la temperatura máxima alcanzada en la combustión del ace¬ 
tileno C 0 H 2 con la cantidad estequiométrica de oxígeno puro teniendo 
en cuenta que el calor de combustión del acetileno C 2 H 0 a 25 °C es 
A H° = —300 kcalmoP 1 y las capacidades caloríficas del C02,02, H 2 O y 
N 2 pueden ponerse de la forma ^ = A + BT con valores de A y B dados 
en la siguiente tabla. 



C0 2 

0 2 

H 2 0 

n 2 

A 

5.457 

3.639 

3.470 

3.280 

10 ¿ B 

1.045 

0.506 

1.450 

0.593 


(8.39) 


* Método: Considerar la variación de entalpia en el proceso de reacción a 
temperatura constante y luego la variación de entalpia experimentada por 
los productos de la reacción hasta una máxima temperatura. Tener en 
cuenta que en el proceso experimental la variación global de entalpia será 
nula. La presencia de nitrógeno no influye en los valores de la reacción, 
aunque si en el aumento de temperatura de los productos. 


• La reacción 


C 2 H 2 + -0 2 H 2 0 + 2 C0 2 


(8.40) 


tiene un calor de reacción A H° = —300 kcalmoP 1 a T = 25°C lo que significa 
que si la reacción tiene lugar en dichas condiciones, se libera una cantidad de 
300 kcal por cada mol de C 2 H 2 consumido. 

• Sin embargo en la práctica la reacción no tiene lugar a temperatura constante 
sino que las condiciones se asemejan más a un proceso adiabático. 

• Nosotros podemos imaginar un proceso en dos etapas que contemple: 

— a) B\ —> B 2 : Reacción a temperatura constante T = 25°C 

— b) B 2 —► B 3 : Aumento de temperatura de los productos hasta la tempera¬ 
tura final. 


De manera que el proceso B\ —> B$ reproduzca el proceso experimental de un 
proceso adiabático a P cte que en la hipótesis de que es reversible tenemos: 


H, = ff, 


(8.41) 


• Así pues la variación de entalpia en el proceso —► B 2 (reacción a T = cte) 

H 2 — H\ = A H° = —300 kcal mol ^ 1 (8.42) 


120 











• Variación de entalpia en el proceso B 2 —> B 3 (Aumento de temperatura de un 
número de moles de los productos igual a los coeficientes estequiométricos). 


H 3 - H 2 = / 3 V v jCp dT 

Jt 2 j- r+ i 

(8.43) 

que está dado por 


H 3 — H 2 = í (2(cp)co2 + ( cp)h20 ) dT 

Jt 2 

(8.44) 

que teniendo en cuenta que 


(c P )co 2 = i?(5.457+ 1.045 x 10" 3 T) 

(8.45) 

(c p )h 20 = R( 3-470 + 1.450 x 10^ 3 T) 

(8.46) 


nos lleva a 


H 3 -H 2 = í 3 (2i?(5.457 + 1.045 x 10~ 3 T) + i?(3.470+ 

J 298 

+1.450 x 10" 3 r)) cLT 

H 3 ~H 2 = —8829. 5 + 28. 581T 3 + 3. 517 x 10^ 3 T 3 2 (8.47) 

de donde haciendo 

H 3 -H 2 = - (H 2 - Hí) (8.48) 

se obtiene 

-8829. 5 + 28. 581T 3 + 3. 517 x 10~ 3 T| = 300 x 10 3 caZ mor 1 (8.49) 
cuya solución es 

T 3 = 6150. 5 K (8.50) 

• Sin embargo en la práctica esta temperatura no es tan alta debido a que habi¬ 
tualmente la combustión del acetileno es alimentada con aire y no con oxígeno 
puro. En este caso suponiendo que añadimos la cantidad exacta de aire para 
provocar una combustión de todo el oxígeno tenemos 

C 2 H 2 + — O2 + 9.4 N 2 H 2 0 + 2C0 2 + 9.4./V2 (8.51) 

donde hemos tenido en cuenta que la relación nitrógeno oxígeno en el aire es 
y por tanto por | moles de oxígeno debe haber || x | = 9.4 moles de nitrógeno. 

• En el proceso B\ —► B 2 de reacción a temperatura constante el nitrógeno no 
tiene ningún efecto por lo que H 2 — H 3 = —300 kcal. 
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• Sin embargo en el calentamiento de los productos de la reacción sí que es sig¬ 
nificativa la presencia de nitrógeno cuya capacidad calorífica molar escribimos 
como 

{c p )n2 = R{ 3.280 + 0.593 x 1(T 3 T) (8.52) 

lo que nos lleva a 

H 3 -H 2 m í 3 (2 R (5.457 + 1.045 x 1(T 3 T) + R (3.47 + 1.45 x 1(T 3 T) 


J 298 

+ 9.4 (R (3.28 + 0.593 x 1(T 3 T)) dT (8.53) 

H :i -H 2 = -27578. + 89.844T 3 + 9.055 x KT 3 T 3 2 (8.54) 

de donde haciendo 

H 3 -H 2 = - (H 2 - H 3 ) (8.55) 

se obtiene 

-27578. + 89.844T 3 + 9. 055 x 10 _3 2f = 300 x 10 3 (8.56) 

T 3 = 2835. 7 K (8.57) 


• Incluso en el procedimiento habitual la cantidad de oxígeno no es exacta y por 
tanto también habría que añadir dicha corrección, así como que el proceso no es 
completamente adiabático, por lo que en la práctica la temperatura es significa¬ 
tivamente menor. 

• Ejercicio adicional. 

— Obtener la máxima temperatura en la combustión del butano (C 4 H 10 ) sa¬ 
biendo que la entalpia propia de su reacción de combustión es A H c = 
-2881.9 kJ mol- 1 
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8 -6) Analizar las condiciones de equilibrio de la reacción de descomposición 
de la piedra caliza ( carbonato cálcico [ CaCO 3 ]) para producir cal viva 
(óxido de calcio [CaO(s)]) 

CaC0 3 {s) <=t CaO{s) + C0 2 {g) (8.58) 

teniendo en cuenta que la presión parcial del dióxido de carbono en 
la atmósfera es Pco 2 = 32 Pa. 



A H^kJmoR 1 ) 

s°{JmoR 1 R- 1 ) 

c^JmoRtR- 1 ) 

CaC0 3 (s) 

-1206.92 

92.9 

81.88 

CaO(s) 

-635.09 

39.75 

42.80 

C0 2 

-393.51 

213.74 

37.11 


*Método: Relacionar la energía libre con las entalpias y entropías de re¬ 
acción en condiciones estándar. Analizar la variación con la temperatura 
a presión estándar. Para analizar la variación de la energía libre con la 
presión tener en cuenta que para la sustancia en fase gaseosa el potencial 
químico depende logarítmicamente de la presión 

• Cualquier reacción puede ponerse como combinación lineal de las reacciones de 
formación de los compuestos respectivos con lo que siguiendo la ley de Hess la 
entalpia propia de esta reacción está dada por 

A H° = A H° Ca0 + A H° C02 - A H° CaC03 (8.59) 

numéricamente 

A H° = (-635.09 + (-393.509) - (—1206.92)) kJ mol" 1 = 178.32 kJ mor 1 

(8.60) 

donde hemos de matizar que las unidades mol están referidas a la descompo¬ 
sición de 1 mol de carbonato cálcico. 

• La entropía propia de reacción, dada por AS° = JA Ujsf donde v : ¡ son los coefi¬ 
cientes estequiométricos y s° las respectivas entropías molares podemos ponerla 
como 

A S° = (39.75 + 213.74 - 92.9) J mol^R- 1 = 160. 59 J mol^R- 1 (8.61) 

• La energía libre o función de Gibbs propia de reacción (AG° = JA u.j /i" ) está 
dada por 

A G° = A H° - T°AS° (8.62) 

con lo que 

A G° = 178.32xl0 3 J mo/" 1 -298.15Ax 160. 59 J moR 1 R^ 1 = 1.3x10 5 J moR 1 
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• La reacción tiene lugar entre fases puras, y la función de Gibbs propia es positiva 
A G° » 0 ( A° << 0) y por tanto la función de Gibbs es una función creciente 
con el grado de avance por lo que en condiciones estándar el equilibrio se produce 
para grado de avance cero, es decir la caliza no se descompone espontáneamente 
en óxido de calcio (cal viva) y dióxido de carbono. Por el contrario si el óxido de 
calcio se pone en contacto con el aire (y por tanto con el CO 2 ) espontáneamente 
se transforma en caliza. 


• Para analizar la variación con la temperatura de las condiciones de equilibrio 
analizaremos en primer lugar la simplificación de que la entalpia y la entropía 
de reacción no dependen de la temperatura. Así pues 

A G°(T, P°) = A H° - TAS° (8.63) 

por lo que el equilibrio se invertirá para aquella temperatura en la que 


es decir 


A G°(T) = A H° - TAS° = 0 
A H° 

T = 


AS' 0 


que tomando los valores calculados anteriormente (8.60),(8.61) 

T = 178.32 x 10 a J mor 1 = mft 4g 
160.59 J mol~ 1 K 1 


(8.64) 

(8.65) 

( 8 . 66 ) 


• Sin embargo liemos de tener en cuenta que en contacto con la atmósfera la 
reacción no tiene lugar en condiciones de presión estándar ya que la presión 
parcial del CO 2 es Pco 2 = 32 x 10” 5 6ar. Es necesario pues analizar la variación 
de la energía libre con la presión. Para ello tenemos en cuenta que el potencial 
químico del CO 2 que es la única sustancia que se encuentra en fase gaseosa, a la 
presión dada puede ponerse de la forma 

úco 2 (T°, Pco 2 ) = lP(T°. P°) + RT° ln ^ (8.67) 

donde P° se refiere a la presión estándar de P° — 1 bar, y que para los restantes 
compuestos al ser sólidos el potencial químico varía muy poco con la presión, 
con lo que la energía libre de reacción podemos escribirla como 

c P 

AG°(T°, P) = £ Vj¡jL°(T°, P°) + RT° ln — (8.68) 

3 =1 

es decir 

A G°(T°, P ) = AG°(T°, P°) + RT° ln (8.69) 

donde hemos de tener en cuenta que en este caso P se refiere a la presión parcial 
del CO 2 en la misma. Sustituyendo numéricamente A G°(T°,P°) = —1.3044 x 
10 5 Jmol- 1 ; R = 8.314 J mol^R- 1 ; T° = 298.15 K;P C o 2 = 23 x 10“ 5 5ar 

A G°(T°,Pco 2 ) = -1.09 x 10 5 Jmor 1 (8.70) 
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• La afinidad propia en este estado sigue siendo negativa, por lo tanto la reac¬ 
ción tampoco tiene lugar espontáneamente a la temperatura estándar y presión 
parcial del CO 2 en la atmósfera. 


• Nos planteamos ahora estimar a qué temperatura se invierten las condiciones de 
equilibrio cuando la presión es la presión parcial del CO 2 en la atmósfera lo cual 
tiene lugar cuando 

AG°(T, P C o 2 ) = 0 (8.71) 

• Para ello tenemos en cuenta que según el razonamiento anterior 

A G°{T,P C o 2 ) = AG° (T, P °) + RT ln (8.72) 


Si utilizamos para A G°(T,P°) la aproximación (8.63) planteada al inicio del 
problema A G°(T, P°) = A H° — TAS° podemos poner 

A G°(T, P C o 2 ) = A H° - TAS° + RT ln (8.73) 


y por tanto el equilibrio se invierte para la temperatura 


T = 


A H° 


AS° - Rln~^ 


(8.74) 


que nos ofrece el resultado 


178.32 x 10 3 Jmol- 1 

160.59 J moZ” 1 /!' -1 — 8.314 J mol~ 1 K~ 1 ln 23 x 10 -5 


774.5 K (8.75) 


que es aproximadamente la situación en un horno de la cal. 

Ejercicio adicional. 


— Realizar el cálculo de la temperatura de inversión del equilibrio incorporan¬ 
do la dependencia de la entalpia de reacción con la temperatura evaluada 
a partir de los valores de las capacidades caloríficas dadas en el problema. 


El proceso de fabricación de la cal es conocido desde tiempo inmemorial. 
La piedra caliza (carbonato cálcico [ CaCO 3]) es calentada en un horno para 
producir cal viva (oxido de calcio [ CaO ]). La cal viva tuvo diferentes usos, 
por ejemplo como desinfectante en momentos de epidemias, pero su mayor 
utilización ha sido como como material de construcción mezclándola con 
agua para producir cal muerta (liidróxido de calcio). El hidróxido de calcio 
reacciona con el dióxido de carbono del aire volviendo a dar nuevamente 
carbonato cálcico. Si el hidróxido de calcio se mezcla con arena, el sólido 
mantendrá cohesionado los granos de ésta, resultando en un material más 
resistente para su uso en la construcción. 
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Germain Henri Hess (1802-1850) nacido en Ginebra (Suiza) pero que desarrolló su prin¬ 
cipal trabajo en San Petersburgo (Rusia) sentó las bases de la Termoquímica. Estudió 
medicina en Dorpat (hoy Tartu (Estonia) si bien se trasladó a Estocolmo donde fue dis¬ 
cípulo de Berzelius, con quien profundizó en sus conocimientos de química. Realizó dis¬ 
tintos experimentos midiendo el calor liberado en las reacciones químicas estableciendo 
en 1840 la ley que lleva su nombre: 

“El calor absorbido o liberado en una combinación química es siempre el mismo bien 
sea la combinación realizada directamente o bien si tiene lugar indirectamente a través 
de diferentes etapas." 

La ley de Hess permitió la determinación indirecta del calor absorbido en un gran núme¬ 
ro de reacciones inaccesibles de forma experimental, acumulándose pronto un gran 
número de datos termoquímicos. Estos datos fueron fundamentales para el estableci¬ 
miento en la segunda mitad del siglo XIX de una interacción directa de la Termodinámica 
con la Química lo que condujo al establecimiento y aceptación de las condiciones de 
equilibrio químico y del Tercer Principio de la Termodinámica. 
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Capítulo 9 

Tercer principio. 


9.1 Objetivos 

• Comprender el enunciado del tercer principio. 

• Comprender el procedimiento de cálculo de entropías absolutas. 



Walther Nemst (1864-1941) enunció con su teorema del calor lo que posteriormente sería considerado como el 
Tercer Principio de la Termodinámica. Su interés estaba centrado en el estudio de las reacciones químicas a alta 
temperatura con el objeto de fabricar un diseño novedoso de bombilla eléctrica. Paradójicamente el teorema del 
calor derivó su interés al estudio de las bajas temperaturas al ponderar las consecuencias que su teorema tendría en 
dicho rango, en lo cual tuvo un papel importante la predicción de Einstein de que las capacidades caloríficas debe¬ 
rían anularse en las proximidades del cero absoluto. 
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9.2 Problemas 


9-1) Probar que en un sistema magnético la variación de la magnetización 
con la temperatura manteniendo el campo constante es nula en el cero 
absoluto. 


*Método: Relacionar dicha variación con la variación de entropía utilizan¬ 
do las relaciones de Maxwell. 

En efecto, la variación de la magnetización con la temperatura a intensidad 
de campo constante podemos ponerla de la forma: 


dM = 


ídM\ 

\dTjB. 


dT 


(9.1) 


ahora bien, para un sistema en el que la magnetización sea la única coor¬ 
denada de trabajo relevante, podemos definir la función de Gibbs de igual 
manera a como se define para un sistema expansivo, salvo que en este caso 
las variables naturales serán la temperatura y la inducción magnética en 
el vacío B a . 


G = G(T, B 0 ) (9.2) 

La derivada (^f) b corresponde a una de las relaciones de Maxwell ya 
que 


ídM\ í d 2 G \ 
\dT ) H ~ \dT8B 0 ) 


d 2 G 

0B o dT 



(9.3) 


Aplicando el Tercer Principio de la Termodinámica, en las proximidades 
del cero absoluto un proceso isotérmico es asimismo isoentrópico, por lo 
que 


y por tanto 


dS\ 

dB 0 ) T ~ 

(9.4) 

Cal 

U¡ 

II 

O 

(9.5) 
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9-2) El estaño es susceptible de presentarse en fase sólida en dos fases 
diferentes denominadas estaño gris y estaño blanco. A la presión 
ordinaria la temperatura de transición es T = 13.2 °C Medidas de 
las capacidades caloríficas de ambas fases ofrecen que para el estaño 
blanco ^-dT = 46.76 Jmol^R -1 mientras que para el estaño gris 

I? T¡kdT = 38.62 Jmol~ 1 K~ 1 . Determinar el calor latente de la transi¬ 
ción. 


*Método: Considerar la variación de entropía en un ciclo entre la transi¬ 
ción T = 286.4 K y T = 0 A'. 

• Denotemos por s*- 21 la entropía del estaño blanco a la temperatura de transición 
y por la correspondiente al estaño gris. El calor latente está dado por l\ 2 = 
T(s^ 2 - ) — sW). Ahora bien, ambas entropías pueden relacionarse considerando un 
ciclo en el que ambas sustancias son llevadas hasta el cero absoluto donde la 
entropía debe ser la misma para ambas fases 

S C) __ s (!) — (g( 2 ) _ gÓ 2 )) q. S C) _ s (!) _ (gW — gW) (9.6) 

donde el subíndice s 0 representa la correspondiente entropía a la temperatura 
del cero. 

• Ahora bien en virtud del tercer principio s^ — s^ = 0 y por otra parte 

(s( 2 ) - sí 2) ) = / ( f e ^dT = 46.76 Jmol^R- 1 
(s« - s { o ] ) = /,f e %dT = 38.62 Jmol ” 1 A'~ 1 
con lo que 

g(2)^s (1 ) — 46.76 Jmor 1 K- 1 -38.62 Jmol- 1 R- 1 =8.14 Jmol^R- 1 (9.7) 
por tanto el calor latente está dado por 

l V2 = 286.4 A' x 8 .14 J moR 1 R- 1 = 2331.3 J mol" 1 (9.8) 


El estaño blanco es la fase estable a mayor temperatura mientras que por debajo 
de 13°C la fase estable es la del estaño gris. Dado que el estaño gris se desme¬ 
nuza fácilmente, pierde su utilidad como elemento metálico y los instrumentos o 
soldaduras realizadas con el mismo se deterioran. Si bien el estaño blanco puede 
permanecer en estado metaestable durante un tiempo prolongado por debajo de 
la temperatura de transición, ésta se acelera por la presencia de una pequeña 
cantidad de estaño gris, que puede moverse con el polvo del ambiente. Esto hace 
que una vez que un elemento en un recinto se deteriora convirtiéndose en estaño 
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gris, éste se transmita rápidamente entre los restantes por lo que antiguamen¬ 
te cuando el fenómeno era conocido pero no su interpretación termodinámica, 
éste se interpretaba como una enfermedad denominada la "peste del estaño”, 
dándosele incluso una interpretación diabólica cuando el fenómeno ocurría en un 
reciento religioso como por ejemplo en los tubos de órgano de una catedral. 

Este fenómeno lo sufrieron las tropas de Napoleón cuando invadió Rusia en 1812 
ya que los botones de sus abrigos eran hechos de estaño por lo que se deterioraban 
fácilmente en el frío invierno ruso. A este fenómeno también se le atribuye buena 
parte del fracaso de la expedición de Robert Scott en 1912 a la Antártida quien 
falleció con varios de sus compañeros en su regreso al campamento base después 
de haber alcanzado el polo sur y comprobar que el noruego Roald Amundsen 
se les había adelantado por sólo unas semanas, siendo éste el primer hombre en 
alcanzar dicho punto. Los detalles de la expedición, conocidos al encontrarse el 
diario que Scott mantuvo escribiendo hasta poco antes de su muerte, permitieron 
conocer que en el regreso a su base se quedaron sin combustible para calentarse ya 
que éste se había derramado por el deterioro de las latas en que lo almacenaban 
al estar soldadas con estaño. 


La enfermedad del estaño puede minimizarse y la temperatura de transición 
disminuirse utilizando aleaciones de estaño con otros metales. Durante muchos 
años se ha utilizado en la soldadura como por ejemplo la realizada en dispositivos 
electrónicos una aleación de estaño y plomo. Sin embargo ésto genera en el 
medio ambiente problemas con los residuos cada vez más abundantes, por lo que 
recientemente se ha incrementado la investigación en la búsqueda de elementos 
de soldadura alternativos más respetuosos con el medio ambiente. 



Las transiciones, estructura y diversas utilidades del estaño fueron ampliamente estudiadas por Emest Julius Cohén 
(1869-1944) profesor de la Universidad de Utrech, que realizó grandes esfuerzos por reestablecer relaciones amis¬ 
tosas entre científicos de países separados por las guerras principalmente en los años posteriores a la primera güe¬ 
ña mundial. Lamentablemente su reconocido pacifismo no le evitó morir en la cámara de gas de Auschwitz en 
1944, tras haber sido encarcelado con otros muchos judíos holandeses. 
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Capítulo 10 

Conducción del calor 

10.1 Objetivos 

• Comprender la aplicación de la ley de Fourier en diversas situaciones. 

• Analizar el flujo de calor en distintos cuerpos dependiendo de su geometría. 

• Comprender la aplicación de la ley de Stefan Boltzmann. 

• Comprender la aplicabilidad conjunta de la ley del enfriamiento de Newton y la 
ley de Fourier para aplicarla a problemas en los que se den los mecanismos de 
condución del calor y radiación. 


10.2 Problemas 


10-1) Supongamos una pared que está constituida por dos materiales de 
diferente conductividad térmica k\ y k' 2 - el primero de espesor x\ y el 
segundo de espesor X 2 ■ Sean Ti y T 2 las temperaturas a ambos lados 
de la pared. ¿Cual es el flujo de calor a través de la misma?. Evaluar 
la relación entre el flujo de calor si se sustituye un vidrio plano de 
una ventana de 6 mm de espesor, por dos vidrios de 3 mm cada uno, 
separados por 3 mm de aire. En este caso el aire puede considerarse 
como un cuerpo estático suponiendo que el calor que fluye a través del 
mismo obedece a la ley de Fourier. klvidrio] = 75 Wm^K -1 , klaire] = 
2.5 Wm^K” 1 
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*MétocLo: Calcular los flujos en ambos casos y aplicar la ley de las proporcionali¬ 
dades 

• Sea Ti la temperatura en la separación de las dos superficies. El flujo de calor a 
través de la primera pared será 


<Pi=-ki\^ A 

Ax-\ 


( 10 . 1 ) 


mientras que a través de la segunda 


x _ , T 2 - Tí 

<P2 — ^2—T- A 

Axo 


( 10 . 2 ) 


En régimen estacionario ambos flujos son idénticos por lo que 


Ti — Ti 


Axi 

ki 


A 


Ti-T 

A X‘2 
k'2 


A 


(10.3) 


Aplicando la propiedad de las proporciones f = ^ = llegamos a 

, Ti - Ti 


Atci i Ax- 
fci ^ ko 


-A 


(10.4) 


• En el caso de que fueran tres las capas de la pared el problema se generalizaría 
de forma similar dando lugar a 


Ti-Ti 


Ax\/k\ + Axi/ki + Axz/kz 


A 


(10.5) 


que para compararlo con el flujo en el caso de que hubiera sólo un capa pon¬ 
dríamos 


T 2 -T x 

Ax/ki 


A 


( 10 . 6 ) 


La relación entre ambos flujos está dada por 


<j>' _ Axi/ki + A xi/ki + Ax^/kz 
(j> Ax/ki 


(10.7) 


• En el problema que nos ocupa como k\ = y Aaq 
que 


Ax¡ = 1/2 Ax llegamos a 


£ = i + !^ = i. 


-i V -1 


2 fe 2 


75 Wm^K 
2.5 Wm~ 1 K~ 1 


31 


( 10 . 8 ) 


el flujo de calor sería 31 veces mayor en el caso de que la ventana tuviese un 
único cristal. 
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10-2) El crecimiento de la capa de hielo al descender la temperatura es 
un factor relevante en la navegación por los lagos y ríos de las re¬ 
giones nórdicas. Evaluar la variación con el tiempo de la capa de 
hielo en la superficie de un lago en el que la temperatura ambien¬ 
te cae bruscamente de 0 °C a —20°C'. Conductividad térmica del hielo 
k = 2.2 Wmr l K ^ 1 , calor latente de fusión del agua lf = 80 cal <7 1 y den¬ 
sidad del hielo p = 0.971 g cm~ 3 . (Despreciar el efecto de la variación 
de la presión con la profundidad). 


Considerar el flujo de calor cuando se forma una capa infinitesimal de hielo 
e integrar la expresión sobre un intervalo de tiempo finito. 

• Supongamos que sobre el estanque se encuentra una capa de hielo de espesor 
h. Si consideramos que su espesor aumenta una cantidad dh , la masa de agua 
congelada será 

dm = pAdh. (10.9) 

donde p es la densidad del hielo y A la superficie del estanque. 

• El calor cedido por dicha masa al congelarse será 

SQ = Ifdm = IfpAdh (10.10) 

donde lf es el calor latente de fusión. 

• Por tanto el flujo de calor a través de la superficie del estanque deberá ser 

Q = l fP A d ± (10.11) 

• Ahora bien la ley de Fourier nos permite poner 

Q = -k Te ~ Ti A (10.12) 

donde k es la conductividad térmica del hielo y T e y denotan la temperatura en 
el exterior y en la interfase líquido sólido respectivamente. Por tanto igualando 
ambas expresiones llegamos a 

< 1013 > 

Una ecuación que podemos integrar de la forma 

[ h^-dt = - í k Te ~ Ti dt (10.14) 

Jo dt J o Ifp 
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que conduce a 


h(t) 



(10.15) 


es decir que el crecimiento de la capa de hielo es proporcional a la raiz cuadrada 
del tiempo. 


• Numéricamente con los datos aportados por el problema 
k = 2.2 J m~ 1 s~ 1 K ~ 1 

lf = 80 cal g x 4.18 JcaZ -1 x 10 3 <7 kg = 3.344 x 10 5 J fcg -1 
p = 1000 kgm ~ 3 


h(í) = t/2 x 2. 2 x Jm~ 1 s~ 1 K 1 


20 /\ 


3. 344 x 10 5 J kg x 1000fcp?n -3 


í 


que nos ofrece un valor de 


h(t) = ( 


5.13 x 10 _4 m s = 


)V~t 


(10.16) 


El efecto del crecimiento de la capa de hielo en los mares polares sobre la na¬ 
vegación tuvo su punto culminante en la epopeya de la expedición de 1914-1916 
de Ernest Shackleton, quienes en su intento de ser los primeros en cruzar la 
Antártida, vieron como su barco, ”E1 Endurance” era atrapado por el hielo, hun¬ 
diéndose posteriormente por la presión de éste. Shackleton y sus hombres se 
dejaron arrastrar sobre bancos de hielo, y consiguieron llegar en dos botes a la 
Isla Elefante, un lugar inhóspito y deshabitado desde donde Shackleton navegó 
con seis hombres en uno de los botes más de mil quinientos kilómetros hasta la 
isla de San Pedro en las Georgia del Sur, una travesía de 16 días soportando 
un temporal en las embravecidas aguas del Atlántico Sur. Una vez en la isla, 
arribaron por el lado contrario al que se encontraba la estación ballenera, por lo 
que con dos de sus hombres hubo de cruzar las nevadas cumbres de la isla en 
una travesía de 36 horas sin descanso. Finalmente y no sin severos contratiem¬ 
pos Shackleton consiguió rescatar a los 28 componentes de su expedición sanos 
y salvos. Paradójicamente después de haber sobrevivido durante dos años en los 
hielos y mares antárticos muchos de ellos perecieron en las trincheras europeas 
en la Primera Guerra Mundial, iniciada antes de su partida y que aún no había 
terminado a su regreso. Shackleton regresó a la Antártida al mando de otra 
expedición falleciendo de un ataque al corazón durante la misma el 4 de Enero 
de 1922. 
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10-3) Calcular el flujo de calor a través de las paredes de un tubo cilindrico 
de acero de diámetro interior rq = 2 cm y exterior r 2 = 2.5 cm cuyo 
interior se mantiene a una temperatura Ti = 100°C' y su exterior a 
una temperatura T -2 = 25 °C. (Conductividad térmica del acero k = 

16.3 Hm-'r 1 ) 


*Método: Calcular el flujo integrando la ley de Fourier y considerando que el sis¬ 
tema tiene simetría radial y que el flujo a través de cualquier superficie concéntrica 
con la superficie interior es constante. 


• Sea dA un elemento infinitesimal de superficie en el interior de las paredes del 
cilindro a un radio r perpendicular a éste. El flujo de calor a través de dicha 
superficie está dado por 

d$ = -k^dA (10.17) 

Or 


• Si consideramos un cilindro a un radio r de longitud l, el flujo de calor a través 
de dicho cilindro, supuesta la temperatura interior y exterior uniforme a lo largo 
de l está dado por 

dT 

$ = —k-^2nrl (10.18) 

Or 

La ecuación anterior podemos escribirla de la forma 

1 dT 

-$ = -k2nl— 10.19) 

r or 


que podemos integrar entre los límites r\ y r 2 radio interior y exterior respecti¬ 
vamente. 



-<f>dr = — 
r 



dT 

k2nl--dr 

Or 


( 10 . 20 ) 


• Teniendo en cuenta que en régimen estacionario el flujo a través del cilindro es 
independiente de r, la expresión anterior queda de la forma 


de donde 



= —/c27rZ (T> — Ti) 


$ = 


—k2nrl (T 2 — T \) 



( 10 . 21 ) 

( 10 . 22 ) 


• E 11 nuestro caso T 2 — T\ = 75 K; r 2 = 2.5 x 10 2 m; r\ — 2 x 10 2 m; k = 
16.3 por lo que sustituyendo numéricamente llegamos a 


$ 

T 


= -8.3 x 10 5 Wm~ l 


(10.23) 


es decir, se pierden 8.3 x 10 5 W por cada metro de tubería. 
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10-4) Supongamos un cilindro de radio interior rq y radio exterior r 2 de 
manera la temperatura en el interior es T\ y el exterior del cilindro 
está en contacto con un fluido a temperatura T a , y sea h el coeficiente 
de convección entre el cilindro y el fluido. Calcular el flujo de calor. 


*Método: Calcular el flujo de calor del cilindro hasta la atmósfera utilizando la ley 
del enfriamiento de Newton y el flujo de calor desde el interior al exterior. Ambos 
flujos son iguales 


• En este caso en régimen estacionario, el flujo de calor que atraviesa la superficie 
exterior es el mismo que el que es intercambiado con el fluido. 

• Según la ley del enfriamiento de Newton este flujo podemos calcularlo como 

= h(T 2 - T a )2nr 2 l (10.24) 

donde el término 2nr 2 l representa la superficie de contacto entre el fluido y el 
cilindro. 


• Por el resultado del problema anterior tenemos pues 

—k2irl (' T 2 — Ti) 


$ = 




= h(T 2 - T a )2nr 2 l 


que podemos ponerlo de la forma 

-2nl (T 2 - Ti) _ (T 2 - T a )2nl 


$ = 


i 

hr 2 


Ms) 

que aplicando la ley de las proporciones (f = f = f fff) conduce a 

2tvI (Ti - T a ) 


$ = 


Í ln (n ) + 


(10.25) 


(10.26) 


(10.27) 


que sólo depende de la temperatura en el interior Ti, y la del fluido exterior T a 
pero no de la temperatura en la superficie exterior T 2 ■ 


• Si suponemos fijo el radio interior, así como las temperaturas exterior e interior, 
la dependencia del flujo con el radio exterior r 2 está dada por una función no 
lineal de la que podemos calcular el valor mínimo evaluando 


_d_ d 2nl (Ti - T a ) Ti-T a (\ _1_ 

in dr2 i !■■(?) + * (í>"(s)+ ¿) 2 ' tr2 hTl 
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• Dicho valor mínimo corresponde pues a 



= 0 


(10.28) 




(10.29) 


es decir el flujo es mínimo cuando el radio exterior es igual al cociente entre la 
conductividad térmica y el coeficiente de convección. 

• Suele ocurrir que por el interior del cilindro circula asimismo un fluido a tem¬ 
peratura Ti y denotemos entonces por hi al coeficiente de convección con dicho 
fluido y por I 12 al coeficiente de convección con el fluido externo. 

• En este caso el flujo de calor entre el fluido interno y el cilindro viene determinado 


por 


$¿ = -fti(Ti - Ti)2-Knl 


(10.30) 


donde el signo negativo establece el convenio de que el flujo es positivo si tiene 
lugar hacia afuera, en la dirección radial del cilindro. 

• E 11 este caso, en régimen estacionario, el flujo del fluido interior al cilindro, a 
través del cilindro y del cilindro al exterior son el mismo, de manera que podemos 
poner 



(10.31) 


de manera que aplicando la ley de las proporciones llegamos a: 

» _ 27tZ (Tí — T a ) 



(10.32) 






10-5) Determinar la temperatura de un satélite artificial en la órbita te¬ 
rrestre suponiendo que se comporta como un cuerpo negro, y consi¬ 
derando el sol como un cuerpo negro esférico de radio R = 6.96 x 10 8 m 
cuya temperatura en la superficie es de T a ss 6000 K [radio medio de la 
órbita terrestre d = 1.49 x 10 n m] 


*Método: Considerar que la energía emitida por el sol se distribuye en una 
esfera de radio la órbita terrestre y evaluar la fracción incidente sobre la 
superficie de un objeto situado sobre la misma. 

• La energía emitida por el sol por unidad de tiempo y superficie es: 

Tí = ctTq (10.33) 

por lo que la energía total emitida por el sol por unidad de tiempo en todas 
direcciones es 

£ soi = 4uR 2 aT* (10.34) 

• La energía emitida por el satélite que se encuentra en régimen estacionario a una 
temperatura T e estará dada por 

£ S at = 4nr 2 aT F (10.35) 


• La fracción de energía emitida por el sol que es absorbida por el satélite ( Tí sa t ) 
está dada por la fracción de superficie que ocupa el satélite en una esfera de 
radio la distancia del satélite al sol 

•2 


Tísat = 47 rR 2 aT n 


7r r 
ótt d? 


(10.36) 


• En estado estacionario la energía absorbida por el satélite es igual a la emitida 
{£sat = Tísat ) de forma que 


4tt r 2 oTj = 4nR 2 aTf 


7rr 

And 2 


T -= T - ,| S 


que numéricamente nos ofrece: 

T f = 6000 K i 


6.96 x 10 8 m 
2 x 1.49 x 10 n m 


289 1< 


(10.37) 

(10.38) 

(10.39) 


• Ejercicio adicional: 

— Calcular la constante solar, es decir la energía proveniente del sol que incide 
sobre la tierra por unidad de tiempo y superficie. 
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10-6) Un componente de un criostato está formado por un recipiente cilindrico 
que contiene helio líquido en su punto de ebullición normal a 4.2 A'. Ro¬ 
deando al mismo se encuentra otro recipiente cuyas paredes se man¬ 
tienen a la temperatura de ebullición normal del nitrógeno (77.3 K ) 
habiéndose hecho vacío en el espacio intermedio. Suponiendo que 
la superficie de ambos cilindros se comporta como un cuerpo negro, 
determinar las pérdidas de helio por unidad de tiempo por evapora¬ 
ción debido al suministro de calor por radiación. [Altura del crios- 
tato h=10 cm, diámetro d=6 crn , calor latente de vaporización del helio 
l v = 5.0 cal/g\ 


*Método: Evaluar el flujo de calor utilizando la ley de Stefan-Boltzmann 


• El calor absorbido por unidad de tiempo por el helio por radiación de las paredes 
exteriores Q es igual a: 

Q = HTk - T He)A (10.40) 

siendo A la superficie y Tjv 2 y T¡j e las temperaturas respectivas del nitrógeno y 
helio en ebullición. 


• Por tanto la pérdida de masa del helio por evaporación por unidad de tiempo rh 
será igual a 


. Q 

m= T 

ly 

donde l v es el calor latente de vaporización del helio. 
• La superficie del recipiente está dada por 


A = 7T hd + 2 X 7T ( ^ 


donde d es el diámetro del cilindro y h es su altura. 
• Las pérdidas de helio pueden evaluarse pues como 


lv 


-A 


(10.41) 


(10.42) 


(10.43) 


Con los datos del problema 


A = 3.1416 x 10 x 10~ 2 x 6 x 10~ 2 + 2 x 3.1416 


(sxu^y 


= 2.45 x 10“ 2 m 2 
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a = 5.67 x 10 8 Js 1 m 2 A' 4 
l v = 5.0 cal g~ l x 4.18J cal~ l — 20.9 J 
- T Áe = 77 4 - 4.2 4 = 3.5 x 10 7 A' 4 

• Sustituyendo numéricamente 


5.67 x 10 -8 J s -1 m" 2 A' -4 x 3. 5 x 10 7 A' 4 


m = 


20.9 Jg 


-i 


2.45x10” 2 to 2 = 2. 33xl0 _3 gs” 1 


es decir se pierden rh = 8.4 gramos de helio por cada hora. 



La determinación numérica de la conducción del calor está basada en los trabajos del francés Joseph Fourier (1768- 
1830) quien publicó en 1822 su teoría analítica del calor. Las herramientas matemáticas utilizadas han dado lugar 
a una rama de las matemáticas de amplia aplicación. Algunos autores citan también a Fourier como el primero en 
intuir en un ensayo publicado en 1827 el posible calentamiento de la atmósfera debido a la presencia de gases en 
la misma, un efecto posteriormente conocido como efecto invernadero y que es la base sobre las teorías actuales 
sobre el cambio climático. 
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10-7) El gradiente geotérmico de la tierra, es decir la variación de la tempe¬ 
ratura con la profundidad en las proximidades de la superficie es de 1 K 
cada 33 m. Suponiendo que dicho gradiente es consecuencia de la gene¬ 
ración de calor de forma uniforme en el interior de la tierra, calcular 
el calor generado por unidad de masa y unidad de tiempo asumien¬ 
do que la tierra es una esfera homogénea de densidad p = 5.5 gcm~ 3 
y conductividad térmica k = 2.5 Wm~ 1 K~ 1 con un radio r = 6378 km. 
¿Cuál sería la temperatura en la superficie suponiendo que no exis¬ 
tiera aportación externa de energía proveniente del sol?¿Cuál sería la 
temperatura a 10 km de profundidad? 


*Método: Considerar el flujo hacia el exterior considerando el término 
del calor generado en cada elemento de volumen e integrar la ecuación 
resultante en coordenadas esféricas. Para evaluar la temperatura en la 
superficie considerar que todo el calor generado es emitido en forma de 
radiación y relacionarlo con la temperatura en la superficie mediante la ley 
de Stefan Boltmann. Para calcular la temperatura en el interior integrar 
la expresión del gradiente geotérmico. 

• Para calcular la temperatura en el centro de la esfera y por tanto el gradiente 
geotérmico tenemos en cuenta que la disminución de la energía por unidad de 
tiempo en cualquier volumen V interior a la tierra es igual al flujo de energía que 
sale por la superficie del mismo. Por una parte siguiendo el desarrollo teórico 
utilizado en la deducción de la ecuación de Fourier este flujo está dado por la 
expresión 

<f> = - í kV 2 TdV (10.44) 

Jv 

donde V 2 T = Por otra parte dicho flujo de calor corresponde 

a la disminución de energía consecuencia de la variación de temperatura de cada 
elemento de volumen más el calor generado en el mismo. En estado estacio¬ 
nario la variación de temperatura con el tiempo es cero y esta disminución de 
energía corresponde al calor generado que podemos ponerlo como el calor gene¬ 
rado en cada elemento infinitesimal de volumen integrado sobre todo el volumen 
considerado. 

<j> = í pqdV (10.45) 

Jv 

Por q hemos denotado el calor generado en el interior de la tierra por unidad 
de masa y unidad de tiempo, pdV = dm. corresponde a la masa de un elemento 
infinitesimal de volumen dV, siendo p la densidad. 

• Igualando ambas expresiones del flujo llegamos a 

[ pqdV = f kV 2 TdV (10.46) 

Jv Jv 
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Como esta relación debe verificarse para cualquier volumen V se tiene que los 
integrandos son iguales y por tanto 

pq = -kV 2 T (10.47) 


• Como el problema tiene simetría esférica y T sólo depende de r podemos poner 
utilizando coordenadas esféricas 


que conduce a 


1_d_ 

r 2 dr 


que podemos integrar escribiendo 


[ r d_ 

Jo dr 


lo que lleva a 


2 dT 
dr 

dT 

dr 


1 d i 

r 2 dr ' 


(10.48) 

:-£) 

II 

(10.49) 

- ) dr = 

J 

= - í T~r 2 dr 

Jo k 

(10.50) 

1 

3 

PQ 
k r 

(10.51) 


que representa el gradiente geotérmico o variación de la temperatura con la 
profundidad en el interior déla tierra, de donde deducimos que el calor generado 
en el interior de la misma por unidad de masa está dado por 


3 k dT 
rp dr 


(10.52) 


• Para calcular la temperatura en la superficie hemos de suponer que en estado 
estacionario el calor generado en el interior de la tierra por unidad de tiempo 
es igual al perdido por ésta desde la superficie por radiación. El calor generado 
puede ponerse como el calor generado en cada elemento infinitesimal de masa 
integrado sobre toda la esfera, y el calor perdido por radiación como la suma del 
calor perdido en cada elemento de superficie, es decir 


í pqdV 
Jv 


I 

Jn 


aTfdfl 


(10.53) 


donde la primera integral corresponde a una integral de volumen sobre toda la 
esfera y la segunda a una integral de superficie sobre ésta y donde T s denota la 
temperatura en su superficie. 


La expresión (10.53) conduce a 


pq f dV = crT s 4 f dfl 

Jv " Jn 


que nos lleva a 


pq-nr 3 = aT^Anr 2 

ó 


(10.54) 

(10.55) 
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de donde deducimos que la temperatura en la superficie está dada por 


1 PÚr\ 

3a) 


1/4 


(10.56) 


• La temperatura en función de la profundidad podemos calcularla integrando 
(8.63) desde la superficie r hasta una profundidad z, n = r — z 


Jr dv ,/r 


-—rdr 
3 k 


lo que nos lleva a 


T(ri) = T(r ) - ~ ((r - z) 2 - r 2 ) 
• Numéricamente en nuestro caso 


(10.57) 


(10.58) 


p = 5.5 g cm 3 x 10 3 kgg 1 x 10 6 cm 3 m 3 = 5. 5 x 10 3 kgm 3 
k = 2.5 W m~ 1 K~ 1 
-0.03 K m _1 


dT 

dr 


• El calor generado internamente (10.52) 
3 x 2.5 VEm“ 1 A'~ 1 


q = ~: 


(—0.03A' m” 1 ) = 6.4141 x 10 ~ 12 W kg~ 


6378 x 10 3 m x 5.5 x 10 3 kgm~ 3 
• Temperatura en la superficie asumiendo única fuente el calor interno (10.56) 


T=, 3 


1 5.5 x 10 3 kgm~ 3 x 6.4141 x 10 ~ l2 Wkg~ x x 6378 x 10 3 m\ 1/4 


5.67 x 10~ 8 J s~ 1 m~ 2 K~' í 


= 34 A' 


Temperatura a 10 km de profundidad (10.58) 

1 5. 5 x I0 3 kgm~ 3 x 6.4141 x 10 -12 VE kg~ l 

T(0) = 34A H- - ------— x 

6 2.5VTm~ 1 A 1 


x ((6378 x 10 3 ) 2 - (6368 x 10 3 ) 2 ) = 333 A' 
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